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news and views 


PGAP SurRvVEY 


The membership was surveyed by a post card ques- 
tionnaire to obtain better information on needs and in- 
terests. The number of returns was gratifyingly high, 
but, unfortunately, a good many of the group did not 

complete the questionnaire in the desired form so that 
_ they could not be tallied. A sufficient number of members 
did respond correctly and the tally of the correctly 
filled-in questionnaires probably gives a statistically 
_ satisfactory result. 

The post card questionnaire is reproduced here for 

convenience in checking the tabulated returns: 


PGAP Survey 
(Number Replies in Order of Importance) 


Field of Interest: 
Antennas 
Propagation 
Microwave Tech. 
Microwave Optics 
Radio Astronomy 


Paper Pref. for Transactions: 
Experimental Type 
Analytical Type 
Applications Type 


Meeting Preferences: 

{ More-Less National Meetings 

_ Prof. Group Membership: More-Less Local Chapter Mtgs. 
Antennas & Propagation 

_ Microwave Theory & Tech. 
List Other Prof. Groups 

~ Comm. Systems 
Telemetry & Remote Control 


Local Chapter Membership: 

Do you support your Local 
PGAP Chapter? Name Local. 
Do you support any Local IRE 
Group? 


The results of the survey are shown in Table I. Under 
- Field of Interest, a tally for first through fifth place is 
shown for each field. Thus under Radio Astronomy, five 
people indicated this field as their first preference, thirty- 
one people indicated this to be their second preference, 
and so forth. Ideally, if everyone filled in all numbers 
from one to five on the questionnaire, the vertical and 
horizontal totals should be the same. This is not so be- 
cause of the interesting display of individuality in the 
replies, where some filled in perhaps two out of the five 
boxes, etc. 

However, it does appear that Antennas, Propagation, 
and Microwave Techniques, in order, represent the 
greatest interest to our group with an interesting de- 
emphasis on Radio Astronomy. 


TABLE I 
PGAP SuRVEY 


Field of Interest 


Antennas 249 | 161 42 8 2 
Propagation SEN) TISKe 98 68 22 
Microwave Tech. 84} 98 | 90] 44 24 
Microwave Optics 18 50 | 104 90 56 
Radio Astronomy 5 31 60 | 67 130 
Paper Preference 
Experimental 113 | 190 | 59 
Analytical 178 | 90 | 164 
Applications 155 | 153 | 148 
Meeting Preference 
More National Meetings 157 
Less National Meetings 94 
More Local Chapter Meetings | 189 
Less Local Chapter Meetings 42 


Under Papers Preference, there is no pronounced 
trend. While analytical papers appear to lead as first 
preference, about as many people also put such papers 
in third place. One can at least say that people are not 
neutral about analytical type papers; they are either for 
them or against them. 

Finally, there is a fairly definite preference for more 
national meetings and more local chapter meetings. 

In a general way, the results are in line with the 
thinking of the Administrative Committee as to their 
preconceived feelings of the interests and needs of the 
membership. The Administrative Committee is thus en- 
couraged to proceed in its previously set general direc- 
tions, particularly insofar as more national and local 
chapter meetings are concerned. 


LocAL CHAPTERS 


An application for a new section of PGAP at San 
Diego has been approved by the Administrative Com- 
mittee. On approval of the IRE Executive Committee 
this group can officially get under way. 

Formation of a PGAP Chapter in Columbus, Ohio is 
under consideration. Dr. D. D. Knig of Johns Hopkins 
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University is the member of the National Administra- 
tive Committee responsible for the liaison and promo- 
tion of local chapters and he will be glad to assist in any 
way possible. 


Washington Chapter—The Washington Chapter looks 
as though it will become quite active. A regular meeting 
was held on March 6, 1956 at George Washington Uni- 
versity with speaker Dr. D. D. King of Johns Hopkins 
who presented a paper on recent work at his laboratory 
on millimeter wave research. The agenda for the meet- 
ing appears quite interesting. Unfortunately, a report 
of the discussion is not available at this writing. The 
agenda, however, is reproduced below: 


Business 


1) A nominating committee will be appointed to con- 
sider candidates for chapter offices next year. 

2) A symposium is being planned for this Fall and 
would possibly be organized along the same lines as the 
Scatter Symposium which was held last November. 
Some of the groundwork for this has been started and 
with the cooperation of chapter membership there is no 
reason why this area cannot be the site for another out- 
standing symposium. 

3) There are a number of items relating to the over- 
all policies of PGAP, particularly with regard to meet- 
ings and publications, which should be of interest to all 
members. This might well be the source of a lively and 
beneficial discussion. 


A Reminder 


Our meeting this Tuesday night is next to the last 
which will be held this year. There is much to be done 
toward outfitting the chapter for its intended purpose, 
which is, primarily, to serve the needs of the member- 
ship in advancing our technical knowledge in the field 
of antennas and propagation. 

Problems always exist relative to the procurement of 
interesting speakers, balancing of the different interests 
among members, timing of the meetings, and coopera- 
tion with other PG chapters. I most sincerely urge that 
you make it a definite point to attend this next meeting 
and express your views so that sound policies can be 
formulated for direction of our efforts next year. 

—Coleman Goatley, Chairman 


MEETINGS 
WESCON, August 21-24, 1956 


The 1956 WESCON will be held in Los Angeles, Cali- 
fornia this year and as in the past the PGAP will spon- 
sor the technical sessions on antennas and propagation. 

Paper summaries are required by April 15, 1956 and 
should be addressed to: 

Di Res. Ehiott 
Hughes Research Laboratories 
Culver City, California 
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Authors will be notified of the disposition of their pa- 
pers by May 15. 

Since a Convention Record will be published along 
lines parallel to the [RE CONVENTION RECORD, it will be 
necessary to have reproducible copies of all papers ac- 
cepted in Dr. Elliott’s hands by June 15, 1956. 

The recommended format for the reproducible copy 
will be furnished at the same time the disposition of the 
papers is announced. 


WESCON Schedule for Papers 


April 15—Summaries to Dr. R. S. Elliott. 
May 15—Author notification. 
June 15—Reproducible copies to Dr. Elliott. 


CORRESPONDENCE 


The following is an excerpt from a letter from John 
Smyth to George Sinclair containing some thoughts 
worthy of discussion: 

“At no time has a reviewer been identified to an au- 
thor. When the paper is acceptable, with or without 
some revision, the reviewer’s remarks are sent along to 
the author as a service to him. When the paper appears 
in print the reviewer as well as anyone else can send 
along detailed comments for publication. The reviewer’s 
comments are the basis for accepting or rejecting papers 
for publication; they are also an aid to the author in 
that they frequently afford him an opportunity to im- 
prove on his paper, correcting errors, filling in desired 
details and the like. Rumsey’s statement, ‘This ma- 
terial may be appended to the published paper (as com- 
ments by the referee) if the author so desires’ is quite 
acceptable. Or, in the event the author wishes to ignore 
some of the referee’s comments, the reviewer can always 
write a letter to the editor about the paper once it is in 
print. 

I think it would improve our TRANSACTIONS if more 
people sent along their remarks about published papers 
for the Communications section. One might go so far as 
to say the English system of simultaneously publishing 
the comments of the referee and rebuttal by the author 
is better than the way we handle our reviews. In fact, I 
personally like it. 

Rumsey’s idea of handling reviews is certainly worthy 
of considerable thought and discussion and I recom- 
mend that it be placed high on the agenda of the next 
Administrative Committee meeting.” 


CALENDAR OF EVENTS 


Apr. 30-May 3, 1956: URSI Spring Meeting, National 
Bureau of Standards, Washington, D.C. 


May 1-3, 1956: IRE-RETMA-AIEE-WCEMA Electron- 
ic Components Symposium, U. S. Department of In- 
terior, Washington, D. C. 


May 14-16, 1956: National Aeronautical and Naviga- 
tional Conference, Hotel Biltmore, Dayton, Ohio. 


1956 


May 21-22, 1956: Symposium on Reliable Applications 
of Electron Tubes, University of Pennsylvania, Phila- 
delphia, Pa. 


June 4-6, 1956: Second Annual Radome Symposium, 
Ohio State University, Columbus, Ohio. 


Aug. 20-21, 1956: National Telemetering Conference, 
Biltmore Hotel, Los Angeles, Calif. 


Aug. 21-24, 1956: IRE-West Coast Electronic Manu- 
facturers’ Association, WESCON, Biltmore Hotel, 
Los Angeles, Calif. 


Sept. 11-12, 1956: Second RETMA Conference on Re- 
liable Electrical Connections, University of Pennsyl- 
vania, Philadelphia, Pa. 


News and Views 
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Sept 17-21, 1956: Instrument-Automation Conference 
and Exhibit, Coliseum, New York, N. Y. 


Sept. 24-25, 1956: Industrial Electronics Symposium, 
Manger Hotel, Cleveland, Ohio. 


Oct. 1-3, 1956: National Electronics Conference, Chi- 
cago, Ill. 


Oct. 1-3, 1956: Canadian IRE Convention and Exposi- 
tion, Automotive Building, Exhibition Park, Toronto, 
Canada. 


Oct. 8-9, 1956: Second Annual Symposium on Aero- 
nautical Communications, Utica, N. Y. 


Oct. 10-12, 1956: Symposium on Applications of Optical 
Principles to Microwaves, Washington, D. C. 


OSS Ce 
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Variational Principles for Electromagnetic 
Resonators and Waveguides* 
A. D. BERK 


Summary—Variational expressions are presented for the propa- 
gation constants of a waveguide and for the resonant frequencies of a 
cavity directly in terms of the field vectors for situations not covered 
in the literature. These situations occur when the electromagnetic 
problem is not susceptible of a scalar formulation and are typified 
by the presence of inhomogeneous or anisotropic matter. The varia- 
tional expressions are applied to several illustrative examples and 
the results are compared to known exact solutions. The variational 
expressions are also applied to the derivation of certain perturbation 
formulas, some of which were previously obtained by different 
methods. 


INTRODUCTION 
NOWN variational expressions for the resonance 


frequencies of a resonator are restricted to the 

special case when the electromagnetic field is 
derivable from a single scalar function which satisfies 
the Helmholtz equation.! When the substance within a 
cavity is inhomogeneous or anisotropic, such a scalar 
formulation is inadequate. The need for vector varia- 
tional principles is thus apparent. It is the purpose of 
the first part of this paper to present such variational 
formulas for resonant frequencies directly in terms of 
the field vectors. In problems of propagation through 
anisotropic or inhomogeneous media, variational ex- 
pressions for the propagation constant are also of inter- 


* Manuscript received by the PGAP, August 11, 1955; revised 
manuscript received January 27, 1956. The research reported in this 
paper is based on a chapter of a doctoral thesis submitted by the 
author to the Massachusetts Institute of Technology in September, 
1954 and was supported in part by the Signal Corps, the Office of 
Scientific Research, Air Research and Development Command, and 
the Office of Naval Research. 

t pee Aircraft Co., Culver City, Calif. 

M. Morse, and ‘i. Feshbach, “Methods of Theoretical 
ey ” McGraw-Hill Book Co., New York, Ne Ye; pe ti 12 19539 


est. The second part of this paper is concerned with the 
derivation of such expressions. The variational expres- 
sions for resonant frequencies and those for propagation 
constants appearing in this paper are to be distinguished 
from variational expressions for reflection coefficients, 
scattering amplitudes or impedance matrices derived by 
other authors.!~* In the third part, the variational for- 
mulas are applied to illustrative specific examples and 
to the derivation of some perturbation expressions. 


VARIATIONAL PRINCIPLES FOR 
RESONANT FREQUENCIES 


E-Field Formulation 


Consider a resonator with perfectly conducting walls 
enclosing a medium of permittivity € and permeability 
uw. Both w and e may be tensors‘ and functions of po- 
sition. Let a resonant angular frequency be w and let the 
corresponding electromagnetic field be characterized by 
the vectors E and H. The following is then asserted to 
be a variational expression for w, provided € and yw are 
hermitian; 7.e., provided no losses are present: 


i (curl E*) -y-!- (curl E)dv 
Oe a (1) 


| Br Bde 


2 D. S. Saxon, Notes on Lectures by J. Schwinger, “Discontinui- 
ties in Waveguides, ” Radiation Laboratory, M.I.T.; February, 1945. 

3H. Levine and J. Schwinger, “On the theory of ‘diffraction by an 
aperture in an infinite screen plane,” Phys. Rev., vol. 74, pp. 958-974; 
October, 1948. 

4 For simplicity, no special symbols will be used to indicate the 
tensor character of these quantities. 


1956 


The integrals are over the volume of the resonator, y-! 
is the inverse of u, E* is the complex conjugate of E. To 
prove this assertion, we must show that those field con- 
figurations E and E* which render w stationary are 
solutions of 


curl (u-1-curl LE) — we- E = 08 (2) 


and of the complex conjugate of (2) and have vanishing 
tangential components at the boundary. This is indeed 
the case, for, on varying E and &* in (1) we obtain, after 
utilizing the hermitian character of € and yu, the follow- 
ing expression for the variation of w: 


2) 


(69) = 
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i} (curl E*)-y-!- (curl E)dv -$ n-|E X (u-*-curl E*) |ds — § n- |E* X (u-?-curl E) |ds 
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defined by (2) and the boundary condition »XE=0. 
Admissible trial fields must have vanishing tangential 
components at the boundary, must be continuous to- 
gether with their first derivatives and must possess finite 
second derivatives everywhere in the cavity except at 
surfaces where e€ and wp are discontinuous. At such sur- 
faces n XE and »X(u-!-curl E) must be continuous. 
Eq. (1) can be modified so that trial vectors E are not 
required to satisfy the boundary condition nXE=0 at 
the wall of the resonator. This can be achieved following 
a known general method® by adding appropriate terms 
to the numerator of (1) as follows: 


rv 


, (S) 


[ eka 


[ Be cun (u-!- curl £)dv — gn [EX (uw *-curl E*) |ds 


or, on combining the first and third terms in the nu- 
merator, 


(6) 


[ Bene 


(f E* «: Bdo) so? 


= fom. [curl (u-!-curl E) — we E]dv 


— f 5B (n Cee curl eds 


+ f dE: [curl (u-*- curl E*) — we*- E* |dv 


= § En Monee wcurl ds. (3) 


The second and fourth integrals are over the boundary 
of the cavity. Their appearance is due to the utilization 
of the vector identity 


f Accu Bdv = [ Becurt ade $$ wB x A)ds, (A) 


where z is the outward normal unit vector. The varia- 
tion of w? will vanish, provided E satisfies (2), E* satis- 
fies the complex conjugate of (2) and provided the 
surface integrals in (3) vanish. The latter condition can 
be satisfied only if »xX6£ and its complex conjugate 
vanish over the boundary, since ” X(u~!-curl £), being 
proportional to the tangential component of the mag- 
netic field, cannot vanish over the complete boundary. 
Eq. (1) is thus a variational formulation of the problem 


5 This equation is the result of eliminating the magnetic field 
between Maxwell’s equations. 


When the distribution of matter within the cavity is 
discontinuous, (1) can be further modified so that trial 
fields will zot be required to have continuous tangential 
components of — and of (u7!-curl £). The modification 
consists of adding to the numerator of (1) the term 


— i) n- | E.* & (uy}-curl Ey) — B_* X (u-?-curl E_) |ds 


— complex conjugate, 


where the subscripts +, — refer to values on opposite 
sides of the surface of discontinuity and the integrals are 
over such a surface. 

The passage from (1) to (6), enables one to expand 
the class of admissible trial functions. The value of such 
a procedure will be demonstrated in the section on Ap- 
plications, where trial functions outside of the class as- 
sociated with (1) will be utilized. 

Finally, (1) can be modified to apply when the bound- 
ary condition at the walls is 


nm X (uw-!-curl E) = — jwoY- Ey (7) 


where (jwY) is a hermitian two dimensional dyadic 
and F;, is the tangential component of /. The term 
to be added to the numerator of (1) is, in this case, 
f E.*: (jwY)-E.ds. A practical example where (7) applies, 
is at the junction of a cavity and a tuning stub. 


6 R. Courant and D. Hilbert, “Methods of Mathematical Physics,” 
First English Translation, Interscience Publishers, pp. 208-211, 1953. 
See also reference 1, pp. 1131-33. 
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H-Field Formulation 


A variational expression similar to (1) but in terms 
of the magnetic field can be obtained by interchanging 
e and uw and by replacing E by H: 


if (curl H*)-e—!- (curl H)dv 


[tan 
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but they are not required to satisfy any particular 
boundary condition. If matter is discontinuously dis- 
tributed within the cavity, then trial vectors E and H 
must have continuous tangential components at the 
surfaces of discontinuity. The latter restriction and the 
restriction of vanishing trial tangential E at the walls of 
the cavity can be eliminated by the addition of appropri- 
ate terms to the numerator of (9). For example, in the 
following variational expression 


Oe 


3 (10) 


[ Pre nant f He Has 


The distinguishing feature of (8) is that trial vectors H both trial vectors F and // are unrestricted at the wall. 


are not required to satisfy the proper boundary condi- 
tion, m X(e~!-curl H) =0, at the wall. However, the dif- 
ferentiability and continuity conditions are the same as 
those for E in connection with (1). 

Eqs. (1) and (8) reduce properly to the variational 
principle of reference 1 in the special case when the 
electromagnetic problem can be expressed in terms of a 
single scalar field. 


Mixed-Field Formulation 


In contrast to the preceding formulas which are in 
terms of either the electric vector or the magnetic vec- 
tor, the following variational expression is in terms of 
both field vectors. 


[econ Rav = [ B*-cur Hdv 
w=] ) (9) 
[ Brena + [a -Hde 


ay 


VARIATIONAL FORMULAS FOR THE 
PROPAGATION CONSTANT 


Mixed-Field Formulation 


Consider a waveguide with perfectly conducting walls, 
possibly enclosing anisotropic matter whose distribuion, 
may be a function of the transverse coordinates but not 
of the coordinate along the direction of propagation. 
If z is this coordinate, the field vectors may be expressed 
as E(x, y)e-i7, 3C(x, ye’, where y is the propagation 
constant. E and 3 are three dimensional vectors de- 
pending on x and y only and they satisfy the following 
relations obtained by substituting the field vectors in 
Maxwell’s equations: 


curl € + jou-5C = yja, X E, 
curl 50 — jwe-E = yja. X 3, 


(11) 
(12) 


l| 


where a, is the unit vector in the z-direction. 

Premultiplying (11) by 3¢*, (12) by &*, integrating 
over the cross section of the waveguide and subtracting 
we obtain the following variational expression: 


wf Ere bis +e fae u-seds +7 f &*-cun Keds = jf se*- curt Eds 


[ s0-*a. x Eds — | €*-a, x aeds 


where e¢, are again assumed to be hermitian. That (9) 
is indeed a variational expression can readily be verified 
by computing the variation of w and observing that the 
latter vanishes, provided E, H satisfy Maxwell’s equa- 
tions, provided E*, H™ satisfy the complex conjugate of 
Maxwell’s equations, and provided the trial electric 
vectors have vanishing tangential components at the 
boundary. Thus, admissible trial E-vectors must be 
continuous, must possess first derivatives and their 
tangential components at the boundary must vanish. 
Admissible trial H-vectors are subject to the same con- 
tinuity and differentiability conditions as the E-vectors, 


(13) 


That (13) is indeed a variational expression can be 
shown by evaluating the variation of y and observing 
that the latter vanishes provided €, 3¢ satisfy (11) and 
(12) and provided the tangential component of € van- 
ishes at the walls of the waveguide. Thus, trial fields 
€, 5 must be continuous, and must possess first deriva- 
tives throughout the cavity. At the boundary, the 
tangential component of € must vanish, but 3 is arbi- 
trary. When discontinuities are present in the distribu- 
tion of matter within the cavity, the tangential com- 
ponents of both € and 3¢ must be continuous at the 
surfaces of discontinuity. 
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In the following modified form of (13) both E and # 
are arbitrary at the boundary. 
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the propagation constant of a waveguide. We shall now 
illustrate them by a number of examples. For a particu- 


wf Ete€ds te fat n-seas +7 [ EF cur seds — jf secur Eds — ih n(E x seh)al 


y 


(14) 


J seta. x Eas — f Oa. x sed 


The last integral in the numerator of (14) is over the 
boundary of the cross section of the waveguide. 


€-Field Formulation 


If 3 is eliminated between (11), (12), the result pre- 
multiplied by €*- and integrated over the cross section, 
the following is obtained: 


vf XE) 1 (0.x Eds 
—in f [(curl E*) -w- (a,x E) — (a, X E*) -p - (curl €) |ds 


+ [ (cur E*) -u-)- (curl Eds—u* [ E*« Eds=0, (45) 


Varying €, E€* and y in (15) it can be shown that y¥ re- 
mains stationary provided €, E€* satisfy 


2a X (wu? aX €) +7y [curl (u-a.& E) +a, X (u--curl €) | 
—curl (u7!-curl €)+we-E=0 (16) 


and the conjugate of (16). Eq. (16) is the result of elimi- 
nating 3C between (11) and (12). Hence, it is the proper 
“Euler” equation. Solving (15) for y we obtain a varia- 
tional expression for the propagation constant. 

Admissible trial vectors E must have vanishing tan- 
gential components at the walls, must be continuous 
together with their first derivatives and must possess 
second derivatives except at surfaces of discontinuity of 
the medium enclosed within the waveguide. At such 
surfaces, © and (u-!-curl E—jyy--a,X€) must have 
continuous tangential components. 


5C-Field Formulation 


If in (15) w and e¢ are interchanged and € is replaced 
by 3, a new variational equation results. Admissible 
trial vectors 3C must satisfy the same continuity and 
differentiability conditions as does € in connection with 
(15), but they are arbitrary at the boundary of the 
waveguide. 

The variational principles referred to in the preceding 
sections can be modified so that trial functions will be 
less sensitive to boundary conditions. The method is 
similar to that employed in deriving (5) and (10) and 
will not be repeated here. 


APPLICATIONS 


We have, thus far, offered several variational expres- 
sions for both the resonant frequencies of a cavity and 


lar problem, the choice of a certain variational formula 
in preference to others largely depends on physical con- 
siderations. For example, if the configuration of the 
electric field can be “guessed” more readily than that 
of the magnetic field, then it is sensible to use a formula 
in terms of the electric vector only. 

The successful choice of trial fields also depends on 
our familiarity with the physical aspects of the problem 
at hand. It is important to remember that the conditions 
of admissibility of trial functions are minimum require- 
ments and that one should make every effort to select 
trial fields which satisfy as many of the known features 
of the solution as possible. In particular, one should at- 
tempt to devise trial fields which, at surfaces of discon- 
tinuity, not only have continuous tangential compo- 
nents, but continuous normal components of (e€:£) and 
(u-H) as well. This is because the last set of boundary 
conditions does not follow from the first set unless the 
trial fields satisfy Maxwell’s equations. 

The first three of the following examples are intended 
to illustrate the advantages of approximate calculations 
based on the variational expressions of this paper over 
exact treatments. In the remaining examples, the varia- 
tional expressions are applied to the derivation of cer- 
tain perturbation formulas, some of which have already 
appeared in the literature by different authors. 


ne 


Purser 


Fig. 1—Cross section of a waveguide with a dielectric slab. 
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Specific Examples 


Let it be required to find the cut-off frequency w, 
of the fundamental mode of a rectangular waveguide of 
width a in which a dielectric slab of thickness 6 is placed 
adjacent to one of the walls as shown in Fig. 1. The slab 
has a permittivity e different from €, the permittivity 
of air. Its permeability, however, is assumed to be wo, 
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the permeability of air. Now, the cut-off frequency of a 
waveguide is the resonant frequency of a two dimen- 
sional cavity defined by the cross section of the guide. 
Let us then use (1) with E =a, sin (7x/a) as a trial field. 
The result is 


¢ 7 : 
We Moeo = | — 
a 
(e are €0) 6 1 i 276 ie 
i — sin : 
€0 a Qa a 


Upon introducing the cut-off wavelength i, related to 
we by w.2uo€o = (27/X-)? and computing a/d. for various 
thicknesses of the dielectric slab and for a dielectric 
constant of 2.45 (polystyrene) Table I is obtained. The 
results are compared to the exact values given in a sepa- 
rate column and read from published curves.’? Note the 
simplicity of (17) and its possible use as a design equa- 
tion. The exact treatment requires the solution of a 
transcendental equation. 


(17) 


DABLE 


Cut-OFF WAVELENGTH OF A RECTANGULAR WAVEGUIDE 
WITH A DIELECTRIC SLAB 


6/a a/X, (exact) a/X- (approximate) 
0.000 0.500 0.500 
0.167 0.485 0.486 
0.286 0.450 0.455 
0.500 0.375 0.383 
0.600 0.355 0.352 
1.00 0.318 0.320 


Next we shall calculate the propagation constant of 
the fundamental mode of a rectangular waveguide with 
a dielectric slab located symmetrically as shown in Fig. 
2. The slab has a dielectric susceptibility x. = (€/eo) —1 
and a magnetic susceptibility zero. Utilizing the varia- 
tional expression (16) and taking E=a, sin (7x/a), asa 
trial field, we obtain 


yu? = 1 «(2 et 44 ~) — 25(a/r-2, (48) 
a 7 a 


where X, A, are the free-space and waveguide wave- 
lengths, respectively. Eq. (18) is plotted for various 
thicknesses of the dielectric slab and the results com- 
pared in Fig. 2 to exact solutions given in reference 5. 
A rectangular waveguide with a ferrite slab off center 
as shown in Fig. 3 is important as a differential phase 
shifter. The differential propagation constant is defined 
as the difference of the propagation constants in the two 
opposite directions, y4-y-. This difference is a conse- 
quence of the tensor susceptibility of the ferrite, which 
for the configuration shown in Fig. 3, is of the form: 


He eK 
Xe OL AO Outs 
KI Oa lee 


7C. G. Montgomery, R. H. Dicke, and E. M. Purcell, “Principles 
of Microwave Circuits,” M.I.T. Rad. Lab. Series, vol. 8, p. 387; 1948. 
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Fig. 2—Variation of propagation constant vs frequency for 
slabs of various widths in a rectangular waveguide. 


An important quantity is the displacement of the slab 
from the wall, d, for which the differential propagation 
constant is a maximum. To determine this quantity, we 
must first determine the propagation constant. Let us 
utilize, for this purpose, (13) with a trial field selected 
as follows: For the E-field we take, 


Re ene 
€ =, (sin —+A sin =), (19) 


a a 


— —— EXACT 
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Fig. 3—Differential propagation constant as a function of the 
location of a ferrite slab in a rectangular waveguide. 


where A is a variational parameter. We determine 3¢ by 
substituting (19) into (11) where for this equation only 
Mis assumed to be a scalar and equal to po and ¥ is as- 
sumed to equal the propagation constant of the funda- 
mental mode with the ferrite slab replaced by a dielec- 
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tric one® having the same permittivity as that of the 
ferrite, but with a permeability equal to that of the free 
space. Substituting E and 3¢ in (13), determining A from 
the condition dy/0A =0, and substituting the value of 
A, thus determined back in (13), y is obtained. Note 
that the immediately preceding procedure involves sev- 
eral simplifying assumptions and it amounts essentially 
to utilizing the variational principle (13) subject to the 
constraint (11). Details have been omitted because they 
are lengthy. The differential propagation constant 
plotted as a function of the distance d and for specific 
values of slab thickness and tensor permeability is 
shown in Fig. 3 and compared with the exact solution 
obtained by Lax, Button, and Roth.® Note the agree- 
ment in the location of the maximum. The ratio of slab 
thickness to the width of the guide is 0.044. For ratios 
appreciably less than 1 per cent, we can simply substi- 
tute for € and 3¢ in (13) the expressions valid for an 
empty waveguide taking into account, however, the 
continuity of the normal component of the magnetic 
flux density at the surface of the slab. The result is 
aw O K 2nd 
Ve yD sin ) 
oe a 


(20) 


where 6 is the thickness of the slab. 
This is the expression Lax, Button, and Roth ob- 
tained from an expansion of a transcendental equation.® 


Fig. 4—Cavity with deformed boundary. 


Perturbation Formulas 


Cavity Surface Perturbation: For a general application 
of (6) consider the following problem: A resonator 
bounded by conducting surfaces S, S; is deformed so 
that surface S; is replaced by S2: as shown in Fig. 4. If 

8 A rectangular waveguide with a dielectric slab can be approxi- 
mately determined by the methods discussed in this paper. _ j 

9 B. Lax, K. J. Button, and L. M. Roth, “Ferrite phase shifters in 


rectangular waveguide,” J. Appl. Phys., vol. 25, pp. 1413-1421; No- 
vember, 1954. 
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the field configuration and the resonant frequency for 
the undeformed cavity are characterized by E, H, w an 
approximate expression for the frequency of the de- 
formed cavity is obtained by substituting E as the trial 
field in (6). Note that £ violates the boundary condition 
at S: but it is an admissible trial field because of the 
presence of the surface integral in (6). The result is 


f [H*-y-H — E*-e-E|dv 


| [ Beka 


where v’ is the volume bounded by S; and Sy. When pu 
and e€ are scalars (21) becomes identical to a formula 
derived by Slater by a different method.” 

Waveguide Wall Perturbation: Let E, 3C, y represent 
the known solution for a certain waveguide. Suppose, 
now, that the cross section is deformed so that s is the 
cross-sectional surface between the original and final 
cross sections and let y,; be the new propagation con- 
stant. Substituting €, 3 as trial vectors in (14) we 


obtain 
f set-us seas — f E+ Eas 
Ss § 


| gct-a. x Eds — [ 6-0, x 5Cds 
sy Ss 


The integrals of the denominator are over the cross- 
sectional surface of the deformed waveguide. 

Cavity Volume Perturbation: Let a resonator be char- 
acterized by w, E, H, when the permittivity and the 
permeability of the medium are e¢, u. Let w; be the reso- 
nant frequency when e, » change to &, fu, the boundary 
of the resonator remaining the same. Substituting /, // 
as trial fields in (9) we get 


[ wewHtde + [ Be Bae 


(21) 


v—-y=o (22) 


(oh == ©! (23) 
[ we u-tao + [ Beta 
or, forming the relative frequency shift, 
Aw 01, — 
Wy Wy 
fn — w)-Hdv + [ BG — ¢): Edv 
= — ~(24) 


[we w-tas + [ Be Bao 


Eq. (24) can also be obtained from a well known for- 
mula.?? 


10 J. C. Slater, “Microwave Electronics,” D. Van Nostrand 
Company, New York, N. Y., p. 81, 1950. 

11H. B. G. Casimir, Philips Res. Rep., vol. 6, pp. 162-182, 1951. 
Essentially the same formula was also derived by Bethe and Schwinger 
in a Rad. Lab. Report of M.I.T., March 4, 1943. 
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Waveguide Volume Perturbation: Let a known solution 
be y, ©, 5 for a certain waveguide when the permit- 
tivity and permeability of the medium are e, uw. Let us 
find y,, the propagation constant when e, wu are changed 
to €, w1, the boundary of the waveguide remaining the 
same. Substituting €, 5 as trial fields in (13) we obtain 


Pea 1—€) eds+ f 5*- (ui — 


fx A ChayOs eis [a @xatds 


- Keds 


# i aes 9 f - (25) 


Eq. (25) can also be derived from a formula first de- 
rived, but not published, by B. Lax of Lincoln Labora- 
toryee Mele Te? 
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Discussion 
COMMENTS BY THE REFEREE 


It should be noted that formulas similar to those given 
by Berk, but applicable to the general lossy case, can be 
derived by using the reaction concept.'* The method 
can be illustrated by considering the problem of calcu- 
lating the resonant frequency (and bandwidth, if the 
loss is not zero) in terms of an assumed distribution of E 
in a closed cavity. To apply the reaction concept we 
have to find some source which would generate the as- 
sumed distribution. Let the source distribution be taken 
as a volume distribution of electric dipoles J plus a 
surface distribution of magnetic dipoles K, where 


VX H= (joe toE+J (26) 
VX E= — jopH (27) 
; el 
. J = — Gwe t+ o) E+ —V X (VX LE) (28) 
Ww 
and 
IKe = 72 d< 12'S) (29) 


where S is the boundary of the cavity and 7 a unit nor- 
mal vector pointing outwards. Inspection of the above 
equations shows that the combination of J and K, gen- 
erates the assumed E£ in the cavity, for Maxwell’s equa- 
tions are satisfied at every interior point and the 
boundary conditions on S are satisfied. 

Using the method and notation of the reaction con- 
cept!’ this leads to the equation 


(a,0) = 0 (30) 


for the resonant frequency (and bandwidth) where “a” 
represents the approximate source consisting of the 


2 Private communication. 

3 V. H. Rumsey, “Reaction concept in electromagnetic theory,” 
Phys. Rev., vol. 94, pp. 1483-1491; June, 1954, and vol. 95, p. 1705; 
September, 1954. 
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combination of J and K,. [It can be verified that the 
value of w obtained from (30) is stationary for varia- 
tions of the assumed E about the correct distribution 
(cf. ref. 1).] Written out in detail (30) gives: 


oa ff [yn ff wonas 

Vv S 

— [[ f #-Gee + 0) Ban 
V 

+ ff LEW x WV x Bao 
Vw 

== ff ex wv x Bnd. 

w® Ss 


If c=0 and E£ is real, (31) reduces to Berk’s (6). If 
o £0, (31) gives a complex value of w. The real part is 
the resonant frequency and the imaginary part is simply 
related to the bandwidth. 

To derive a formula based on an assumed distribution 
of E and H we take the source distribution as a volume 
distribution of electric dipoles J and magnetic dipoles K 
plus the surface distribution of magnetic dipoles A, 
given by (29). In this case J is given by (26) and K is 
given by 


(31) 


Kile jo eC (32) 
Substitution in (30) gives now 
o= iste) Us Ke Wp ff Ke nas 
v s 
= eieh [-E(jwe + E+ EVXH 
-- 
+ HV XE + joHpl \de 
_ ff 2x Horas. (33) 
s 


If £ is pure real, H is pure imaginary and ¢ =0, (33) re- 
duces to Berk’s (10). 

The corresponding formulas for the propagation con- 
stant of a waveguide follow from (30) in like manner, 
except that the formula for the reaction is different. 3:4 
Starting from an assumed distribution of E and H over 
the cross section of the waveguide, (30) gives 


0 = Jf cove + Ks3x)dS +f Ksseal 


= {fi- Es(jwe + cE — sjou 5 


+ EsV xX 3 — 5esV & EldS 


+ Rs(n X E)dl 
G 


ins f { [Es(@ x 30) — ses(@ + €) Jas (34) 


“V. H. Rumsey, “Travelling wave slot antennas,” J. Appl 
Phys., vol. 24, pp. 1358-1365; November, 1953. 
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where 
E = e-iv€(x) ete. 
s represents the matrix 


a0 
Orman 0); 
OneEO 


n = unit outward normal to the waveguide wall, 2 = cross 
section of waveguide, C=cross section of wall, and 
@=unit vector in the z direction. Eq. (34) corresponds 
to Berk’s (14). Equations corresponding to Berk’s (15) 
and (16) can be derived in a similar manner. 


AUTHOR’S REPLY 


I was pleased to see that some of the formulas appear- 
ing in “Variational Principles for Electromagnetic 
Resonators and Waveguides” (VPERW) can be ob- 
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tained from V. H. Rumsey’s “reaction concept.” (See 
Referee’s Comments.) In this connection the following 
remarks are in order. 

Professor Rumsey’s formulas are variational provided 
the dielectric constant and permeability are, at most, 
symmetric tensors. This means that isotropic materials 
and materials with crystalline anisotropy, lossy or non- 
lossy, are allowed while gyrotropic media such as mag- 
netized ferrites and magneto-ionic media are excluded. 

The formulas of (VPERW), on the other hand, are 
restricted to nonlossy substances but they are appli- 
cable to gyrotropic media or media, in general, whose 
dielectric constant and permeability are hermitian 
tensors. Thus each set of formulas has its own sphere of 
applicability. If the complex conjugate sign is dropped 
the formulas of (VPERW) still remain variational pro- 
vided their application is restricted to those cases where 
the expressions derived from the “reaction concept” are 
applicable. 


Nonlinearity of Microwave Ferrite Media* 
N. G. SAKIOTIS}, H. N. CHAITT, anp M. L. KALES} 


Summary—Existing theories of propagation in magnetized ferrite 
media predict propagation constants which are independent of the rf 
field strength only if a number of restrictive conditions are satisfied. 
In general, however, the propagation constants are functions of the 
rf field strength and nonlinear propagation is to be expected. One of 
the conditions for linearity is that the rf magnetic field intensity be 
small compared to the static magnetizing field intensity. This condi- 
tion can be violated when the peak power level of the wave incident 
on the ferrite medium is sufficiently high and indeed, it has been 
observed by the authors that ferrite loaded waveguides can become 
nonlinear at peak power levels as low as 1 kw. 

Results are presented of a study of the behavior of ferrite loaded 
waveguides at 9,400 mc over input peak power levels from 0.1 kw to 
100 kw which includes power levels commonly encountered in radar 
applications. The dependence upon input power level of the ferrite 
losses, phase shift, and rotation of the plane of polarization is de- 
scribed for a number of ferrite materials. The effect upon the degree 
of nonlinearity of the intensity of the static magnetizing field is dis- 
cussed as well as the dependence upon the dimensions of the ferrite. 


OLDER! and others have developed a theory 
Precciin with the propagation of electro-magnetic 

waves in a magnetized ferrite medium. According 
to this theory, the propagation will be linear with re- 
spect to the rf field strength only if a number of restric- 
tive conditions are fulfilled. Under the most general 
conditions, however, nonlinear propagation is to be ex- 
pected. Among the assumptions which lead to a linear 
theory is one which requires that the rf magnetic field 


* Manuscript received by PGAP August 17, 1955. ‘ , 

+ Microwave Antennas and Components Branch, Electronics 
Div., Naval Res. Lab., Washington, D. C. 

1 DPD. Polder, “On the theory of ferromagnetic resonance,” Phil. 
Mag., vol. 40, pp. 99-115; January, 1949 


should be small compared to the static magnetic field. 
Another is that the damping constant is independent of 
power level. The first of these assumptions will clearly 
be violated if the peak power level is sufficiently high. 
The theory in its present form does not indicate how 
large the rf magnetic field must become before an ap- 
preciable departure from linearity can be observed. For 
this reason the authors have conducted tests to deter- 
mine if any nonlinearity is evident at power levels which 
are normally encountered in radar applications. Actu- 
ally it turns out that the propagation becomes non- 
linear at relatively low power levels. 

This situation is very much of practical interest since 
the authors have observed that the nonlinear behavior 
can occur at peak power levels as low as 0.7 kw and the 
absorption loss can become large enough at peak power 
levels less than 100 kw so as to make ferrite devices un- 
usable in many radar applications, even though these 
devices appeared low-loss when tested at klystron power 
levels. 

An example of the increase in absorption can be seen 
in Fig. 1. Here we have a ferrite cylinder, 0.0225 inch 
in diameter, 6 inches long, mounted on the axis of an 
X-band circular waveguide and magnetized by a static 
magnetic field in the direction of the axis. The wave is 
circularly polarized and the loss due to absorption of 
energy by the ferrite is plotted, for each of the two 
senses of polarization, as a function of the peak power 
incident upon the ferrite. 

The positive sense of circular polarization refers to 
that polarization whose configuration rotates, with 
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time, in the direction that would advance a right- 
handed screw in the direction of the applied magnetizing 
field. 

It can be seen from the curve for the wave circularly 
polarized in the negative sense that as the peak power 
level is increased, the absorption loss remains constant 
at 1.2 db, which is the absorption loss observed at klys- 
tron power levels, until the peak input power reaches 
about 1 kw. As the input power is increased beyond this 
level the fraction of the input power absorbed, in db, 
rises, reaching 18 db at the 100 kw level. The wave 
polarized in the positive sense typically exhibits the 
same behavior except that the power level at which the 
absorption begins to increase is higher, approximately 
10 kw, and the rate of increase with input power level is 
smaller. 
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The linear theory for a ferrite medium magnetized 
along the direction of propagation predicts a different 
effective permeability for each of the two senses of circu- 
lar polarization. This difference in permeability gives 
rise to a difference in the rf field strength in the ferrite 
for the two senses of circular polarization. Therefore the 
nonlinear behavior should be different for the two cases. 
It is clear from the curves that the nonlinear effects do 
indeed depend on the sense of circular polarization. 

It is interesting that the Faraday rotation, which is 
due to the unequal phase velocities of the two senses of 
circular polarization, has been observed to be independ- 
ent of the power level up to a peak power level of 100 kw 
for the materials and dimensions we have tested. 

During this experiment the magnetizing field was 
kept constant at 480 oersteds which is the magnitude of 
applied field necessary to produce the maximum amount 
of Faraday rotation for the ferrite samples that will be 
discussed. 

For the data to be discussed, the power source was a 
magnetron, operating at 9,400 mc, modulated by 1 mi- 
crosecond pulses at a repetition rate of 1,000 cycles. Thus 
the average power is one-thousandth of the peak power. 
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Since the propagation characteristics of a ferrite 
medium are also temperature sensitive, care had to be 
taken to insure that the observed variation of attenua- 
tion with power level was not due to a rise in tempera- 
ture of the ferrite, since the ferrite samples had to dissi- 
pate 100 watts at times. This was checked in a number 
of ways. One such check involved a comparison of the 
high power data with temperature data. The attenua- 
tion and rotation were first measured as a function of 
power level. Then with the power level held at a low 
constant value the attenuation and rotation were meas- 
ured as a function of temperature. It was found that 
when the temperature was raised to a level high enough 
to obtain losses comparable to those observed at high 
power levels, the change in rotation was very large com- 
pared to the change in rotation resulting from an in- 
crease in power level. If the heating of the ferrite at high 
power had been the primary cause of the observed 
changes in properties, one would expect the changes in 
rotation to be comparable. For example, it was found 
that for a ferrite rod 0.225 inch in diameter and 2 inches 
long the loss increased from 0.5 db to 3 db, and the ro- 
tation changed at most 5° when the peak power level 
was increased from 0.4 kw to 10 kw. In order to get a 
3 db loss at the low power level the temperature had to 
be raised to more than 550°F. At this temperature the 
change in rotation was 120°. 
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For a given input power level, the magnitude of the 
rf magnetic field intensity in the ferrite will also be a 
function of the ferrite diameter. As the diameter is de- 
creased, the magnitude of the rf field intensity decreases. 
Thus one would expect the power level at which the. 
ferrite becomes nonlinear to be dependent upon the 
ferrite cylinder diameter. Fig. 2 shows the observed de- 
pendence on diameter for both senses of polarization 
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when the ferrite is magnetized by an applied field of 480 
oersteds and also for the unmagnetized case. It is seen 
that as the diameter of the ferrite cylinder decreases, the 
amount of rf power required to bring about the non- 
linear behavior of the absorption loss increases. This is 
especially evident in the case of the negative circular 
polarization. The 0.150 diameter cylinder maintains its 
low-power absorption characteristics up to 60 kw while 
the 0.225 diameter cylinder becomes nonlinear in the 
vicinity of 0.7 kw and at 100 kw the absorption loss has 
increased by 16 db. 

The tendency toward more linear propagation as the 
diameter decreases is also exhibited in the curves for the 
positive sense of polarization and in the unmagnetized 
case. 

Another factor which affects the nonlinear behavior 
of the ferrite is the intensity of the applied static mag- 
netizing field. In fact, if the ferrite diameter is small, it 
can be shown that the rf magnetic intensity in the un- 
magnetized ferrite decreases upon application of a mag- 
netizing field, for the negative sense of polarization, and 
increases, for the positive sense. It might therefore be 
expected that, for rods of sufficiently small diameter, the 
input power level necessary for the incidence of the non- 
linear behavior is higher for waves circularly polarized 
in the negative sense than it is for waves circularly 
polarized in the positive sense. This behavior can be 
seen in the curves for the smallest diameter rod. 
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The effect of the applied static magnetizing field in 
the case of a larger diameter rod is shown in Fig. 3. 
Each of the curves represents the absorption loss as a 
function of the applied magnetizing field strength with 
the incident peak rf power level kept constant. It can 
be seen from these curves that the change in absorption 
loss between the lowest and highest power levels in- 
creases with increasing applied field, reaching a maxi- 
mum at approximately 750 oersteds, then decreases 
until at an applied field of 1,750 oersteds it disappears. 
Above this value of applied field strength there appears 
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to be a small amount of change, the absorption loss de- 
creasing with an increase in the rf power level. Thus the 
effect of nonlinearity on absorption seem to be greatest 
in the region of applied field strength below 1,000 oer- 
steds, which is the region most commonly used in the 
design of microwave ferrite phase shifters and rotators. 

Fig. 4 shows the same type of data taken for the wave 
circularly polarized in the positive sense. Qualitatively 
the behavior is much the same as for waves circularly 
polarized in the negative sense for applied fields below 
2,000 oersteds. 
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Ferromagnetic resonance absorption occurs for this 
sense of circular polarization, and we see the character- 
istic resonance loss peak at the higher power level is 
smaller and narrower. This decrease in absorption at 
the resonance field strength as the rf field intensity is 
increased has also been observed by Bloembergen? and 
his co-workers using cavity techniques. 

An estimate of the dependence of the effective attenu- 
ation constant upon the microwave power level can be 
obtained by studying the absorption loss as a function of 
the length of the sample, maintaining the sample diame- 
ter and the magnetizing field strength constant. 

In Fig. 5 are shown curves of absorption loss vs ferrite 
cylinder length, for the negative sense of circular polar- 
ization, the diameter held fixed at 0.225 inch. For 
each curve the peak power incident on the ferrite cylin- 
der is constant. Looking at the 20 kw curve, the absorp- 
tion in db per unit length is almost constant up toa 
length of about 2 inches. Beyond this length, there is an 
interval in which the absorption loss in db per unit 
length decreases rapidly with length. As the length of 
the sample is further increased, or if the incident power 
level is decreased, the absorption per unit length is 
again constant. Since the power level has decreased by 
about 8.5 db in traversing the first two inches of the 
ferrite, the power level at that point is about 3 kw. 
Therefore for power levels between 20 kw and 3 kw it 


2 N. Bloembergen and S. Wang, “Relaxation effects in para- and 
ferromagnetic resonance,” Phys. Rev., vol. 93, pp. 72-83; January, 
1954. 
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appears that the effective attenuation constant does not 

change very much with rf power level. It can be seen 

from the curves of Fig. 5 that the effective attenuation 

constant is also constant at a small value for power 

levels below 1 kw. Between 1 kw and 3 kw there appears 

large change in the effective attenuation constant. 
Neg. Cir. Pol. Wave 
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Fig. 6 shows a curve of the effective attenuation con- 
stant for the 0.225-inch diameter ferrite rod computed 
from the slopes of the absorption vs length curves. The 
transition region from the low loss to the high loss be- 
havior between 1 kw and 3 kw can be seen quite clearly. 

Thus we can visualize what happened in the rod as 
follows. If the incident power is more than 3 kw peak, 
the attenuation constant of the material is at a high 
value over that length of the material which is necessary 
to drop the power level to 3 kw. This is followed by a 
relatively short length of material in which the attenua- 
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tion constant is decreasing rapidly. After the power 
level had dropped to about 0.7 kw peak, the attenuation 
constant in the remainder of the length is again appre- 
ciably constant at a relatively small value. Therefore as 
the input power level is changed, keeping the ferrite 
dimensions the same, the percentage of the total length 
which has a high attenuation constant changes, result- 
ing in a total absorption which varies with input power 
levels in the manner of the experimental curves. 
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Qualitatively, much the same nonlinear behavior is 
exhibited by ferrites in rectangular waveguide. Fig. 7 
shows the absorption loss of a ferrite slab in X-band 
rectangular waveguide plotted as a function of the ap- 
plied magnetizing field for a number of input peak 
power levels. The range of applied field strength repre- 
sented here is from —1,500 oersteds to +1,500 oersteds. 
The increase in absorption loss with increasing power 
level is quite evident. As in the case of the circular wave- 
guide configuration, it was found that as the slab is 
made thinner, the incidence of the nonlinear behavior 


1956 


occurs at higher power levels. This is shown in Fig. 8. 
Here the applied magnetizing field is held constant at 
750 oersteds and the loss is plotted as a function of the 
input peak power level for both directions of magnetiza- 
tion. It is seen that the curves for the thinner piece be- 
gin to rise at a higher value of incident power. 

During the past few years, considerable effort has 
been made to produce better low-loss ferrite materials 
for microwave applications. Much of the evaluation of 
new materials has been made at low klystron power 
levels. The presence of the nonlinear effects makes it 
necessary to produce materials that are not only low- 
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loss at very small power levels but low-loss at those 
power levels at which they will ultimately be used. 

Similarly, any design of a ferrite component should 
be tested at the power level and with the waveguide 
geometry of the final product since the rf field intensity 
will in general depend upon the waveguide-ferrite geom- 
etry as well as the power level of the incident wave. 

Although ferrites exhibiting the nonlinear behavior 
which has been described would not be suitable for many 
applications, it may well turn out that new applications 
which make use of the nonlinear behavior will be de- 
vised. This possibility is under study. 


Diffraction of Electromagnetic Waves Caused by 
Apertures in Absorbing Plane Screens* 
HANS E. J. NEUGEBAUER{ 


Summary—The perturbation field caused by holes in a plane, 
infinitely thin screen of arbitrary material on which electromagnetic 
radiation is incident, can be split up in two fields, one of which is 
symmetrical in the tangential components of the electric, the other 
in those of the magnetic vector. Both fields satisfy equations within 
the holes which are generalizations of Bethe’s conditions for the 
perfectly conducting screen and which can be used to derive approxi- 
mate solutions of Kirchhoff type. The solution for a screen which 
absorbs practically all of the incident radiation is essentially different 
from Kirchhoff’s result for scalar radiation incident on a so-called 
perfectly absorbing screen. 


SHORT REVIEW OF THE THEORY FOR THE 
PERFECTLY CONDUCTING SCREEN 


HIS PAPER is concerned with relations for the 
"| eesneagente diffraction field caused by holes 

in a screen of any material, which are generaliza- 
tions of Bethe’s! conditions for the perfectly conducting 
screen. 

Bethe’s conditions constitute a statement about the 
perturbation field and, consequently, about the total 
field in the aperture. Being rigorous for the ideal case of 
a perfectly reflecting screen of zero thickness, they can 
be used to derive an integral equation for the diffraction 
field.? Another application is to the problem of calculat- 
ing the perturbation field which can be formulated as a 
mixed boundary value problem.* Rigorous solutions are 


* Manuscript received by the PGAP, August 4, 1955. Revised 
manuscript received January 12, 1956. The basic ideas of this paper 
are contained in this author’s Technical Report No. 25, November 
1952, Eaton Electronics Res. Lab., McGill University, Montreal, 
Se under USAF Contract AF 19(122)-81 through Cambridge 

es, Ctr. 

} Formerly, Adalia Ltd., Montreal; now, Res. Labs., RCA Victor 
Co., Ltd., Montreal. 

1H. A. Bethe, “Theory of diffraction by small holes,” Phys. Rev., 
vol. 66, pp. 163- 182; October, 1944. 

2B. B. Baker and E. T. Copson, “The Mathematical Theory of 
Huygens’ Principle,” Clarendon Press, Oxford, Eng., p. 164; 1950. 

3G. Bekefi, “Diffraction of electromagnetic waves by an ‘aperture 
in an infinite screen,” J. Appl. Phys., vol. 23, p. 1403; December, 
1952, “Diffraction of electromagnetic waves by. an aperture ina large 
screen, ” vol. 24, 1123; September, 1953. 


known for the infinite half-plane* and for a circular 


aperture’ in a perfectly reflecting screen. Yet, the greater 
value of these conditions seems to lie in the fact that 
they can be used as basis for approximate solutions.® 
Such solutions which can also be derived by different 
methods’~® are of a similar nature as Kirchhoff’s solu- 
tions of scalar diffraction problems and can be obtained 
for various shapes of the holes. An extensive experimen- 
tal study under G. A. Woonton has shown that theory 
and experiment are in good agreement providing the 
apertures are larger than a wave length.!° Therefore, 
this approach seems to have definite advantages. 

The generalized Bethe’s conditions to be derived in 
this paper hold rigorously for screens of any material 
providing screen thickness is small compared to free 
space wave length. They can be used, in a very similar 
way to the one mentioned above, for obtaining approxi- 
mate solutions for the perturbation fields caused by 
holes in a screen. 


4 A. Sommerfeld, “Ueber die ausbreitung der wellen in der draht- 
losen telegraphie,” Ann. Physik, vol. 28, p. 665; March, 1909; vol. 81, 
p. 1135; December, 1926. 

5 J. Meixner and W. Andrejewski, “Strenge theorie der beugung 
ebener elektromagnetischer wellen an der vollkommen leitenden 
kreisscheibe und an der kreisfoermigen oeffnung im vollkommen leit- 
enden ebenen schirm,” Ann. Physik. (6), vol. 7, p. 157; April, 1950. 
W. Andrejewski, “Strenge theorie der beugung ebener elektromag- 
netischer wellen an der vollkommen leitenden kreisscheibe und an 
der kreisfoermigen oeffnung im vollkommen leitenden ebenen schirm,” 
Naturwiss., vol. 38, p. 406; September, 1951. 

6S. Silver and M. J. Ehrlich, Antenna Lab. Reports No. 181 
(1951) and No. 185 (1952), Dept. of Engrg., University of California. 

7H. Severin, “Beugung elektromagnetischer zentimeterwellen an 
metallischen blenden,” Z. Naturforsch., vol. 1, p. 487; September, 
1946. 

8H. E. J. Neugebauer, Eaton Electronics Res. Lab. Report No. 
b5 on contract to the USAF, Cambridge Res. Ctr. (1952); “A new 
method of solving diffraction problems,” J. Appl. Phys., vol. 23, p. 
1406; December, 1952. 

9S. J. Buchsbaum, Eaton Electronics Res. Lab. Report No. b9 
on contract to the USAF, Cambridge Res. Ctr., 1953. 

10S. J. Buchsbaum, iene Milne, D.C. Hogg, G. Bekefi, and G. A. 
Woonton, “Microwave diffraction by apertures of various shapes,” 
J. Appl. Phys., vol. 26, p. 706; June, 1955. 
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Some introductory remarks on the theory for the per- 
fectly reflecting screen may be given before the general 
case will be discussed in the succeeding sections. 

If electromagnetic radiation of harmonic time de- 
pendence, exp (jwt), is incident from the left side, z <0, 
on a plane and perfectly conducting, infinitely thin 
screen with holes of any shape and size, the total field 
before and behind the screen can be described as the 
sum of the incident, reflected and perturbation field. 


=> = => 
Since incident and reflected field, E‘, H* and E’, H’, are 
those which would exist if the holes were closed, they 
are zero behind the screen, 2>0. If the holes are open 


the perturbation field, E>, HP, is added to the incident 
and reflected field, and the tangential components E,? 
Ey,” of the electric vector as well as the normal compo- 
nent H,” of the magnetic vector are shown to be sym- 
metrical with regard to the plane of the screen, z=0. 
Other components E,?, H,?, H,? are antisymmetrical. 

These symmetry relations can be derived from Strat- 
ton and Chu’s vector Green’s theorem.!! The simplest 
way, however, to establish them is to add fictitious 
sources lying symmetrical to the real sources with re- 
spect to the screen and radiating with opposite phase.” 
The incident and reflected fields of the fictitious sources 
are the continuations, through the plane z=0, of the 
reflected and incident fields respectively of the real 
sources so that both types of sources together do not 
produce any perturbation field. 

Therefore, the perturbation field of the fictitious 
sources at any point is opposite that of the real sources 
at the same point. In addition, due to the phase relation 
mentioned above, the perturbation field of the fictitious 
sources is antisymmetrical to that of the real sources in 
the E,, E,, and H, components, but symmetrical in the 
E,, H,, and H, components. Hence follow the symmetry 
properties of the perturbation field. 


> o> lO 
The further condition that the total field, ‘+ 47+ E?, 
—-> = 


ois 
H'+H'’+H?, must be continuous in the holes, together 
with its first derivatives, leads to simple equations relat- 
ing the perturbation field to the incident and reflected 
field. They are known as Bethe’s! conditions and are a 
special case of (5). 


SYMMETRICAL AND ANTISYMMETRICAL 
PERTURBATION FIELD 


Let it now be assumed that the screen consists of any 
material and that its thickness is negligible compared to 
free space wave length. The further assumption that the 
screen may be opaque is made only to simplify the for- 
mulas since this theory can be carried through in essen- 
tially the same way for a transparent screen. The required 
additions for a transparent screen will be mentioned. 

The total field is the sum of the incident, reflected and 
perturbation field. 

1 J. A. Stratton and L. J. Chu, “Diffraction theory of electromag- 


netic waves,” Phys. Rev., vol. 56, pp. 99-107; July, 1939. 
12 B. B. Baker and E. T. Copson, loc. cit., pp. 159-160, 163-164. 
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> lS 


\_S Es + Ee + Bp (1) 
H = Hi+ H+ Ae. 


Incident and reflected fields may be known. 

Only the coordinate z of the field point will be written 
in the following equations while « and y will be sup- 
pressed as they are always the same for all field func- 
tions occurring in the same equation. 

The perturbation field can be split up into the sum of 


two fields, Es ie and Es, Hr, by means of the following: 


# Ez y'(@) = 4 [Beve) eee 

E,*(2) = $1 £,@) = E29) 

isha tay opel) 912i) 
Le (2)e= 2 LG) ee ao 2) 
Exy°(2)-= 4|E2a(2)— Ex,s(—2) | 
‘| £7@)— 2 |G ea 
« Hz) = 4|H. 2) + Heal 2 

Hyt@) =| 18) oi 


It can easily be verified that 
> —> > > => => 


a) #°=E*:+E* and He = +E; 

b) both partial fields, Es i as well as E*, Ht, satisfy 
Maxwell’s equations in free space; 

c) the field components marked with an asterisk in 
(2) are symmetrical, the others are antisymmetri- 
cal with respect to the plane of the screen, z=0. 


As regards the nomenclature, the tangential field 
components of the electric vector are considered repre- 
=> 
sentative. Hence, the partial field H*, H* will be called 
=> => 


the symmetrical, £*, H* the antisymmetrical perturba- 
tion field. 

For future applications it may be noted that the de- 
rivative, with respect to z, of a symmetrical field com- 
ponent is antisymmetrical and vice versa. 

It may also be noted that Kirchhoff’s expression for 
the perturbation field in the case .of scalar waves is 
known to be the sum of a symmetrical and an antisym- 
metrical function.¥ 


GENERALIZED BETHE’S CONDITIONS 


Applying the continuity condition of the field within 
an aperture to the x-component of the electric vector, 


E.(—0) + E,"(—0) + E.*(—0) + E,*(—0) 
= E,\(+0) + E.(+0) + £.(+0) + E*(+0). (3) 
Because of 
E,*(+0) = £,"(+0) = 0, £,*(—0) = £,'(+0) 
and 
—E(—0) = E,*(+0), 


13 C. J. Bouwkamp, ee ia theory,” Phys. Soc., Prog. Rep. 
in Phys. X VII, p. 52; 1954 
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(3) reduces to 
E.%(+0) = $[E.(—0) + E,"(—0)]. (4) 
Correspondingly, 
Bee 0) = H(—0) -- H,°(—0) + H.*(—0) 
= H,(+0) + H,"(+0) + H.°(+0) + H,%(+0) 
which reduces, because of H,'(+0)=H,’(+0)=0, 
H,§(—0) = —H,'(+0) and H,*(—0) = +H,*(+0), to 
H.*(+0) = 3[H.(—0) + H.*(—0)]. (5) 
All the symmetrical and antisymmetrical field com- 
ponents satisfy similar equations to (4) and (5) respec- 
tively. This holds true also for the derivatives with re- 


spect to 2. 
In the case of a perfectly conducting screen 


eee ges 
E* = He =0 
a) Fay —0), £(—0) = £4 ~0) 
0) = Hog (0), > H"(—0). = — 1, —0) 
and (5) simplifies to one of Bethe’s conditions, 
H,*(+0) = H,*(—0). 


For a transparent screen the unperturbed field con- 
consists of the incident and reflected field, 


=> lS > => 
Be eke ahve be Tah 


in front of the screen, z>0, and of the transmitted field, 


> 
t t 
es, 


behind the screen, z>0. Hence, the following changes 
must be made if the screen is transparent: 


— => 

Et and H'* respectively are added to the right sides 
of (1). 

E,'(+0) is added to the right side of (3) and is sub- 
tracted in the square bracket of (4). Similarly, 

H,,‘(+0) is subtracted in the square bracket of (5). 


It may be noted that the derivations of this section 
are rigorous providing screen thickness and multiple 
scattering between sources and holes can be neglected.** 


APPROXIMATE SOLUTIONS 


Generalized Bethe’s conditions (4) and (5) can be 
used to calculate approximate solutions of Kirchhoff 
type for the diffraction field in cases where the reflected 
field can be determined. The general idea is to assume 
the field in the shadow region near the screen as negligi- 
bly small so that the tangential magnetic components 
of the symmetrical field and the tangential electrical 
components of the antisymmetrical field can be deter- 
mined by means of Green’s theorem, 


1 (6) —jkrop 
eo if Haa'(Q) < [ee wus "ase (6) 


Tv TQP 
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1 ) —4k 
BoP) =— [ E,,*(Q) = J ease (7) 
aperture 


Tv 0ZQ YQP 


where P is a field point behind the screen, gp >0, Q isa 
point within the aperture, rgp is the distance between 
Q and P. Reference is made to Fig. 1. Field components 
Fz,y7(Q) and Hz,y?(Q) are found from (4) and (5). 


Fig. 1 


Although more general cases can be treated, the equa- 
tions for a plane wave normally incident on a circular 
hole may be given here as an example. Let the incident 
field, without the time dependent factor, be 


Hz = Ee exp (—7ke), = 0 
2 (8) 
HH, = — exp (—jkz), H,' = HH; = 0. 
wu 
Then, the reflected field is 
EE,” = —ZEq exp (+jkz), E,” = £,,=0 
ZREo (9) 
Hi = exp (+7kz), HH,’ = dH," = 0, 
we 
where Z is the reflection coefficient. 
From (4) to (9) it follows 
RE (A +Z 0 Texp (—7kr 
ie o( ) [ee as (10) 
A4rwp aperture Oz if 
E,\(1 —Z 0 [exp (—7kr 
an f Pos a (11) 
Aa aperture Oz r 


Apart from the factor $(1+2Z) (10) is the expression 
yielding the total perturbation field for the perfectly re- 
flecting screen. In addition, both factors $(1+Z) and 
4(1—Z) being disregarded, the antisymmetrical field, 
(11), can be derived from the symmetrical, (10), by 
proper interchanging of E and H. Hence, the total per- 
turbation field is a linear superposition of the perturba- 
tion fields caused by a perfect electrical and a perfect 
magnetic conductor. 

One limiting case of special interest is constituted by 
the perfectly reflecting screen, Z =1. The other limiting 
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case is constituted by a screen that absorbs so much of 
the incident energy that the reflected and transmitted 
fields can be neglected, Z=0. Although it is difficult to 
make such a screen whose thickness is negligible com- 
pared to free space wave length, it is of interest to in- 
vestigate the nature of the field in this limiting case. The 
screen would be perfectly absorbing in the sense of 
Kirchhoff’s theory as practically no energy is reflected 
from, nor transmitted through, the closed screen, al- 
though it may not be perfectly absorbing for grazing 
incidence. Since, however, we confine these derivations 
to the case of normal incidence and since we consider the 
perturbation field near the screen as negligibly small, 
this shortcoming of realizable screens should not weigh 
too heavily. 

A screen that transmits or reflects only negligible 
amounts of normally incident radiation may be called 
“zero reflection screen.” It is characterized by Z =0. The 
perturbation field is the sum of the two fields given by 


Eok ) [= ee] 
je diesen 
Aro 02 r 


E oe — 7k 
eas f | 520) (a) 2 as. 
Ao Oz Yr 


This result is essentially different from Kirchhoff’s solu- 
tion for the perfectly absorbing screen in the scalar case 
where also two integrals must be calculated, one of 


which is 
0 [= se) 
poner, 
Oz if 
as above, but the other one is 


ii exp (—jkr) a s 
Yr Oz 


Eqs. (12) and (13) are in agreement with an earlier solu- 
tion derived by this author with the help of the Hertz 
vector II for the symmetrical, and the II-vector for the 
antisymmetrical field. However, they are at variance 
with an earlier solution® based on II alone and which does 
not satisfy the generalized Bethe’s conditions unknown 
at that time. 

Since the total perturbation field is a superposition of 
two fields, one caused by a perfect electrical, the other 
one by a perfect magnetic conductor, and since the latter 
follows from the former by interchanging of E and H, it 
is sufficient to evaluate only (10) for the case of a per- 
fect conductor, Z=1. 

From (10), because of 


H,’ (12) 


and 


(13) 


0 te) 
Ee eo Pe 
1 ) exp (—7kr) 
HP esa ay 
2m OZp aperture r 


“4H. E. J. Neugebauer, Symposium on Microwave Optics, Mon- 
treal, Canada, Paper No. 51, p. 4; June, 1953. 
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By means of Green’s theorem the integral over the 
aperture can be transformed in an expression containing 
an integral over the cylindrical surface whose genera- 
tices are the rays, in the geometrical optics sense, through 
the rim of the hole. The integration along a generatrix 
can be performed leading to a contour integral of the 


form 
i a? + p* — 2ap cos? 

-exp [—jk(a? + p? — 2ap cos 6 + 2%)!/?]d0 
where p= + 4/x?+ 7? and tan 0=y/x. Up to this point 
the calculations are almost the same as those published 
by Rubinowicz. For further evaluation the field points 


are confined to the space near the perpendicular axis 
through the center of the hole, 


(p/a)? <1. 


With a?+2?= A? and terms of higher order in p/a being 
neglected, integral (14) becomes 


2n @ — kA a 
i DGD. d (1 + Uacns ) exp (2 in cos ”) d0 
0 a 


a 


a — pcos 8 


(14) 


which can be solved by means of the Bessel functions 
Jo and J; of first kind. 

Finally the electric field is derived from the tangential 
magnetic components by means of 


_ 
ae 


WE 


as as — 
V‘-H=0 and £.= — << He 

The result in a form suitable for experimental tests, 
shall be quoted for a screen of any material with a large 
circular aperture, ka>>1, and for field points behind the 
screen, z>0, which lie either in the H plane, x =0, or in 
the £ plane, y=0, near the axis through the center of 
the hole, p?<a?. 

a) H-plane (x=0): 


1+Z 1-Z k 
E,? = rh éz + —— — hk, 
2 we 
Leis : 
H-,? = 0 
1+Z 1—Z k 
peel? We ern sce ats 
2 wp. 
b) E-plane (y=0): 
14+Z 1—Z Rk 
sip Boh el 
y, 2 we 
E,? = 0 
H,? =0 
17 1—Z 0k 
Hy? = >= hy tem anita 
2 WU 


* A. Rubinowicz, “Die Beugungswelle in der kirchhoffschen 
theorie der beugungserscheinungen,” Annalen der Physik (4), vol. 53, 
pp. 257-278; August, 1917. 
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with 


€z = Eo {exp (—jkz) — exp (—7kA) 


| sotee) -< a“ } 


és = Ey fexp (—jkz) — exp (—7kA) 


a? a® J,(kB) 

45 (ae Sk a Ny a 

( ) JobB) ie kB } 
RE 


ps ho = —*Sexp ape serp (ik) 
ay) A 


| soca) + ij (1 2s “) a (a) |b 


xa 
xi for the H plane 


a 
B —— Lied = 
A xa 
— for the £ plane, 
A 
k? = wen. 


FURTHER APPLICATIONS 


The theory briefly described in the preceding section 
has been shown to be in good agreement with experi- 
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ments for different sizes and shapes of holes in a per- 
fectly conducting screen.!° Experiments on screens of 
zero reflection have been started by several workers. No 
results have been published yet, probably due to the 
difficulties of obtaining a sufficiently thin screen of high 
absorbing power. 

The present theory can be applied to screens of any 
reflection coefficients. Hence, it could be tested without 
making a perfectly absorbing screen. The theory can also 
be applied to metal screens whose absorption cannot be 
neglected because of poor conductivity. 

In addition it seems possible that the scattering by 
absorbing solid bodies whose dimensions are large com- 
pared to free space wave length can be treated in a 
similar way by first solving the problem for the perfect 
electrical and perfect magnetic conductor and finally 
adding both solutions together. One simple case where 
all the required calculations for the perfect conductor 
have been done is the scattering of a plane wave by a 
circular cylinder. 
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Side-Lobe Suppression by Pattern Multiplication* 
RAYMOND JUSTICE} 


Summary—It is shown that the properties of the radiation pat- 
terns of broadside and end-fire line radiators can be used to mini- 
mize the side-lobe level of the radiation patterns of uniformly excited 
rectangular arrays. 


I. INTRODUCTION 
RSs CERTAIN conditions the radiation pat- 


tern of a rectangular end-fire array can be repre- 
sented as the product of the pattern of an end-fire 
line source and the pattern of a broadside line source. 
This well-known fact has some interesting consequences 
which may not be generally recognized.! If an end-fire 
line source with phase velocity equal to the velocity of 
light is uniformly energized, the highest side lobe in the 
resulting pattern is about 20 per cent of the maximum. 
Similarly, a uniformly energized broadside line source 
* Manuscript received by the PGAP, June 26, 1955. _ 
+ Dept. fo Electrical Engineering, Ohio State University, Colum- 
bus, Ohio. 
1 Final Engrg. Rep. 400-11, December 15, 1951, Antenna Lab., 
Ohio State Univ. Res. Found.; prepared under Contract AF 33(038)- 


9236, Air Res. and Dev. Command, Wright Air Dev. Center, Wright- 
Patterson AF Base, Ohio. 


has a radiation pattern with the same side-lobe level. It 
follows that if the dimensions of a rectangular array are 
chosen so that these two patterns have the same beam 
width between first nulls, the product pattern will have 
side lobes that are not greater than 4 per cent of the 
maximum. The purpose of this paper is to illustrate the 
relationships between the length and width of such a 
rectangular array in order to best utilize the side-lobe 
suppressing characteristics of this multiplication process. 


II. THE END-FIRE AND BROADSIDE PATTERNS 


If the rectangle of the dimensions and orientation 
shown in Fig. 1 is covered by a continuous distribution 
of isotropic point sources that are uniformly energized 
in amplitude, in phase in the Y-direction and with a 
phase velocity in the Z-direction not greater that c, the 
speed of light, then the normalized radiation pattern in 
the Y-Z plane is given by? 


2 Cf., for example, J. D. Kraus, “Antennas,” McGraw-Hill Book 
Co., New York, N. Y., 1950. 
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| £0) | sin X,; sin =| XG | (1) 
* Xi Xo sin Xo 6—0 
where 
sin X1 


is the pattern of a broadside line source of length W’ 
1 


sin Xo , . , } 
= is the pattern of an end-fire line source of length] L 
2 
NG Salhi’ Sin G 
X_ = tL(c/v — cos 8) 
W 
W=— 
r 
iG 
L=— 
r 


\ = free space wavelength 


Cc velocity of light in vacuum 


v phase velocity in Z-direction 


Fig. 1—Dimensions and orientation of array. 


Curves of the individual patterns are given in Fig. 2 
and Fig. 3 for a pair of values of W and L. Fig. 4 com- 
pares the end-fire patterns obtained from a given length 
for two values of c/v. Certain observations about these 
types of pattern can be made immediately on the basis 
of {these examples. As mentioned in the introduction, 
the side lobes in the broadside pattern are at most 20 per 
cent of the maximum. For the end-fire pattern the same 
result holds if c/y=1, but the side lobes increase rapidly 
in relative magnitude as the phase velocity is decreased. 
From the expression for the end-fire pattern and the fact 
that the function sin X/X has its largest secondary 
maxima at X =37/2, it follows that the side-lobe level 
of the end-fire pattern for values of c/v2=1 is given by 


Side-lobe level 


sin Xo ») 
= | | - | 100 percent ofmain beam (2) 
X2 lpaomLsin X2_lo—0 
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Fig. 3—End-fire patterns L=9, L=19, c/v=1. 
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Fig. 5—Side-lobe level of end-fire pattern as a function of ¢/v. 


Examples of this side-lobe level as a tunction of c/v for 
fixed L are given in Fig. 5. The rapid increase in side- 
lobe level for increasing c/v will ultimately limit the ad- 
vantages of using such values of c/v for beam sharpening 
purposes (see Section IV). 
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The first zero of sin X/X occurs at X =7; hence the 
broadside and end-fire patterns have their first zeros at 
fiw=sin (1/W) and @:,=cos (c/v—1/L), respec- 
tively. This latter result illustrates the beam sharpening 
effect of increasing ¢c/v, and it might be anticipated that 
a judicious use of this effect, while recalling the results 
of the preceding paragraph, could result in a narrower 
main beam for a given side-lobe level than could be 
achieved with c/v=1. Curves of 6:7 as a function of W 
and @;, as a function of Z and c/v are given in Fig. 6. 
Observe that for L in excess of about 10, 6, does not 
decrease rapidly with increasing L. For later use a curve 
of 6m=cos [1—3/2L] is included. 6m is the location of 
the maximum of the first side lobe of the end-fire pattern 
for c/v=1. 


III. THE CasE For c/v=1 


The discussion in this section will be restricted to 
consideration of the patterns in the plane of the array, 
the Y-Z plane, when the phase velocity in the Z-direc- 
tion is equal to c and the only parameters are the array 
dimensions, L and W. For this case both the end-fire and 
broadside patterns have side-lobe levels approximately 
equal to —13.2 db (20 per cent of main beam maxi- 
mum). The condition of coincidence of first nulls, 
6i:w =6i1, assures a pattern of the array with a side-lobe 
level of at most 4 per cent of the main beam. One might 
expect that proper choice of dimensions would give an 
even lower side-lobe level. That is, if the length of the 
array is fixed, the value of W for which minimum side 
lobes are achieved would be related to the fixed lengths. 
Analytical attempts to determine this relation have not 
been successful. However, by calculating large numbers 
of patterns and drawing graphs of side-lobe level as a 
function of one dimension for various values of the 
other dimension, it was possible to obtain sufficient data 
to estimate this relation, at least over the range of 
values for which patterns were calculated. 

The accuracy of the result is limited by several fac- 
tors. Since it was necessary to calculate large numbers 
of patterns it was not deemed practical to plot each 
pattern. The number chosen for the side-lobe level cor- 
responds to the highest calculated value following the 
first minimum. The patterns were calculated at two- 
degree intervals of 6 and the side-lobe levels determined 
from these calculations may contain considerable devia- 
tions from the true values. Furthermore, for a given 
value of one dimension, only a few values, up to seven, 
of the other dimensions were chosen for purposes of 
calculating the array pattern. As can be seen from the 
examples given in Fig. 7, it was frequently quite diff- 
cult with the available data to estimate with accuracy 
the relation between array dimensions for a minimum 
side-lobe level. However, considering all these things, 
the best estimate of the relation between LZ and W for 
minimum side lobes appears to be 


(3) 
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Fig. 7—Side-lobe levels as a function of W. 


for the case c/v=1. This relation corresponds to the 
first null of the broadside pattern occurring at the same 
angle at which the end-fire pattern has its first (and 
highest) secondary maximum. In general it may be 
qualitatively stated that for this case (c/v=1) the side- Fig. 8—Effect of W on rate of change of side-lobe level with L. 
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lig. 9—End-fire and broadside patterns. 


lobe level does not vary rapidly with L in the neighbor- 
hood of the value which satisfies (3). The variation with 
L is slower for the larger values of W as exemplified by 
comparing the curves of Fig. 8 for W=2 and W=4. 
This is to be expected since the value of Z which satis- 
fied (3) for the larger values of W is in the region where 
the curve of 6, as a function of length is varying quite 
slowly. An example of the results obtainable by choosing 
Land W to satisfy (3) is given in Figs. 9 and 10 where 
the two factors of the array pattern and the array pat- 
tern itself are given. 


IV. THE CASE FoR c/v21 


In this case the array dimensions are fixed and the 
variable parameter is the phase velocity v. If this ve- 
locity is decreased to a value less than c, the end-fire 
pattern has a narrower main beam but higher side lobes. 
The problem is to determine when the disadvantages of 
the latter effect essentially outweigh the former. 

The analytical difficulties encountered in the case for 
c/v=1 are compounded by the introduction of phase 
velocity as an additional variable. In this case it was 
necessary to compute even larger numbers of patterns 
than previously in order to obtain sufficient data to 
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Fig, 10—Pattern of an array with dimensions satisfying (3). 


indicate what relations may exist between the three 
parameters for a minimum side-lobe level. Since no 
definitive result was obtained, a qualitative discussion 
must suffice. 

For the narrower widths, say W<2, it generally 
turned out that the optimum phase velocity for a given 
length Z in the sense of minimum side lobes-was about 
the same as that required by the Hansen-Woodyard 
condition for maximum directivity of an end-fire line 
source of length L.* In the present notation this phase 
velocity satisfies 


1 
one atecra (4) 


For wider arrays the broadside pattern may have suffi- 
cient influence on the end-fire pattern to cause the 
optimum c/v to be greater than that given by (4). 

The limitations on the accuracy of the determination 
of side-lobe level as a function of one of the three pa- 
rameters are the same here as in the case for c/v=1. 
Samples of side-lobe level as a function of c/v for a few 
dimensions are given in Figs. 11 and 12. 


3 W. W. Hansen and J. R. Woodward, “A new principle in direc- 
tional antenna design,” Proc. IRE, vol. 26, pp. 333-345; March, 
1938. 
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Fig. 11—Side-lobe levels as functions of ¢/v. 


V. CONCLUSIONS 


It is not intended to imply that the optimum combi- 
nations of L, W and c/v will give a particularly low side- 
lobe level in all cases. The conclusion to be drawn here 
is that for c/y=1 and one dimension fixed, there is an 
optimum value of the other dimension. For this case the 
side lobes under optimum conditions are not more than 
4 per cent of the main beam maximum and are usually 
appreciably smaller than this. If the array dimensions 
are fixed, then there is an optimum value of c/v which 
will give minimum side lobes for an array of these di- 
mensions. In this case it is not possible to give an exact 
upper bound to the side-lobe levels, but on the basis of 
calculations made to date it may be asserted that opti- 
mum conditions will assure side lobes less than 10 per 
cent of the main beam maximum. 
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Fig. 12—Side-lobe levels as functions of ¢/v. 
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On the Conductance of Slots* 
JAMES R. WAIT} 


Summary—The external radiation conductance G of narrow half- 
wave axial and transverse slots on a circular cylinder, a ribbon, a 
knife edge, and a right-angle wedge are computed. The method is 
based on finding the total power radiated assuming a sinusoidal 
variation of voltage along the slot. The theoretical values of G for 
the ribbon of width 2d are applied to the case of a half-wave slot cut 
in the narrow face of a rectangular waveguide which has a flush 
mounted plate or flange of width 2d. 


* Manuscript received by the PGAP, August 25, 1955. 
} National Bureau of Standards, Boulder, Colo. 


N THE DESIGN of slot antennas it is usually as- 
if sumed for purposes of calculating the conductance 

that the surface on which the slot is cut is flat and 
infinite in extent. In practical radiating systems the sur- 
face is often curved to an appreciable extent such as 
when the slot is located in the leading edge of an aircraft 
wing. In other situations, the slot may be cut in the 
broad or narrow faces of a rectangular waveguide which 
have dimensions comparable with the wavelength. 
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It is the purpose of this paper to show how the con- 
ductance is affected by the curvature and the finite ex- 
tent of the surface on which it is cut. To facilitate the 
solution of the problem it is necessary to idealize the 
shape of the surface which is assumed to be perfectly 
conducting. The simplest case is probably the circular 
cylinder of infinite length of radius a. 

The conductance G is now obtained by calculating the 
total power P radiated from the slot and then dividing 
by V?/2 where V is the voltage at the center of the slot. 
Furthermore, it is assumed that the voltage distribution 
along the slot is sinusoidal. Choosing a spherical coordi- 
nate system (R, 0, ¢) such that 6=0, a corresponds to 
fe ee of the cylinder, the far zone fields are given by 


e@tkR 


Ee VP.(8, >) (1) 


and 


etkR 


Es 


I 


VP5(, $) (2) 


where k = 27/wavelength and Ps, and P, are pattern fac- 
tors expressible in terms of a Fourier Series as follows: 


Po(6, 6) = >) Am cos mb + am sin md (3) 


m=0 


and 


P5(9, 6) = >) Bn sin md + bn cos md (4) 
m=0 
where Am, Gm, Bm and bm are coefficients which are a 
function of @. For a narrow half-wave axial slot centered 
at ¢=0, @=7/2 and located on the surface of the 
cylinder [2], 


Am = ha Jess = 0 (5) 
and 
wv 
cos (F cos 0) 
2 erie 
bm = : (6) 
wkasin?@ 4H,,°’(kasin @) 
where 
con Cn = 2m = 0) 


and where H,, is the Hankel function of the second kind 
and the prime indicates a derivative with respect to the 
argument. In the case of the narrow half-wave trans- 
verse slot centered at the same point [3], 

ka Emt™— cos (mr/2ka) 


Am = ; (7) 
rsin§@ H,,2)(ka sin 8) [(ka)? — m?] 


kacot@ - 2mi™ cos (mr/2ka) 


rsin@ Hy» )’(kasin 6) [(ka)? — m?| 


) (8) 


and d@m=bm=0. The conductance is now obtained from 
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Dieses 2?" EE 
Nov o 0 


VY? 
R-o (9) 
where y9=1207 and the asterisk indicates a complex 


conjugate. The integration with respect to ¢@ can be car- 
ried out to yield 


Ga De Ar (22 /enyr | am 2 


no 0 m=0 


>+ | Bm 


+ | bn |?(2/em)] sin 640. (10) 


In the case of the axial half-wave slot 


2 ™ cos? |(1/2) cos A} 2 Ga 
eee | L(r/2) DS Path 
not k?a” 0 sin? 0 m=0 | H,,2”" [2 
and in the case of the transverse half-wave slot 
y) Tee COSHH IG a 
ee peer 
nows Jo sin mao [(ka)? — m?]? 
(ka)’ m? cot@ 7 ; 
| TEoP i: ance (12) 


where the argument of the Hankel functions is under- 
stood to be ka sin 6. Because of the presence of the square 
of the moduli of the Hankel functions in the denomina- 
tors, the infinite series converge quite rapidly. The 
terms in the summation with order m greater than about 
ka are negligible. The integration with respect to @ can- 
not apparently be carried out in closed form. The inte- 
grand, however, turns out to be slowly varying function 
of 6 and the conductance G can therefore be obtained 
readily by a graphical integration. A slight difficulty 
occurs in (12) for the case of the transverse slot since the 
integrand has a logarithmic singularity at 6=0 and 7. 
The contribution to G in this case in the regions 6=0 to 
6 and 6=z—6 to zm is denoted by G; and for 6 small is 
given by 


1 dé 
a~— | - . (13) 
3004 Jo sin A(ka)? | H, (ka sin 6) |? 


where only the term for m=O is significant if 
ka sin 6kab<1. Using the small argument approxima- 
tion for Hankel functions, it follows that 


1 fe do 
~ 30r'(ka)? Jo {1 + (4/n?) [log (Tkad/2) |? 0 


(14) 


Gs 


with log I'=0.5772. The latter integral can be evaluated 
in closed form to yield 


| ne Tkaé ] 
: tan“! Fee = 1/2 | 


where the inverse tangent is a negative angle less than 
7/2. This contribution is then added to the portion of 
the integral obtained by the graphical procedure to ob- 
tain the conductance of the transverse slot. 


(15) 


G; ~ ————— 
6073( ka)? 
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Fig. 1—The computed external radiation conductance for axial and 


transverse half-wave slots on a circular cylinder as a function of 
the circumference of the cylinder in wavelengths. 


The conductance for both axial and transverse half- 
wave slots on a circular cylinder is shown plotted in 
Fig. 1 with an exponential scale for the abscissa. It can 
be seen that the effect of curvature is very pronounced. 
In the case of the axial slot the conductance varies be- 
tween the limits 0.515 and 1.03 millimhos from small to 
large cylinders respectively. It can be readily deduced 
from Babinet’s principle that the 1.03 corresponds to 
the conductance of a thin half-wave slot cut in a flat 
surface of infinite extent and perfect conductivity. In 
the case of the transverse slot the conductance increases 
as ka becomes smaller. 

The above method can also be employed to calculate 
the conductance of slots on elliptic cylinders. Essen- 
tially, the technique is to integrate the total power radi- 
ated from a prescribed slot on a perfectly conducting 
cylinder of elliptic cross-section. The functions involved 
in this case are angular and radial Mathieu functions [4 | 
which play the role of the exponential and Bessel func- 
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tions respectively in the circular cylinder case. Some 
calculations for short slots on elliptic cylinders for spe- 
cial values of the ellipticity have been published by 
Wong [5]. Further results for axial half-wave slots were 
given recently by Wait and Walpole [6] for the case 
when the slot was located at the center of one side of 
ribbon or elliptic cylinder of vanishing minor axis. In 
this latter paper results were also given for a half-wave 
slot located on one side of a half-plane and parallel to 
the edge. The result for the half-plane or knife edge can 
be easily extended to the case when the slot is on the 
surface of a wedge [7]. For example, if the half-wave 
slot is parallel to and distance p from the apex of a right- 
angled wedge, the conductance G is given by 


oe =| ™ cos? [(a/2) cos | 


No” sin 0 


ive} 


> Em [Jom3( kp sin 0) [2de 


m=0,1,2+++ 


(16) 


where Jom3 is the Bessel function of the first type of 
order 2m/3. 

Using the above-mentioned formulas, the conduct- 
ance G for narrow half-wave axial slots is shown plotted 
in Fig. 2 for the cases where the slot is on a knife edge, 
right angled wedge, thin ribbon and circular cylinder. 
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Fig. 2—The computed external radiation conductance of a _half- 
wave axial slot on a knife edge, right angled wedge, thin ribbon 
and circular cylinder. 


In the case of the ribbon the distance p corresponds to 
the half-width of the ribbon. The results for the axial 
slot on the circular cylinder has been replotted in Fig. 2 
for purposes of comparison. 

In Fig. 3 the conductance G is plotted for the case 
where the half-wave slot is moved from one side of the 
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ribbon to the other for a fixed width (kd =8). The corre- 
sponding result for the slot on the knife edge is also 
shown for comparison. The close agreement between the 
results for the knife edge and the ribbon when the slot 
is near the edge is reassuring. Of course as the slot ap- 
proaches the middle of the ribbon, the other edge influ- 
ences the result. : 

To check the results for the ribbon experimentally, 
the slot was cut in the narrow face of a rectangular 
X-band waveguide. Then thin rectangular metal plates 
of various widths 2d were mounted centrally and flush 
with the narrow face of the guide. The equivalent shunt 
conductance g of the slot normalized to the character- 
istic admittance of the guide can then be written in the 
following form 


480 c¢ Ng mh] 1.03 


(kd) = — — “0s"| 
: OXaG 


ily 
737 b X a, 


where 0 and ¢ are the inner narrow and broad dimensions 
of the guide respectively, \, and \ are the guide and free 
space wavelengths, respectively, and G is the external 
or aperture conductance in millimhos of the thin half- 
wave slot. This formula is obtained by a straightfor- 
ward extension of Stevenson’s analysis [7] who assumed 
the external surface of the waveguide to be infinite in 
extent. In other words, when kd approaches infinity, G 
approaches 1.03 millimhos, and the formula for g is 
identical to Stevenson’s result. 
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Fig. 3—The computed external radiation conductance of a half-wave 
axial slot as a function of distance to the edge of a ribbon of fixed 
width with the result for the knife edge shown for comparison. 


The experimental values of g(kd) are shown plotted 
in Fig. 4 expressed as an increment to the conductance 
g(«) for a very large sheet. The details of the technique 
of measurement are standard, and more detailed results 
can be found in a published report [8]. The theoretical 
results shown are based on the calculations for the rib- 
bon which of course neglect the presence of the wave- 
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| 2 3 4 kd 5 6 1 8 
Fig. 4—The theoretical and experimental curves for the normalized 
shunt conductance of the half-wave longitudinal slot in the nar- 


row face of an X-band rectangular waveguide with a flange of 
width 2d. 


guide behind the plate. For larger plates the agreement 
is good, and for small plate widths the inadequacy of the 
theory is evident. 

The above computed data for slot conductances 
should prove useful in the design of slotted antenna 
arrays on cylindrical structures. For surfaces of more 
complex shape the results show at least the error in- 
volved in assuming the surface to be flat and infinite in 
extent. 
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A Method of Analyzing Coupled Antennas of 
Unequal Sizes* 
CJUASUEVIS? ANpa@ aah aan ul 


Summary—The impedance parameters of coupled antennas are 
calculated by a variational method. The unknown coefficients per- 
taining to the trial functions for the current distributions on various 
elements are shown to satisfy a set of linear equations involving the 
excitation conditions, but the final formulas for the impedance pa- 
rameters involve only the trial functions and the geometry. The 
method is directly applicable to the analysis of Yagi-Uda antennas, 
and the technique can be extended to coupled slots inside a wave- 
guide or coupled radiating slots excited by waveguides. 


INTRODUCTION 


IGNIFICANT contributions to the problem of 
S coupled cylindrical antennas of unequal sizes have 

been made by Uda and Mushiake! in their study 
of Yagi-Uda antennas by an extension of Hallén’s meth- 
od for a single antenna.? A similar method was used by 
King and Harrison* in their treatment of two coupled 
antennas of identical size. The method described in this 
paper is an extension of Storer’s variational method, 
originally used for the treatment of a single cylindrical 
antenna.*® An extension proposed by Begovich® gives 
variational formulas for coupled slots, provided they 
are small electrically. Although no numerical values are 
yet available for the solution of the Yagi-Uda antenna 
by this new method, we should like to present here a 
general derivation of the essential formulas involved in 
the method. 


SIMULTANEOUS INTEGRAL EQUATIONS AND 
THE VARIATIONAL EXPRESSION 


We will illustrate the method by considering two 
coupled antennas of unequal sizes. Corresponding for- 
mulas for x coupled antennas are given in the Appendix. 
Fig. 1 shows a sketch of the two coupled antennas. The 
approximate integral equations governing the current 
distributions on the two antennas are given by 


* Manuscript received by the PGAP, August 10, 1955. The re- 
search reported in this paper was sponsored in part by the Air Res. 
and Dev. Command, Wright Air Dev. Ctr., Dayton, Ohio. 

+ Antenna Lab., Dept. of Elect. Eng., Ohio State Univ., Colum- 
bus, Ohio. 

1S, Uda and Y. Mushiake, 
ing Co., Sendai, Japan; 1954. 

2 "Hallén, “Theoretical investigations into the transmitting and 
receiving qualities of antennae,” Nova Acta, Uppsala, ser. IV, vol. 11, 
November 11, 1938. 

3 R. King and C. W. Harrison, “Mutual and self-impedances for 
coupled antennas,” J. Appl. Phys., vol. 15, pp. 481-495; June, 1944. 

4 J. E. Storer, | Variational solution to the problem’ of the sym- 
metrical cylindrical antenna,” Tech. Report no. 101, Cruft Labora- 
tory, Harvard Univ., Cambridge, Mass.; 1950. 

Oe mipalarne variational solution to the problem of cylindrical 
antennas,” Tech. Report no. 12, Aircraft Radiation Systems Labora- 
tory, Stanford Research Institute, Stanford, Calif.; 1950. 

N. A. Begovich, “Waveguide excited slot radiators,” 
Memo no, 219, Hughes Aircraft Co. Res. and Dev. Labs., 


“Yagi-Uda Antenna,” Sasaki Publish- 


Tech. 
1950. 


1(2’)G13(2 a 2) dz" 


ly 
V,6(z — 0) = f Ih 
=p 


le 
= f I2(2')Gi2(z ia 2!) dz’ (1) 
=r 


h 
Vo d(z — 0) = f 
qty 


T1(z’)Gy2(2 — 2’)dz’ 


lg 
+f I2(2’)Gao(z — 2’)dz’, (2) 
ET 


where 


Vi, V2=voltages applied to the two antennas, 
I,(z), I2(z) =current distributions on the two antennas to 
be determined, 
5(g—0) = 


. VV ry2 2 
: jkn 1 e2 e ik (2—2")" +a 
VAs ) = (1+ A <) i aia ea ) 
4a Rk? as?/ \/(z—2/)?+a,? 


m Ufo 4 
= (=) = 1207 ohms, 
€ 


k=2r/d, 


a,=radius of the first antenna. 


Dirac delta function defined at z=0, 


Gilz 


Sees, 


——2: £, 


naar | 
init eas 


| 
—— 


Fig. 1—Two coupled antennas of unequal sizes. 


The functions Go2(z—2’) and Gi2(z—32’) are obtained by 
replacing a; in Gy(z—2’) by a2 and d respectively, where 
a, denotes the radius of the second antenna and d the 
distance between them. 
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From the linearity of Maxwell’s equations, it can be 
shown that the currents and voltages at the antenna 
terminals satisfy the relationships 


Vi = 21,0) + Zi212(0), (3) 
V2. = Z 121 (0) is ZI (0). (4) 
The parameters Zy1, Z12, Z22 characterize the impedance 


behavior of the antennas and are to be calculated. Eqs. 
(1) and (2) may therefore be rewritten as 


[Z1:71(0) a. Z 121 5(0) ]a(z = 0) 


ly ' 
= ih Th(3’)Gu(z == 2')dz' off an, \Gi0(z PAs "\dz', (5) 


1, 


age + Z 22] (0) Ja(z ae 0) 


= “ha! \G2(z ee z') \dz' ae fh o(2')Goo(s =a z')dz'. (6) 
—l, 
If we multiply (5) by Ji(z) and ener from —I, to i, 
multiply (6) by J2(z) and integrate from —/, to J, and 
add the results, we obtain 


Z1:11?(0) + 2Z1011(0) Io(0) ++ Zoo o?(0) 


= ff reuwn+ 2f [ reutet f [ 16.0 


where the specific arguments of the integrals at the right 
have been omitted for convenience; 2.e., 


lm la 
ff 1p CRY le = f iE PACE NG an 
sly YM ly 


Eq. (7) is the basis of a variational solution for the 
parameters Zy1, Z12, Z22. It may be used in the following 
way. We choose the currents corresponding to three ar- 
bitrary but nonidentical excitation conditions (denoted 
by single, double, and triple dots, respectively) and sub- 
stitute them in (7). The result is three simultaneous 
equations, viz., 


Waals (0) + 2Z1211(0)I2(0) at ATE) 


= if an PGub ee iH if Poti F f Teh) 


Z1 (0) a RAC )12(0) Nt Zend 520) 


=f fteri+2 ff toni+ ff ton, oo 


ZuT2(0) + 2Z211(0)T2(0) + Z22127(0) 


= if f NGuli + 2 fi i} LGils + if ii TsGool. (11) 
This system can be solved for the parameters Zu, Z12, 
Z», provided the determinant of their coefficients does 
not vanish, and this is precisely what was meant above 
by “nonidentical” excitation conditions. (It can be 
shown that it is necessary and sufficient to require that 
I,(0)/J.(0) have a distinct value for each of the three 
excitation conditions. In practice it may be convenient, 


z')I,(2')dz'dz. (8) 
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when two such distributions are known, to obtain the 
third by superposition.) The solutions of (9), (10), (11) 
will be the true values of Z, Z12, Zo2 provided the cur- 
rents are true currents; 7.e., currents which satisfy (1) 
and (2). The solutions have the further property of be- 
ing stationary with respect to variations of the currents 
about such true distributions. This property is demon- 
strated below. 

To demonstrate the variational nature of the simul- 
taneous solutions of (9), (10), (11), we take variations 
of the currents about their true values. If we make use 
of the symmetric nature of the Gna we get from (9) 


621 12(0) + 28Z121,(0)L3(0) + dZ 201 22(0) 
+ 26f,(0)|Zi:f1(0) + Zi212(0) | 


+ 2612(0) |Z12f1(0) + Zeof2(0) | = 2) {f if SiG 


i f i! siete if f ahcnee ii f slabs | (12) 


Since J, and I, were to be varied about the true distribu- 
tions, (5) and (6) must be satisfied by these distribu- 
tions. Multiplying (5) by 6/,(z) and integrating from 
—l, to l, gives 


614(0) [Z1111(0) +Z 121 2(0) | 


=f [iGatt ff shut. (13) 


A similar expression can be obtained by multiplying (6) 
by 6/.(z) and integrating from —/, to /,. When these re- 
sults are used in (12) we find 


6Z11112(0) + 26Z 121 1(0)L2(0) + 5Z22f.2(0) = 0. (14) 


The same procedure, applied to (10) and (11), gives re- 
lations identical to (14) except that the single dot is re- 
placed by a double or triple dot, respectively. These 
three equations may be solved simultaneously for 6211, 
6Z12, 6Zo2. Since the determinant of the coefficients is 
the same as that in (9) to (11), the only solution is 


6Z 11 = 6Z 12 — 6Z 99 = 0. (15) 


This guarantees that if approximate currents which dif- 
fer from the true currents by errors of the first order 
are used in (9) to (11), the calculated parameters will 
differ from the true values only by errors of the second 
order. 


APPROXIMATE SOLUTION USING [TRIAL 
FUNCTIONS OF Two TERMS 


Unfortunately it is generally difficult to guess good 
approximations for the currents Lz), In(z), etc., di- 
rectly. Instead it is customary to represent the currents 
as a sum of functions with coefficients to be determined 
by some auxiliary method. To discuss the final step of 
the formulation, which leads to formulas for numerical 
calculation, let us consider the approximate solution 
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which is based on the use of trial functions consisting of 
two terms. This is actually the case that we intend to 
apply to the analysis of Yagi-Uda antennas. We assume 


Ti(z) = Aifilz) + Asfe(z) (16) 
I2(z) ae A3f3(z) re Aafa(z), (17) 


where the f(z) are four known functions which are be- 
lieved to be a good representation of the current dis- 
tributions on the two antennas. For example, we can 
choose fi(z) and f;(z) to be sin[k&(4—|z|)] and 
sin [R(l.— | z| )], which are the distribution functions for 
“infinitely thin” antennas. If we use these functions 
alone, the variational solution would yield the same 
result as that obtained by the induced emf method.’ 
The terms f2(z) and f4(z) are therefore necessary in order 
to obtain an improvement over the old method. As in 
the case of the single antenna‘ we can use functions of 
the form cos (k) —cos (kz) or k(J—|2|) cos [k(7—|z|)] 
for fo(z) and f,(z). It should be noted that the coefficients 


A, will be, in general, functions of the excitation. The ~ 


currents corresponding to particular excitations will be 
characterized by particular values of the A;; for instance 
we may write 


T(z) = Aifi(z) + Aofa(s), (18) 


and so forth. 
When the currents are assumed to be of this form, the 
stationary expression, (7), becomes 


Zu |Aia oh Ara]? a4 2Z12[A 141 = Aas] [A sas 2 A4a4| 
+ Zo2[Asa3 + Asaa]? = De DA Arvin, (19) 
‘Sontag 


where 
ap 
Vik = i F(2)G@mn(z — 2’) f(z’) dz’dz = yx; (20) | 
SPO, 
with 
1 dor e192 
m = 1 ‘ (21) 
2 for 7 = 3,4 
ty forve A120 
ea (22) 
2 for k = 3, 4 


The constants 1, d:, a3, ad, are abbreviations for 
fi(O), fo(0), f3(0), and f4(0), respectively. Because of the 
stationary character of (7), as shown in (15), we shall 
determine the coefficients by differentiating (19) with 
respect to the A; and requiring 


OZ 14 OZ 19 OZ 25 
=——=——= (23) 
OA, OA; OA; 
The result is a set of equations linear in the A;; 
@1(Z11(A 101 + A 22) + Z12(A3a3 + Aa.) | 
= Ayyu ae Aoyie = Asyis Sie Asyia; (24) 


7P. S. Carter, “Circuit relations in radiating systems and ap- 
plications to antenna problems,” Proc. IRE, vol. 20, p. 1006; June, 
1932. 
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ae [Z11(A 101 ae A a2) a Z1(A 33 + A4as) | 


= Ayyis + Avyon + Asyos + Asya, (25) 
a3[Z12(A 141 + A2dz) + Z22(A 303 + Aads) | 

= Ayyis + Avyos + Asyas + Asya, (26) 
a4[Z12(A1a1 + Aode) + Z22(A sds + Asay) | 

= Ayyut Aoyos + Asysa + Asyus. (27) 


\ 
It should be remembered that the variational expres- 


sion, (7) or equivalently (19), was to be used to calculate 
Zu, Zy2, Zo by specifying three excitation conditions, 
thus obtaining three equations of the form of (19) with 
the A; replaced by Ax, Az, and Ax. Thus the set of linear 
equations above is in reality a shorthand for three sets; 
i.e., for solutions to (24) to (27) we seek three sets of con- 
stants Az, Ax, Ax, which together with a single set of Zn, 
Z, Zw will satisfy the equations. Any such solutions 
will also satisfy the variational expression (19), since 
this expression can be obtained from (24) to (27) by 
simple algebraic manipulations. 

To solve for a particular set of the Ax, we must spec- 
ify the excitation, and a convenient way of doing this 
is to set 


Z11(A ya, + Avade) + Z12(A3a3 + Agas) = Vi, 
Z12(A1d1 + Ande) + Z22(Asa3 + Asda) = Vo, 


(28) 
(29) 


and then to assign particular values Vi, VoeVi, Vue 
V2 to the constants V; and V2. When (28) and (29) are 
substituted into (24) to (27), we obtain the following 
solutions for the A;: 


Ax = [VilarAr, + a2Aon) + Vo(asAse + asAgn|/A, (30) 


where A denotes the determinant | vex and Ax the co- 
factor of the element yx. Eq. (30) again represents three 
sets of equations, the Ax, for example, being given by it 
with particular values Vi, V2 assigned to V; and V2. 
The final step is to substitute these values of A; into 
(28) and (29). From (28) we obtain, after collecting 
terms and noting that Ay, =A;, [because of (20) ], 


Vi [Zun(a17An + 2a,a2A12 + a27Ao9) 
ai Z12(d103A13 + d2d3A23 + aya4Aq4 + 204A 24) A| 
ae Vo|Z11(a1a3A13 + aya4Ai4 + adea3Ao3 + 20494) 


42 Z12(a3"A33 + 2a3asAz4 + ax?Aas) | = 0. (31) 


This equation can again be interpreted as a set which 
must hold when the particular values Vi, V2 etc. are 
substituted; this can be the case only if the coefficients 
of V; and V»2 vanish individually. Setting these coeffi- 
cients equal to zero, we obtain two equations in the un- 
knowns Zy, Z12 which have the solutions 


Z11 = A(a3"Ag3 + 2aza4A34 + ad4’Ays)/M, 
Zy.= 


(32) 
— A(aya3Ar3 + ayasArs + Goa3Ao3 + a204Ao4)/M, (33) 
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where 
ay7Ar + 2a,a2A12 
+ de’Ar. 
@103A13 + ayasAy4 
+ dea3Ao3 + aodsAog 


€103A33 + aya4Ayy 

+ a2d3Ao3 + ad2asAo4 
a3°A33 + 2asasAcy 

+ a47Agy 


Similar manipulation on (29) yields the same value for 
Zy, and in addition 


Zog = A(ayAyr + 2aya2Ai2 + 2"Ao2)/M. 


. (34) 


(35) 


These are the formulas from which calculations for the 
Yagi-Uda antenna are to be made. 

It should be noted that instead of the circuit relations 
of (3), (4) we could have written 


I,(0) — Vilu + V2V x2 (36) 

T(0) = Vil x2 + V2 20. (37) 

Parameters Yu, Vie, Yo2 are related to impedances by 
Zu Zia]? Vu Vie 

= , (38) 
Zi2 Zoe Vio Voo 


In our case, the Y’s have a simpler form than the Z’s, 
v1Z., 


Yu= (ay7Any + 2aya2Ay. + A2’Av2)/A (39) 
Vio = (@,03Ay3 + @yasAi4 + doa3Ao3 + d2d4Ao4)/A (40) 
(a3"As3 + 2azasAz4 + aa?Ags)/A. (41) 

Finally we may mention that the linear set of equa- 
tions for the A,, and the final results (39) to (41), can be 
obtained by applying V. H. Rumsey’s Reaction Prin- 
ciple®® to our problem. This method indeed leads to the 


admittance result more directly, but was not chosen 
here because the concepts are not as widely known. 


es 
I 


APPENDIX: EXTENSION TO ~ ANTENNAS 


The extension to 7 antennas with an arbitrary num- 
ber of trial functions for each antenna is straightfor- 
ward. The geometry of two typical elements of the array 
is shown in Fig. 2. Instead of going through the proof in 
detail, the corresponding steps are only indicated here, 
but the numbering is kept the same; e.g., (7a) here cor- 
responds to (7) in the two-antenna argument. The 
definition of G,; is obtained by replacing a; in Gy with 
a, for r=s, and with d,, for rs. In analogy to the two- 
antenna case, we now have 


ee Oey gt LG —.s’\do’) (12) 
tT —l, 
ez TO) (3a) 


8V. H. Rumsey, “The reaction concept in electromagnetic 
theory,” Phys. Rev., vol. 94, pp. 1483-1491; June, 1954 

9C. A. Levis and R. E. Webster, “Mutual Admittances and An- 
nular Slots,” Project Report 486-19, Antenna Laboratory, The Ohio 
State University Research Foundation; prepared under Contract 
AF 18(600)85, Air. Res. and Dev. Command, Wright Air Dey. Ctr., 
Wright-Patterson Air Force Base, Ohio; 1954. 
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n n 1 
D> Zel(O)8(2 — 0) = >> J Ti(2)Galz — 3’)dz’ (Sa) 


r=1¥ —l, 


n n 


we ZOO) 


r=) 3=1 
n n Ll, l, 
=> f 1, I,(2’)Ger(z — 2')I,(z)dz'dz (7a) 
r=1 s=1¥ —1,VU —ls 
subject to Z;,=Z,;. In analogy with(9) to (11), we here 


have to write $n(m+1) linearly independent equations. 


Element r Element 5 


——7:-¢ | 


| 
a 


Fig. 2—Typical elements of the m-antenna array. 


If the Z,, are obtained as solutions of these equations, it 
follows that 


bZar = 0, (15a) 


Saale ee ee 


Suppose we want to represent each current J; as a sum 
of trial functions f,, with coefficients A, to be deter- 
mined by some auxiliary method. We write, 

Nk 


T(z) = YD” Aafalz) 


Iam; 


(16a) 


with m,=1 and mz4:=2e+1. This numbers the func- 
tions consecutively: The functions fi, fo, - ++, fry con- 
tribute to’ Jy, 7the functions fay, faa) fap one 
tribute to Jz, etc. We also abbreviate by writing a; for 
fi(0). From (7a) we now obtain 


Dee Agar) Aon Le Aa ey 
7=1 s=1 q=m, pHMs, q=1 p=1 


where 


lq Ip 
VD ie fe OGoole — 2')fp(2')dz’dz = ypq (20a) 
EOL Ae 
with u and v chosen to that muSqSnmu, mSpSm. 

Again we obtain the coefficients by differentiating 
with respect to the A; and requiring 0Z,,/0A;=0 for all 
s, r, and k. Noting the symmetrical properties of the 
Zer and Yqp, we obtain 
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nr nn 


ak Di fag > A gq = Dd Aavaks Rit, 2,038: 
f=1 


q=™M, . q=l 


(24a) 


ln 


with s determined by m,.SkSn,. To specify excitation, 
it is convenient to let 


Si aeweed aes 


r=l q=™M, 


51,2 ere ie ee) 


and then to assign specific values to the V;. Use of (28a) 
in (24a) gives 


tn Resa UD aah, ton 
Piles Di A evar We ne 
This system can be solved for the A, 
1 n Ns 
A,= 5 SRS Hes (30a) 


é—1 k=ms 


where A is the determinant | Yor (assumed nonvanish- 
ing) and A, is the cofactor of the element yx. It is now 
necessary to insure that these values also satisfy (28a). 
Substitution in this equation gives 


n n 
ys v.( Zain — An) a OO se * iy (314) 
t=1 r=1 
where 
nr Tht 
At, = 25 Ag gk 
q=m, k=mt 


IRE TRANSACTIONS ON ANTENNAS AND PROPAGATION 


April 


and 


(erOreti—as: 
= 4 
()) Morr HES S. 


Eq. (31a) must hold for 4u(m+1) sets of nonidentical 
excitations. This is possible only if the coefficients of the 
V, vanish; 2.e., 


> $= 1,2,+--,0 
Lian = Abts x 
r=1 ee eerily Homo 640 


For fixed s, we can consider this a set of m equations in 
the n unknowns Z,, having the solutions 


(32a) 


where | aver denotes the determinant of the a,, (assumed 
nonzero) and A,, the cofactor of the element a, in the 
determinant. Again it turns out that an admittance rep- 
resentation yields a simpler formula; vzz., 


(Ghys nr Ns 
Ys = = ( Sos agndet) / A 
A q=m, k=mg, 
s = 1, ; 
(39a) 
ail ,n 


Radiation Characteristics of the Spherical 
Lunebere Lens* 
EDWARD H. BRAUN} 


metric region in which the dielectric constant 
varies as 2—r? (where 7S1 is the normalized 
radius) will focus rays from a point source at r=1 intoa 
plane phase front having a propagation vector pointing 
along the diameter passing through the source point. 
Construction of such a dielectric region has been im- 
practical in the past because of lack of availability of 
suitable dielectric materials. Recently, however, ma- 
terials have been developed which make it feasible to 
build such a lens. This so-called “Spherical Luneberg 
Lens” has a number of practical applications, based 
principally on its high degree of symmetry. These ap- 
plications as well as the details of construction however 
will not be described in this paper, the principal purpose 


UNEBERG [1] has shown that a spherically sym- 


* Manuscript received by PGAP, May 27, 1955; revised manu- 
script received November 29, 1955. 

+ Microwave Antennas and Components Branch, Electronics Di- 
vision, Naval Res. Lab., Washington, D. C. 


of which is to discuss the far field radiation character- 
istics of the lens. 

Initially we will assume only that ‘the amplitudes of 
the electric and magnetic fields are arbitrary known 
functions of position over the surface of the lens, while 
the phases are such as to result in a plane phase front 
perpendicular to the axis of the lens when the rays are 
extended out through the aperture. 

The pattern is given by the expression 


ERD 6) 2 ae oe 
4r R 
where 
FQ, 4) = Sf x H — (6X H)-RR 


6 de i See eee 
= /~ (p X E) X R | eistfas, 
bh 
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Here p is a vector from the center of the lens to a 
point on the surface, p is a unit vector in this direction, 
R is a unit vector in the direction of the field point (See 
Fig. 1). Taking the dot product with 6 and @ (unit vec- 
tors at the field point) we obtain, 


Fp = Sfle x H)-6 — /~ (p xB) | gare Eas ont) 


re~ [f[oxtr e+ 4/<@xBa]o%es. o 


We take the electric field at any point on the lens to 
be linearly polarized in the x direction and given by 


Cea 
E=E(8"," $’)&, where & and ¢’ are polar coordinates on 
the surface of the lens. 

Hence 


Es, = E&%-8' = E cos & cos ¢’ 
Ey = Ex-¢' = — Esin ¢’. 
The magnetic field at any point on the lens is given by 
H=H(s', ’)%, hence 
He = H3-8' = H cos & sin ¢ 


and 
Hy = Hy-¢' = H cos ¢’ 
ue ees ah 6 
p X H-6 = H cos@ sin ¢’¢’-60 — H cos $0": 8. 
Now: 
¢’-6 = cos 8 sin (@ — ¢’) 
3-6 = sind’ sin 6 + cos W cos 6 cos ( — ¢’) 
px E.é = Ecos¥’ cos ¢'¢'- + Esin 8’ 
¢’-® = cos (@ — ¢’) 
3’. ® = cos sin (¢’ — &) 
p x H-é = H cos sin ¢’ cos (® — ¢’) 


— H cos # cos ¢’ sin (¢’ — ®) 
oe 
(p X E)-0 = Ecos # cos ¢’ cos 6 sin (@ — ¢’) 
+ E sin ¢’[sin & sin 6 
+ cos 3 cos @ cos ( — ¢’)| 
ae 
p-R = acos# cos@ + asinw sin 6 cos (& — ¢’). 
The distance to a plane parallel to the (x, y) plane and 
passing through the point (0, 0, a) is d=a—p'3 
=a—a cos ?"’. 
Hence we may write / and H as 
p= E\(d", gp’ )ethae—ika cos y’ 
H = Hye", @ jer*se—1hs cos 3! 


where the phases of Eo(2’, ’) and Ho(8’, d’) are con- 
stant over the sphere. 
Substituting into (1), and using the identity: 


cos ¢’ cos (® — ¢’) — sin ¢’ sin (@ — ¢’) = cos ® 
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Fig. 1 


we obtain: 


r/2 Qn ee 
— citeat f i, | (a cos 6 + (/~ Fa) 
0 0 Le 


-cos 0 cos ® + Ay sin 6 sin # cos ‘| 


Fo = 


. @-tkal2 sin? 0/2 cos g’—sin 8 sin y’ cos (—-4')] 
- sin 3'dd'dd’. (3) 
Using (2) and the identity 
sin ¢’ cos (® — ¢’) + cos ¢’ sin (@ — ¢’) = sin ® 
we obtain: 


a2 Qn Te 
Fo = ater [ f {| a a {/~ Eo cos | 
0 0 Me 


€ 
-cos &@ sin @ + /— Ey sin 6 sin ¢’ sin ott 
be 
ev ikel2 sin” 6/2 cos y/—sin 6 sin 9’ cos (@—¢’)] 
-sin 8/dd'¢’. (4) 
Letting Fy=a’ei*f, 


m/2 rea 
fo = f ds’ {| cos # + (/~ Fa | 
0 Kb 


< eg See 
-cos ® cos # sin &e—7#42 sin” 6/2 cos 9” 


Qn 
‘ Ir eika sin @ sin 9’ cos wag 
0 


nr /2 
ie f dd! Hy sin 6 sin?’ e-ika? sin® 0/2 cos 9” 
0 


2a 
- {) cos g’etika sin @ sin gy’ cos (©-$ dg, (5) 
0 
The integral with respect to ¢’ in the first term is 
la 
f etka sin @ sin g’ cos (@-4)dg! = IrJo(Rkasin 6 sind’). 
0 


We now specialize to the planes =0 (the x, 2 plane) 
and ®=7/2 (the y, z plane). 
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The integral with respect to ¢’ in the second term of 
(5) is then: 


Plane =0 

if cos gp’ etka sin 0 sin 3’ cos ¢' dg’ = 2rpJ (ka sin @ sin wv) 
0 

Plane P=7/2 


2a 
‘hi cos og’ eke sin @ sin y’ sin odd’ 
0 


o/=29r 
Ss it eika sin # sin 9’2qy — (x = sind’). 
¢ 


T=) 


The latter integral shows, taken together with the fact 
that a factor of cos ® appears in the first term of (5), 
that He vanishes everywhere in the plane ®=7/2, as 
could have been predicted beforehand. 

In the plane 6=0, we have: 


w/2 aes 
fo(0, ® = 0) = an f dd! E cos 8 + {/~ fa | 
0 Lb 
-cos 0 sin f/e—ika2 sin’ 6/2 cos 9’ J)(ka sin 6 sind’) 
w/2 
+ sin 6277 il do! Ho 
0 
“sin? e-ia2 sin’ 6/2 cos oJ, (ka sin 6 sin 0’). 
— 
If we now assume that E and H at the aperture are 


related as in a plane wave, we have, 
letting Ge=fo/2r, 


0 
u = 2kasin®?—, and v= kasin@; 
G,(0, & = 0) 


a /2 
= (1 + cos @) f Hy cos & sin 3’e7i 8 8’ J (vy sin 3) dv’ 
0 


aw /2 
+ sin ai [ Hy sin? 3’ e—7 5 8 Ji (y sin 3’) dd’. (6) 
0 


Considering the first integral: 


a/2 
I; f Hy cos # sin de 8 9’ J, (y sin 8’) dd’ 


0 

1 v x 2 
— Hoe 1 — (=) x Jo(x)dx 
YJ 9 v 


where x = vsin &. 


Integrating by parts: 


Tait (=) a8) 
2 


v 


U a/2 
—j “f sin? 3’J,(v sin 0 )e~™ © 8 Hody’ 
a) 


1 m/2 
— —{ sin 3’ J1(v sin 3’ )e—™ ©5 8’ (dH/dd’) dd’. (7) 
v0 
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To obtain the first term in (6), (7) has to be multiplied 
by (1-+cos 6) =2 cos? 6/2. The second term in (7) be- 
comes 


w/2 
= —jsin of sin? 3/J,(v sin &)e—™ © 8 Hod’ 
0 


and this just cancels the last term in (6). Hence we 
have: 


J \(v) Tv 
G9, ® = 0) = (1 + cos 8) Hy (=) ; 
v 


7 
COW aa 


2 w/2 
{ sin 0/J,(v sin 8’) e7™ © 8 (dHo/dd’)dv’. (8) 
0 


ka 


We now have to introduce an assumption regarding 
the amplitude distribution of the aperture field. The 
geometrical optics field in the aperture will, in general, 
exhibit a singularity at the edge of the lens, since the 
radiation pattern of the feed must be multiplied by 
sec!/? 3’ to obtain the aperture field as a function of the 
central angle ?’. Therefore, unless the feed pattern falls 
to zero (and more rapidly than cos? 2’) at 90 degrees, 
the field intensity at 3’ =90 degrees in the aperture be- 
comes infinite. In practice, of course, this is not the case 
since the reasoning based on geometrical optics is not 
valid in such regions of rapidly varying field strength, 
and, in addition, losses in the dielectric would tend to 
reduce the field strength near the edge of the lens. 

A very flexible aperture distribution to choose is 
1+A cos? 3’. This means that the illumination taper is 
1+A, and the factor Ho(r/2) appearing in the first 
term of (8) is just unity. By assigning various values to 
p and A we have considerable control over the way in 
which the amplitude varies over the aperture. Thus, for 
example, by using negative values of A we have fields 
which increase towards the edge of the aperture; by 
using positive values of A we have fields which decrease 
towards the edge. Different values of » determine how 
peaked the distribution is. 

With this aperture distribution the integral in-(8) be- 
comes (apart from a factor) 


w/2 
Doze f sin? 2 cos?! 3’ J,(u sin 3’) e—™ 8 Hdd", 
0 


Expanding the exponential 


oe Ss ee sin? 3’J,(z sin 0”) cos"+?-1 9 dv’ 
n=0 T(n ae 1) 0 


The integral can be evaluated by Sonine’s formula: 


{2 
i J,(v sin 3) sint! 9 cos2et! odie! 
0 


af 2°T(a + 1) 


niet Fipapa@) 


with w= 1 a= 
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_— & (ju) 
oe ae ERE J (n+p) /241() 


n=0 Tn i 1) 


_ Substituting this into (8), and putting in the factors 
which were omitted previously, the final expression for 
the pattern is: 


ONT 


Go(8, ® = 0) = (1 + cos) ~ 


If we now return to (4), and proceed as before, using 


the fact that E and ah in the aperture are related by the 
factor ./u/e, we obtain exactly the same expression for 
Ga. Thus the patterns in the two principal planes are 
identical. 

As A-—0 (constant illumination) the pattern ap- 
proaches that of a circular aperture having the same 
radius as the sphere. That this must be so can be seen 
from the following argument. We consider a plane wave 
traveling in the z-direction, and excise a hemispherical 
volume (oriented as shown in Fig. 1) from the field. The 
contribution from the sources inside this volume to the 
field may be computed from the fictitious electric and 


Go = (1 + cos 8) Fito) 


RO db mad 


Avr ue 


ae 


ka Vv m=0 
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= 4 
magnetic currents 7 XH and —” XE flowing on the sur- 
face. Since there are no sources in the volume, this con- 
tribution must be zero. Hence the contribution from the 


r(m +=) 


T(2m + 1) 


1 2 pl2 2\ m 
I /2 m=0 v 


circular base (which corresponds to the circular aper- 
ture) must equal the contribution from the spherical cap 


J m1+(p/2)(2) 


; ead 
(» ae 


Sia p/2)(U). 9 
Tom 2D +(3/2)-+ (p/2) (0) (9) 


(which corresponds to the lens). 


u? Dire 
te =) (2m)! 


v(m oie *) J mst (p/2) (0) 


p 1 
(m+ 2+ =) Juramecorn() (10) 


PARTICULAR APERTURE DISTRIBUTIONS 


Let us now consider the explicit expressions obtained 
for specific values of p. For =1 we can transform (9) 
by means of the duplication formula for the Gamma 


function: 
G+) 
Det) Jr 


ih 1 nN 
(n + 1) i (F i 1) 
2 
which leads to the following expression for the pattern: 


uy? ™m 
“(=) sae J m+ (3/2) (2) 


if 


irre 2 (v) 


J (2) Al = 242 Tr 
Gp = (1 + cos A) . | 4/ J3;2(v) — s(t J 5)2(2) 
v ka u 20 DONG? 2v 


il a ae (v) 1 (=) 3 
SS —* —— v 6S Sse —— 
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For the range of @ of interest in practice, the series 
converges extremely rapidly [2]. For the main beam and 
first few side lobes it is only necessary to take the first 
two real terms, and none of the imaginary terms, except 
right in the vicinity of the nulls where one or two addi- 
tional terms in each series may be required for very high 
accuracy. 

Although the radius of the lens appears explicitly in 
the formulas, it can be eliminated for any lens likely to 
be of interest in practice. This makes it possible to 
write everything in terms of the single variable 
v=ka sin 6, which in turn makes it possible to derive 
general formulas for quantities which characterize the 
performance of the lens, independent of the lens size. 

To do this, we note that uw=ka— 1/(ka)?—v? may 
be expanded by the binomial theorem to read 
u=v /2ka+v'/6(ka)?+ ---. The second term is less 
than about 1 per cent of the first for values of ka 
and v which are of interest in practice (ka >~40, 
v<~10 for main beam and first side lobe.) Hence, 
to derive general formulas for beam width, first side 
lobe level, and aperture efficiency, we can write 
u=v"/2ka. 

For p=1, we have 


VAC), eee Zs 
Ge = +1 4/ = tanto | 
v v 2u 


BEAMWIDTH FORMULAS 


To find formulas for the beamwidth, we have to solve 
the equation 


Go(v) = 0.7071G4(0). 


At v=0 the value of Gs may be obtained by noting that 
J32 may be written in closed form as 


/2 /sin v 
ants) = 4/—( ; — cosa), 
Tv 


2J 1(v) 2A /sin v 
Ge(v) & | ( — cos r). 


v v v 


hence 


(11) 


Expanding sin v and cos v, the expression for Gg(v) in 
the neighborhood of v=0 becomes: 


9 


ar 1 v? v2 
Gis) = 142A pi to. 


ee 2A 
G0) =1-+ aes (11a) 


Therefore the transcendental equation to be solved 
for the beamwidths is: 


oT A 
Tao) el VERETO = 0.7071 (0. a <) Te Ges) 
v 
(The quantity W(m/2v)Jnyj2.(v) has been retained 


throughout since it is tabulated as Q,(v) in the WPA 
table of “Spherical Bessel Functions.”) 
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Solution of the transcendental equation for various 
values of A results in the following beamwidth formulas: 


Taper Half-Power Beamwidth 
— 20 db 1/2 = 25.7 — 
a 
Ne 
—10 db END 
a 
0 db 29.6— 
a 
NS 
10 db 310 
a 
NG 
20 db dees 
a 


(Positive values of taper refer to a decrease in field in- 
tensity at the edge of the aperture.) 
SIDE LOBE LEVEL 


To compute the side lobe level as a function of illumi- 
nation taper, we first differentiate (11) with respect to 
v, and set the result equal to zero: 


d (_Jy,(v) 2A /sin v 
an +—( — cos )} =" ()); 
dv v v v 


This results in the expression: 


407J o(v) 
A Fe . . 
sin Vv vsin v 
cos v — — 
v 3 


We would now like to solve this equation for v as a 
function of A, and calculate the side lobe position for 
various values of A. Then we could substitute the side 
lobe position in (11) to obtain the side lobe level. How- 
ever, since we cannot solve for v as a function of A, we 
consider the two equations 


4y°J o(v 
ee TEL 
sin v v sin v (13) 
cos v — : 
v 3 
and 
2A 
— 
3 
L = — 20lo 


12 = 
2J1(v) 2A Tr 
rt LW ee 


as a pair of parametric equations. By substituting dif- 
ferent values of v into these equations, we obtain corre- 
sponding values of taper (1+-A) and side lobe level (ZL). 
(The equation for Z has been normalized to zero at the 
peak of the main beam.) 

Values of v from about 4.6 up to about 6 will give all 
tapers of practical interest. This information on the ap- 
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proximate range of v may be determined from a plot of 
the field (Fig. 2) for no taper, which shows that v must 
be of the order of 5.2; v is then increased from this 
point until the denominator in the expression for A 
changes sign, since the point where it passes through 
zero corresponds to infinite taper, 7.e., zero field at the 
edge of the lens. It is then also decreased until 4 = —1, 
since this includes all cases where the field increases in 
going from the center of the aperture to the edge. 
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Fig. 3 


The curve of side lobe level versus illumination taper 
obtained from these parametric equations is given in 
Fig. 3. 

It is also of interest to examine the distribution 
cos #’, where the field drops to zero at the edge of the 
aperture. The position of the first side lobe is obtained 
from the condition that the denominator of (13) vanish. 
This leads to the transcendental equation 


—3v 


vy — 3 


tanv = 


which has the solution v=5.764. The side lobe level is 
then obtained from the second of (13) by substituting 
A=o, and v=5.764 
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3 
L = — 20 log | ————_— | = 


Tv 
7 J 3/2(2) 


If we now consider the field produced by the distribu- 
tion 1+A cos? #?’, which is a more peaked distribution, 
we have to substitute p=2 in (10). Again we find it 
necessary to retain only the first two terms in the real 
part of the series, and we can neglect the imaginary part 
of the series. The expression for the field in this case is 


2Ji(v) 4A 
Ga= Fic a J2(2). 


v/) 


To evaluate this at v=0, one can apply |’Hospital’s rule 
twice, or else use the series expansion for J,(v) resulting. 
in 
A 
G0) =1+ a (13a) 
To determine the beamwidth, we have to solve the 
equation 


2I(v) 4A A 
v 


which can be written in terms of the A functions tabu- 
lated in Jahnke and Emde: 


A A 
2 2 
For a 10 db taper (A =2.17) the root of this equation 
is approximately 1.8. This leads to the formula for the 


half-power beamwidth 


° 


Xr 
61;2 = 32.8 — (10 db taper, 1 + A cos? #’ dist.). 
a 


The side lobe level vs illumination taper is again de- 
termined in the same way by a pair of parametric equa- 
tions: 


2F1(v) — vJo(v) 


AJ 9(v) 
Fi) + —— = Fs(0) 
ees 
of ae =) 
2J,(2) + 44 — 


v 


A plot of these equations is given in Fig. 3. More 
peaked distributions may be discussed by substituting 
higher values of p in (10). 


GAIN 


The gain of the lens may be computed from the ex- 
pression 
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(14) 


Gain = 


1 € 
te — Ry/— | E(0) |? 
2 Ml 
Pr 


where | | in the numerator is evaluated at 0@=0, and 
Pr is the total power radiated. Pr may be computed by 
integrating the Poynting vector over the aperture of 
the lens; the result is 


e Dain y/ = lieestee shea 
amas are ae PUL) tae et 


Using the expression for E(0): 


—jop e~ihR 


EO) = aei**InGo(0), 


1 


the numerator of (14) becomes 


1 he Twp? 
Ar y/ | 
2 mm 4 
Combining this with the expression for Pr, and rear- 
ranging: 


o| Gua) |] 


| Go(0) |? 

4A Wee 
deena 
Thus the efficiency of the lens, 7, compared to the 


equivalent uniformly illuminated circular aperture, is 


| Ge(0) |? 


4a 
Gain = aa (1a?) 


t- 


NpA = 


For p=1, .e., for the aperture distribution 1+A cos ¥’, 
we have from (ila) that | Gs(0) | =1+2A. Hence 


(1+ 34)" 


ioe ay Bo 
By hr 


This is plotted in Fig. 4. For p=2, | Ge(0)| =1+(4/2) 
from (13a). Hence 
A 2 
Al =, 
Ga 


LAS eae ay od. SMe 


oe ae 
DA te 


UA oa 


This is also plotted in Fig. 4. 
For the distribution cos 2’, we can obtain the effi- 
ciency from ma by allowing A to become infinite: 
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For the distribution cos? J’, we allow A to become 
infinite in ne: 


3 0.75 
N20 = he . . 


Again, efficiency with more peaked distributions may 
be calculated by proceeding in the same way from (10). 
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Fig. 4 


The spherical Luneberg lens has a lower sidelobe level 
and a somewhat wider beamwidth than the correspond- 
ing two dimensional lens having the same illumination 
taper. This is due principally to the effective illumina- 
tion taper produced by the circular aperture shape. 


CONCLUSION 


Although the precise aperture distribution which will 
be obtained with a given primary feed pattern cannot 
be predicted exactly on the basis of geometrical optics, 
these results indicate that satisfactory far field patterns ~ 
are obtainable with almost any reasonable amplitude 
taper. In addition they provide information as to the 
magnitude of the changes required in the amplitude dis- 
tribution to produce significant improvements in side 
lobe level, beam width, and gain, and make it possible 
to arrive at a compromise distribution which will result 
in satisfactory values for each of these three parameters. 
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Antenna Pattern Distortion by Dielectric Sheets’ 
J. H. RICHMONDt 


Summary—The radiation pattern of an antenna covered with a 
radome can be determined from the following data: The antenna 
Fresnel-Zone fields in the absence of the radome, and the fields 
inside the radome produced by a plane-wave incident from an ex- 
ternal source. 

In certain cases the pattern calculations may be simplified with- 
out excessive loss of accuracy by using optical theory to describe 
approximately the plane-wave transmission through the radome, and 
neglecting scattering by the antenna in determining the fields inside 
the radome from the external plane-wave source. 


INTRODUCTION 


N DESIGNING a radome it is often necessary to 
] calculate the far-field pattern of an antenna 

covered with the radome. Such calculations are gen- 
erally based on given data on the Fresnel fields of the 
antenna transmitting without a radome. Ray-tracing 
methods are used to describe approximately the trans- 
mission of these fields through the dielectric radome. 
It will be shown that the problem can be reformulated 
in a manner which requires a knowledge of transmission 
through the radome only for a simple plane wave, rather 
than for the complicated Fresnel-Zone waves of the 
antenna. 


THEORY 


Consider the problem of determining the far-field pat- 
tern of antenna 1 shown in Fig. 1. Let £; and HM; denote 
the fields of antenna 1 when transmitting. Choose a 
reference surface S which must be either an infinite 
plane or any surface which completely encloses antenna 
1. In computing the fields of antenna 1 on the source- 
free side of S, antenna 1 can be replaced by its equiva- 
lent electric and magnetic current distributions J; and 
K,; on S, where 


Jy my xX Ay, (1a) 


and 


Ki =hXm, (1b) 


where 7; is a unit vector in the direction of the outward- 
drawn normal on S. Two vector potentials A and F can 
now be defined in terms of J; and K, in the usual man- 
ner,! and the field intensity can be expressed in terms of 
A and F. If the entire region on the source-free side of S 
is occupied by a homogeneous medium, the potentials 
and the electric field intensity in this region are given 
by: 


* Manuscript received by the PGAP, January 2, 1956. The re- 
search reported in this paper was sponsored by the Air Res. and Dev. 
Com., Wright Air Dev. Center, Wright-Patterson AFB, Ohio. 

+ Antenna Lab., Dept. Elec. Engrng., Ohio State Univ., Colum- 
bus, Ohio. ie 

1A. B. Bronwell and R. E. Beam, “Theory and Application of 
Microwaves,” 1st ed., McGraw-Hill Book Co., New York, N. Y., 
pp. 451-453; 1947. 
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eikr 
Eek f f Te ee (2a) 

Air s iif 

€ Cain 
Dee —ffx ds, (2b) 

dq Ss ip 

and 
1 j 

E=—-—VXF-—VxXVXA. (2c) 


€ Lave 


Fig. 1—Coordinates used in radiation pattern analysis. 


The expressions for A and F become much more compli- 
cated for the case of an inhomogeneous region. There- 
fore, radome calculations are usually divided into two 
steps. First, an auxiliary reference surface S’’ is chosen 
between antenna 1 and the radome. The unperturbed 
fields of antenna 1 on S’’ are used to obtain an approxi- 
mation for the fields on surface S which is outside the 
radome. Eqs. 2 are then used to compute the far-field 
pattern. This type of analysis, which may be called the 
conventional approach to the radome problem at pres- 
ent, requires a mathematical description of the trans- 
mission of energy from antenna 1 through the radome 
to the external reference surface. 

The problem may be reformulated by introducing a 
distant transmitting antenna, antenna 2, as in Fig. 1. 
The far-field pattern of antenna 1 is given by the voltage 
V induced at its terminals as antenna 2 is moved in a 
circular path centered at antenna 1. Antenna 2 is always 
aimed at antenna 1 during this motion. Let E, and H» 
denote the electric and magnetic field intensities on S 
due to antenna 2 transmitting im the presence of antenna 
1. Antenna 2 can be replaced by its equivalent electric 
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and magnetic current distributions Jz and Ke on S, 
where 


= 22 & As, (3a) 


and 


Ks == Es x N2, (3b) 
where 7,2 is the unit vector normal to S in the direction 
of the region occupied by antenna 1. By an application 
of the Lorentz reciprocal theorem, Kouyoumjian? has 
shown that 


V « Jf ore Sky Hes (4) 


As will be shown, the additional information introduced 
by the application of the reciprocity theorem makes it 
possible to simplify the radome analysis by strategic 
use of (4) in place of the conventional approach. 

If S is chosen to be a plane, (4) simplifies to 


V« Jf te206 (5) 


For far-field calculations, antenna 2 is at a great dis- 
tance from antenna 1, the field on S due to antenna 2 is 
a plane wave (if scattering from antenna 1 may be ne- 
glected), and (5) reduces to the usual expression for the 
far-field pattern of antenna 1 in terms of its aperture 
fields. 

To illustrate the application of (4) to the radome 
problem, consider the far-field pattern of antenna 1 in 
the presence of an infinite plane dielectric sheet as in 
Fig. 2. Choose as a reference the surface of the dielectric 
sheet facing antenna 1. Now to use (4) or (5) it is neces- 
sary to know: the fields on S of antenna 2 transmitting 
through the dielectric sheet in the presence of antenna 1, 
and the unperturbed fields on S of antenna 1 when it is 
transmitting. Neglecting the scattered fields of antenna 
1, the fields of antenna 2 on S are found by employing 
the plane-wave plane-sheet transmission coefficient T. 
Using primes to denote quantities in the presence of the 
dielectric and the absence of primes to denote the corre- 
sponding quantities without the dielectric, 


H.’ =~ THe, (6a) 


and 


Saf y IS. (6b) 


V'« iat Jo! - Eyds. (7) 
S 


2R. G. Kouyoumjian, “The Calculation of the Echo Areas of 
Perfectly Conducting Objects by the Variational Method,” Report 
444-13, Contract DA 36-039 sc 5506, Ohio State University Res. 
Found., pp. 29 and 80-84; Noy. 15, 1953. 


From (5), 
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DIELECTRIC 
SRE 


Fig. 2—Transmission through an infinite plane dielectric sheet. 


From (6) and (7), since TJ is independent of position 


on 5S, 
V! = rf [ t.Bas a Yee 
Ss 


V(b) ~ T()V(9). 


That is, knowing the unperturbed far-field pattern V of 
antenna 1, its pattern V’ in the presence of an infinite 
plane dielectric sheet is found merely by multiplication 
by the plane-wave plane-sheet transmission coefficient 
for the polarization and angle of arrival of the plane 
wave coming from antenna 2. This contrasts with the 
usual optical approximations in which the transmission 
coefficient for the angle of incidence y (Fig. 2) of rays 
from antenna 1 is employed. 

Suppose the dielectric surface facing antenna 2 had 
been chosen as the reference surface. The dielectric 
sheet would then be considered a part of antenna 1 and 
it would be necessary to know the fields of antenna 1 on 
this surface in the presence of the dielectric sheet. Thus, 
a surface passing between the radome and antenna 1toG 
perhaps coinciding with the dielectric surface facing 
antenna 1, is preferred since transmission through a 
dielectric body is most easily calculated for plane waves. 

Eqs. (8) are approximate in that they neglect the 
waves scattered by antenna 1. To determine the error 
involved in this approximation, several examples of 
transmission through infinite plane sheets and finite 
plane sheets will now be considered. 


(8a) 


or 


(8b) 


EXAMPLES 


Fig. 3 shows the measured far-field H plane patterns 
of a pyramidal horn at 9,400 mc with a polystyrene 
(€= 2.53) plane sheet held at an angle of 40° across the 
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Fig. 3—Far-field pattern of horn trahsmitting thrausmitting 
a plane polystyrene sheet. 


horn aperture. The unperturbed pattern is also shown 
for reference. The calculated points shown were ob- 
tained using (8b). The dielectric sheet was 0.51 inch 
thick and was large enough so that its edges were well 
outside the region illuminated by the horn. The horn 
had a length of 7.8 inches, flare angles of 20°, and a 
beam width of 23° to the half-power points. The waves 
arriving from antenna 2 were polarized perpendicular to 
the plane of incidence on the dielectric sheet. Hence, the 
transmission coefficient for perpendicular polarization 
was employed in the calculations. 
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Fig. 4—Far-field pattern of paraboloid transmitting 
through a plane polystyrene sheet. 


Fig. 4 shows similar results for a paraboloidal antenna 
with the same polystyrene sheet at 45°. The paraboloid 
had a diameter of 4 inches and a beam width of 23°. 

One more example using the horn is shown in Fig. 5 
for a fiberglas sheet having a thickness of 0.382 inch 
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and a dielectric constant of 4. This sheet had a lower 
transmission coefficient than the polystyrene sheet and, 
therefore, produced a greater reduction in the far-field 
amplitude. The fiberglas and polystyrene sheets were 
each a little over one-half wavelength thick in terms of 
the wavelength in the dielectric. 

To get closer to the radome problem, calculations 
were also carried out for plane dielectric sheets covering 
only a part of the antenna aperture. A reference surface 
S was chosen on an infinite plane coinciding with the 
dielectric surface facing antenna 1. The fields of an- 
tenna 1 were measured on S in the absence of the dielec- 
tric sheet. Eq. (4) or (5) requires a knowledge of the 
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Fig. 5—Far-field pattern of horn transmitting 
through a plane fiberglas sheet. 


fields on S due to antenna 2 radiating in the presence of 
the dielectric. Although at points remote from the edge 
of the sheet the fields are given by plane-wave plane- 
sheet theory, no solution appears to be available for the 
fields near the edge. As a first approximation, the fields 
were assumed to be the incident plane wave up to the 
edge of the dielectric sheet, and were assumed to be 
given by plane-wave plane-sheet theory on the portion 
of S covered with dielectric. The results for this approxi- 
mate calculation, using (5), are shown in Figs. 6 and 7 
for a polystyrene sheet 0.51 inch thick covering one-half 
of the horn aperture at angles of 40° and 0°, respec- 
tively. 


CONCLUSION 


Knowing the unperturbed fields of an antenna, it is 
often desired to calculate the radiation pattern of the 
antenna covered with a radome. Generally, optical ray- 
tracing methods are employed to describe approxi- 
mately the transmission of the antenna Fresnel-Zone 
fields outward through the radome. By applying the 
reciprocity theorem, however, this problem can be 
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Fig. 6—Far-field pattern of horn with aperture partially 
covered with a plane polystyrene sheet. 


traded for that of computing the fields inside the ra- 
dome for a plane wave arriving from an external source. 
This problem is simplified if scattering by the antenna 
can be neglected. To determine the error involved in 
neglecting this scattering, measurements and calcula- 
tions were made of the patterns of small horn and parab- 
oloidal antennas covered with various infinite and finite 
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Fig. 7—Far-field pattern of horn with aperture partially 
covered with a plane polystyrene sheet. 


plane dielectric sheets tilted at various angles from the 
antenna aperture. For the examples which were con- 
sidered, the scattering error was quite negligible and 
satisfactory results were obtained by using an optical 
approximation to describe plane-wave transmission 
through the finite plane sheets. 


Multiple Scattering by Randomly Distributed Obstacles 
—Methods of Solution” 


C..M.-CHUF Anp S.-W. GHURGHIERT 


Summary—Methods of solution of problems of multiple scatter- 
ing of electromagnetic radiation by randomly distributed obstacles 
are evaluated with respect to generality, accuracy, and suitability for 
numerical calculations. Exact solutions are found to be limited to 
infinite regions and are difficult to apply for anisotropic scattering. 
Approximate solutions of the transport equation are limited to simple 
regions and are very complicated for anisotropic scattering. Approxi- 
mate representation of the scattering process with the classical dif- 
fusion equation permits solution for almost all conditions, but the 
solution is inaccurate near the boundaries of the region of scattering. 
Approximate representation of the intensity by a number of discrete 
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components permits solution for many geometries. Several two-com- 
ponent representations have been developed which yield more ac- 
curate results than diffusion theory for anisotropic scattering. A 
six-component model appears to be considerably more accurate for 
most geometries but requires machine computations. 


INTRODUCTION 


HIS PAPER presents the results of an investi- 
Oa seer of methods of solution of problems of 
multiple scattering of electromagnetic radiation 

by randomly distributed obstacles. Proposed methods of 
solution are evaluated with respect to generality and to 
practicality for numerical calculations. Exact and ap- 
proximate methods are both considered. Attention has 
been given to obstacles of uniform size, but results can 
be extended to dispersions of known size distribution. 
Exact solutions have been derived for the scattering 
of radiation by a single obstacle of simple geometry such 
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as a sphere, ellipsoid, or cylindrical rod. For more com- 
plicated shapes, the amount and angular distribution of 
scattered power must be found by approximate methods 
or experimental means. For limited paths of transmis- 
sion such that rescattering is negligible, the transmission 
and angular distribution of scattered power can be de- 
scribed by an exponential decay law. For sufficient 
lengths of transmission, secondary and higher orders of 
scattering become appreciable and description of the 
process is much more difficult. 

The importance of multiple scattering has been recog- 
nized in the transmission of neutrons through dense 
material, of stellar radiation through dispersed matter, 
of light through fogs, of heat through fibrous insulation 
and microwaves through the atmosphere. Previous work 
has consequently appeared over a wide range of the 
scientific literature and has had diverse objectives. 
Most attention in the literature has been given to iso- 
tropic scattering within infinite or semi-infinite regions 
because such conditions are more receptive to mathe- 
matical analysis. In actual applications, the angular dis- 
tribution of scattered power is generally quite aniso- 
tropic and in many practical problems the assumption 
of an infinite medium is not justifiable. Principal atten- 
tion in this study is accordingly given to methods suit- 
able for anisotropic scattering and finite regions. 


MATHEMATICAL FORMULATION 


Single scattering can be characterized by the follow- 
ing parameters: 


1) os, the effective cross section of an obstacle for 
scattering; 

2) o:, the effective cross section of an obstacle for 
‘interaction; and 


> 
3) fi(Q, Q’), the angular distribution function for 
single scattering, which gives the fraction of energy 
— 


scattered from a beam with the direction, Q, intoa 


—> 
unit steradian in the direction, 0’. 


The cross sections and angular distribution function are 
functions of the index of refraction, the dimensions of 
the obstacle and the wavelength. For spherical obstacles 
the angular distribution is symmetrical about the direc- 
tion of the incident radiation and may be expressed in 
terms of a single angle, 6, measured from the incident 
direction; 7.é., AQ, Q’) reduces to fi(6). The combina- 
tions, dg =0;—<s, called the cross section for absorption, 
and w=a;/o;, called the albedo for single scattering, are 
often used for convenience. 


Multiple scattering can be characterized by the above 


di(R—O) 
dSa 
net rate of increase 
in intensity in the 


— 4 a) 
direction Q direction Q 
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parameters for single scattering together with NV, the 
number of obstacles per unit volume. The combinations, 
1/No;, 1/Noa, and 1/No;, called the mean free paths for 
scattering, absorption, and interception, respectively are 
also commonly employed. 

The dependent variable in multiple scattering is the 


— 
specific intensity, 7(R, 2). The specific intensity is de- 
fined as the energy intercepted per unit time by a unit 
—_ 


>) 
area at R, per unit steradian from the direction Q nor- 
—. 
mal to the area, where R is the radial vector from the 
—_ 
origin of the coordinate system and (2 is a unit vector 


> 
representing direction. In general, 7(R, Q) is a function 
of five variables, three for position and two for direc- 
tion. 
Multiple scattering can be formulated from either of 
two points of view: 


1) In the random-walk or Lagrangian point of view, 
the energy in a wave is considered to be in the form 
of discrete e packets (quanta). The specific inten- 


sity, ie. Q) is then interpreted as the density of 


these packets at R traveling in the direction 2 
The mechanism of single scattering can be con- 
sidered as the collision of an energy packet with an 
obstacle. The probability of a packet traveling in 


the direction, 2 ,_ being scattered into dQ’ by a 


collision is AQ, Q’ )dQ’, while the probability of a 
packet being scattered between x and x+dx is 
No,e-N“*dx, where x is the distance traveled into 
the dispersion or after a collision From these two 
elementary probabilities, it is possible to derive an 
expression for the probability of finding an energy 


packet in a unit volume at R traveling in a unit 


— 
steradian in the direction, Q, after one, two..., 
or any number of collisions, and | hence an expres- 


sion for the specific intensity i(R, 0). Exact solu- 
tions have been developed only for an infinite re- 
gion of scattering, but approximate solutions are 
possible for any conditions by the Monte Carlo 
technique. 

2) In the continuous or Eulerian point of view, the 
scattering region is considered to be a macroscopi- 
cally homogeneous media. An energy balance over 
a differential volume of unit cross section and 


length dSg along Q yields the following integro- 
differential equation, generally known as the trans- 
port equation: 


==> 
Q 


\d — Noi(R,2. (1) 


rate of decrease of in- 


tensity in the direction 


a, 2 5 
Q due to interception 
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The integral term in 1 thus introduces the con- 
tribution of multiple scattering. Exact solutions of 
the transport equation by the Laplace transform 
and Fourier transforms have been postulated. Ap- 
proximate solutions have been proposed by quad- 
rature, spherical harmonics, and reiteration. 


Numerous approximate representations have also 
been proposed for the multiple-scattering problem. The 
transport equation can be approximated by the classical 
equation of diffusion. Representation of the specific in- 
tensity by a number of discrete components also leads 
to approximate equations which can be solved more 
readily than the transport equation. 


SOLUTION BY RANDOM WALK 
Exact Solutions 


Exact solutions for an infinite region of scattering 
have been developed by Grosjean [1, 2] based on the 
Lagrangian point of view. The results are expressed in 
integral form. As an example, in the simple case of a 
point source in an infinite, isotropically scattering region 
the probability of finding an energy packet between 7 
and r+dr from the source is $(7)dr with 

a 
do(r) = p(r) + — “[-= — sin rv dv (2) 


where 
p(x)dx = 


is the probability of an energy packet being scattered at 
a distance between x and «+dx from a previous scatter- 
ing, and 


Noe-Nt# dx (3) 


eh © p(x) sin (xv)dx 4) 
0 


x 
Monte Carlo Method 


The increasing availability of high-speed computing 
machinery and current interest in neutron scattering 
has stimulated interest in the so-called Monte Carlo 
technique of numerical analysis [3]. This technique 
might be applied to the multiple-scattering problem as 
follows. The region can be divided into a number of con- 
tiguous cells of any number of shapes. Based upon the 
characteristic properties of the scattering region, 1/No,, 


1/No:, and fi(Q, ©’), probabilities are assigned to each 
cell for interception, absorption, and rescattering out of 
each particular face of the cell for the power entering 
each face. The random walk process is then carried out 
for a very large number of packets of energy entering 
the dispersion, until a statistical pattern is apparent. 
The accuracy of the procedure depends on the number 
of cells and the number of quanta traced. In practice, 
the number of cells and the number of discrete directions 
in which scattering is allowed are limited by the storage 
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capacity of a computer, and economically, by the com- 
puting time required to obtain a statistically representa- 
tive result. Since the solution is numerical instead of 
analytical, the calculations must be repeated for every 
set of conditions. 


Exact SOLUTIONS OF THE TRANSPORT EQUATION 


Exact solutions of the transport equation have been 
worked out for only a few special cases, in all cases for 
infinite or semi-infinite regions and in mast cases for iso- 
tropic scattering. The solutions are usually expressed in 
terms of definite integrals from which numerical results 
may be obtained by standard techniques, such as residue 
theory, series expansion, or numerical reiteration. In 
principle, these solutions can be utilized to construct 
solutions for finite dispersions and anisotropic scattering 
but the numerical procedures become too complicated 
for practical usage. The techniques used in two exact 
methods of solution are outlined below. The details can 
be found in the original references. 


Laplace Transform Method 


The Milne problem in astrophysics; t.e., the determi- 
nation of the radiation escaping from a star surrounded 
by an infinitely thick, isotropically scattering media can 
be solved by application of the Laplace Transform. The 
solution is given in terms of an “H” function which is 
proportional to the intensity in any direction and is de- 
fined by the equation 


a) = +> Sawn f - a as, (5) 


where » =cos 8. The method is discussed by Chandrasek- 
har [4] who calculated and tabulated values of the “H” 
function. 


Fourier Transform Method 


For isotropic scattering in an infinite medium involv- 
ing two degrees of freedom, one for the space coordinate 
and one for direction, the Fourier transform of the trans- 
port equation yields a solution for the specific intensity 
in the form of infinite integrals. For example, for an 
isotropically scattering, infinite dispersion with a unit 
plane source at Z=0, the integral solution is 


; 1 2° A(v)er% 
u(Z, 9) > i ew (6) 
2a J. Nor + jou 
where 
1 
A(v) = : (7) 
No, WN. j 
fs eae 
270, Noy — jv 


Details are given by Glasstone and Edlund [5] and ad- 
ditional examples by Case, et al. [6]. 
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APPROXIMATE SOLUTIONS OF THE 
TRANSPORT EQUATION 


Quadrature Method 


Quadrature has been used extensively to obtain ap- 
proximate solutions of the transport equation in prob- 
lems of astrophysics and neutron scattering [7]. The 
applications have been limited to simple geometries and 
isotropic scattering. The principle involved in this 
method is to approximate the integral in (1) by a finite 
summation, thus reducing the integral equation into a 
corresponding, finite number of simultaneous differen- 
tial equations. The utility of the quadrature method 
therefore depends on the accuracy with which the in- 
tegral can be represented by a limited summation and 
by the complexity and number of the resulting differen- 
tial equations. In so far as actual calculations are con- 
cerned the feasibility of the method has been demon- 
strated only for parallel-plane symmetry. 

For parallel-plane symmetry, (1) reduces to 


di(Z, 
P a = — Nori(Z, p) 
+1 
4 InNo. f Teaver grin (8) 
= 


where Z =distance measured in incident direction. 
The integral can then be replaced by a quadrature 

formula such as that of Gauss, yielding a system of 2 

simultaneous linear differential equations: 


di(Z, Mk) 5 
ees = Noii(Z, Lx) 
+ 2rNo, Ds Cy(Z, wi)filurs os) (9) 
j=1 
with 
k=1,2---n 


The factors C; and y; are defined by the quadrature for- 
mula. The solution of this set of equations gives the in- 
tensity, 7(Z, ux), 7.e., the intensity for any position, Z, 
and for discrete directions. For a highly anisotropic 
distribution function, fi(u, uw’), the number of terms, x, 
and hence the number of equations must be very large 
in order to obtain a satisfactory solution. 

Quadrature appears to yield nonlinear differential 
equations for all other geometries. 


Spherical Harmonics Method 


A standard method of solving the transport equation 
is by the use of spherical harmonics [8]. The specific 
=> 


intensity 2(R, Q) is developed as a series of spherical 
—)> 


harmonics in Q, such as 


3 (2m + 1) Em(R)¥n(Q) (10) 


m=0 


> —> 
(Rk, 2) = 
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where Nios (2) are spherical harmonics of order m and 


E,,(R) are functions of the space coordinates R only. 

Using the orthogonality properties of spherical har- 
monics, it is possible to reduce the integro-differential 
equation into a system of simultaneous ordinary differ- 


ential equations in E,(R). For parallel plane symmetry, 
oma 


Yn(Q) degenerates into the Legendre polynomials, the 
component equations are linear, and a general solution 
may be obtained. With spherical symmetry and appar- 
ently any other geometry the method yields nonlinear 
differential equations of such a complicated form that 
no simple solution is apparent. 

An infinite set of algebraic equations can be obtained 


by taking moments of the space function F,,(R). Ap- 
proximate methods have been proposed to obtain a solu- 
tion from this infinite set [8]. 

It has been pointed out by Chandrasekhar that the 


approximation of 7(R, @) by a series of spherical har- 
monics of 2m terms is equivalent to the method of 
quadrature of order n=2m. The difficulties and results 
of the method of spherical harmonics are therefore the 
same as those for quadrature. 


Reiterative Method 


A formal solution of the transport equation may be 
expressed by a pair of integrals, which can, at least theo- 
retically, yield numerical results by reiteration. The 
necessary order of reiteration, and the complexity of the 
integrals, however, limit the usefulness of such a solu- 
tion. 


Briefly, if a “source function” is defined as, 
—- a a a es 
GUR, 2), = i i(R, 2’) f1(Q, 2) dQ’ (11) 
Qe 
the transport equation may be reduced to 
os 
di(R, Q) —> aed 
See Ne lR pata a RAO) (12) 
Q 
Direct integration of (12) yields 
= =: ae 
i(R, Q) = vo. | a(R, Q)e-NetlR’-RId So (13) 


in which the integration is carried out, along a straight 


line in the direction @ from the boundary to the point R. 
The pair of equations, (11) and (13) formally define a 
complete solution of the transport equation. 

For numerical calculations, it is more feasible to calcu- 


late the contributions to 7(R, 2) due to the various indi- 
vidual orders of scattering up to the order at which the 
contribution is negligible and then to take the sum of all 
orders. Obviously, (11) and (13) can also be interpreted 
for the individual orders of scattering. In general, the 
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required integrations cannot be carried out analytically, 
and the necessary number of reiterations may be very 
large if the dispersion is dense. 

A somewhat simplified form of integral solution may 
be obtained by introducing an approximate distribution 
of the kth-scattered radiation estimated from the dis- 


— 
tribution for single scattering, f;(Q, 2’). This approxima- 
tion, originated by Hartel [9], involves the assumption 
that the dispersion is so dense that only the radiation 
scattered from the immediate neighborhood of an ob- 
stacle; z.e., having the same distribution as from the 
obstacle itself, is of importance. The distribution of kth 
scattered radiation, with respect to direction of the in- 
cident radiation, Q, is then 

> => a a 
Fu(Qo, 2) = ih fir1(Qo, 2’) fr(Q, 2’) dQ’. (14) 
a 

Using these approximate equations, it is only neces- 
sary to evaluate the integrated power due to each order 
of scattering, and multiply it by the proper distribution 
function. Even with this approximation, the required 

calculations are too detailed for practical application. 


APPROXIMATE REPRESENTATION OF MULTIPLE SCATTER- 
ING BY THE CLASSICAL DIFFUSION EQUATIONS 


as 
If the specific intensity, 7(R, Q), is represented only by 
a two-term series of spherical harmonics; 7.e., by an 


ed 
isotropic function of R plus an anisotropic function of 


aad 
R, such as 

ee — => 

(R, Q) = io(R) + A(R) -2 (15) 
the transport equation can be degenerated to the classi- 
cal wave equation, or, in the case of negligible absorp- 
tion, to the Laplace equation. The solution of these 
latter equations, together with the appropriate bound- 
ary conditions, is relatively simple mathematically, 
compared to the general solution of the transport equa- 
tion and has been the subject of considerable study in 
classical physics. The over-simplified representation of 
the specific intensity by (22) yields useful results for an 
infinite dispersion or for several mean free paths from a 
source or the boundaries of a finite dispersion, but it is 
inaccurate near a source or boundary where the specific 
intensity is highly anisotropic. 

Specifically, the introduction of (15) permits reduc- 

tion of (1) to 


= — 
D; div grad ¥(R) + Now(R) = 0 (16) 
where the diffusivity 
D : 17 
for 3No(1 — wim) oh 
with 
a ns 
Ln = ii u(Q, 2’) f(Q, 2’)dQ’ (18) 
Q’ 
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The intensity is obtained from the expression 


(R) [grad y(R)]-2 
os R ra : 
Meat pee e Ba eA) 
Ag 4nrNo(1 — wim) 
For no absorption, ¢.=0, and (16) reduces to 
— 
div grad ¥(R) = 0. (20) 


Complications occur in applying the formal methods 
of solution of (16) and (20). The general solution is ob- 
_ ——s 


tained in terms of ¥(R) or 7(R, Q), but the boundary 
conditions in scattering problems are generally specified 
in terms of the power per unit area crossing a boundary 
in a given direction and cannot be satisfied by (19). For 


example, if no power enters a boundary, 7(R, 2) must be 
exactly zero for all angles over a hemisphere. 

Approximate boundary equations may be derived for 
certain specific geometries as a substitute for (19). For 
parallel-plane symmetry, the flux per unit area in the 
forward direction, Ez,, and the flux per unit area in the 
backward direction, Ez_, are approximately 


¥(Z) _ 1 dy(Z) 


SS a ee ZA 
4 6Noi(1 — wm) dZ (21) 


Ez. = 


if the third and higher derivations of ¥(Z) are neglected. 
It has been suggested that the inaccuracy of the dif- 
fusion approximation near a boundary can be mitigated 


— 
by introducing an imaginary boundary at which yW(k) 
vanishes. For example, from (32), it can be seen that if 

—> 


Ez -=0 at the boundary, ¥(R) would vanish at a dis- 
tance beyond the boundary 


Z 


= 3No(1 — wim) a 


Li 


called the linear extrapolation distance. Glasstone and 
Edlund [5] report that rigorous solution of the transport 
equation for parallel-plane symmetry yields a coefficient 
of 0.7104 rather than the 2/3 indicated in (22). 


DISCRETE FLUX REPRESENTATION OF 
MULTIPLE SCATTERING 


An approximate solution of problems in multiple 
scattering can be obtained if a simplified representation 
is used for the angular distribution of radiation from a 
single scattering. If the actual distribution for single 
scattering is utilized, the transport equation must be 
solved for all directions at all points in the dispersion. If 
the angular distribution is represented by one or more 
discrete components, or by an isotropic component plus 
one or more discrete components, the transport equation 
can be reformulated as a set of relatively simple differ- 
ential equations, and analytical solutions may be possi- 
ble. Such representations have been proposed by Theis- 
sing [11], Fritz [12], and Chu and Churchill [13]. 
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Thetssing’s Two-Flux Model 


Theissing formulated a model for spherical particles 
in which the scattered radiation is confined to two direc- 
tions, forward and backward. The continuous distribu- 
tion and Theissing’s model are sketched in Fig. 1(a) and 
(b). Hartel’s approximate distribution for the kth- 
scattered radiation is utilized, so that the forward and 
backward components of the distribution function for 
any order of scattering, k, are given by 


w/2 
te= anf fx(@) sin 6d0 
0 (23) 
by = 2 i fx(0) sin 000, (24) 
n/2 
i ; 
f (8) 
Re 
Scattering b 
Obstacle ! 


b 
(a) (b) (c) (d) 
Fig. 1—Angular distribution of radiation from a single scattering. 


(a) Idealized sketch of a continuous distribution; (b) Theissing’s 
model; (c) Fritz’s model; (d) six-flux model. 


Energy balances are then written for £;4, the forward 
power per unit area, and F,_, the backward power per 
unit area, for each individual order of scattering. If 
Theissing’s analysis is extended to include absorption, 
the equations for parallel-plane symmetry are: 


dE oy 
dz + NorEn, = 0 >) 
dE yy 
+ NorvEu, = Nosfikoy eo 
dZ 
dE\_ 
tts + NopE, = Nobis Foy Cy) 
or in general 
dE 
opt NarB sn = Newfie Burt Nobis Bs ue 
and 
GE (rp y— 
rf pt NesEusy-= Eo, rp Er-+NosbepiExy (29) 


where Z is again the distance through the dispersion. 
This set of equations is readily solvable by a stepwise 
procedure, but the number of orders of scattering which 
are appreciable, and hence the number of equations that 
must be solved may be very high for dense dispersions. 
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The accuracy of the method is limited by the simplicity 
of the model and by the accuracy of Hartel’s approxi- 
mation. The results should be better than those obtained 
by two-term quadrature, since the weighting factors are 
evaluated from physical considerations. 


Fritz’s Two-Stream Model 


Fritz represented the radiation from a single scatter- 
ing by a forward component and an isotropically distri- 
buted component as indicated in Fig. 1(c). The medium 
is divided into layers having a thickness of one-quarter 
of a mean free path for interception. The diffuse radia- 
tion (isotropically distributed) generated in each layer 
is obtained by numerical integration, as a parallel beam 
is traced through the medium, layer after layer. An 
analytical equation is then fitted to a plot of the diffuse 
radiation generated vs distance. Solutions are then 
sought for the classical diffusion equation with this em- 
pirical expression as a source function. The accuracy of 
the method should be somewhat better than diffusion 
theory for which it represents a first order correction. 


Six-Flux Model of Chu and Churchill 


The representation of the angular distribution for 
single scattering by six discrete components, one for- 
ward, one backward, and four sidewise, as shown in 
Fig. 1(d), permits reduction of the integro-differential 
transport equation into a set of six differential equa- 
tions in terms of the corresponding discrete components 
of the specific intensity. This set of equations can be 
solved for many geometries. The solutions are generally 
detailed, but can be evaluated readily on an automatic 
computer. The specific intensity can in turn be approxi- 
mated from the six discrete components. 

For randomly-oriented or for symmetrical obstacles, 
the four sidewise components of the angular distribution 
function are equal and the forward, backward, and side- 
wise components can be defined as 


i = fA cos? @ sin 6d0 (30) 
b =2r ; 1(0 2 @ sin 6d0 31 
le fi(@) cos? @ sin (31) 
and 
3=(1—f — 5/4. (32) 


The six discrete components of the specific intensity, 


——— 

1(R, Q) are chosen in the direction of and opposite to the 
three coordinate axes of any orthogonal system of co- 
ordinates as shown in Fig. 2. One coordinate is chosen 
in the direction of the incident wave. 

By introduction of Ez,, Ez, Er4, Er-, E,4, and £,_, 
the six components of the specific intensity, and f, 5, and 
5, the three components of the angular distribution 
function into the transport equation, (1), or by writing 
an energy balance around an elementary volume, a rela- 
tively simple set of differential equations are obtained. 
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Solutions were reported for a parallel-plane wave 
oblique to a parallel-plane region of scattering and for a 
point source outside a parallel-plane region. The solu- 
tions appear to be more accurate than any of the ap- 
proximate representations previously considered. 

For geometries having a high degree of symmetry 
such as a parallel-plane wave normally incident upon a 
parallel-plane region the six-flux model degenerates to a 
simple two-flux model. This two-flux model is suitable 
for hand calculations and the solution compares favora- 
bly in accuracy with the results of the diffusional repre- 
sentation. 
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Fig. 2—Representation of specific intensity by discrete components. 
(a) Actual distribution; (b) discrete representation. 


CONCLUSION 


Multiple scattering can be represented in an Eulerian 
system by an integral equation, called the transport 
equation, or in a Lagrangian system by a random-walk 
process. Approximate representation is possible by the 
classical diffusion equation or by a set of discrete com- 
ponents of the intensity. 

Exact solutions for the random-walk formulation can 
be found only in infinite regions and are difficult to ap- 
ply for anisotropic distributions. The Monte Carlo 
technique can be used to obtain particular solutions for 
any geometry and distribution of scattered power, but 
for a reasonably accuracy, the required calculations are 
very extensive, even for a modern automatic computer. 

Exact solutions of the transport equation have been 
found only for isotropic scattering in infinite regions. 
Approximate solution by quadrature, spherical har- 
monics, or reiteration can be carried out for parallel- 
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plane boundaries as well as for infinite regions but prac- 
tical applications are limited to isotropic scattering. 

Approximate representation by the diffusion equation 
permits analytical solution for many geometries and 
numerical solution for any geometry. This representa- 
tion is inaccurate near boundaries and for highly aniso- 
tropic distributions. 

The two-flux model proposed by Chu and Churchill is 
more accurate than diffusion theory but is practical only 
for simple geometries. The discrete flux models of Theiss- 
ing and Fritz are perhaps more accurate but require very 
extensive calculations. The six-flux model of Chu and 
Churchill is considerably more accurate and is adapta- 
ble to many geometries. However, the analytical solu- 
tion obtained by the six-flux method is so detailed that 
calculations are feasible only on an automatic computer. 
All of the discrete flux methods can be used for aniso- 
tropic scattering. 
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Diffraction of Plane Electromagnetic Waves by a 
Rectangular Aperture* 
MICHIO SUZUKI} 


INTRODUCTION 


HE STEADY state problem of the diffraction by 

an aperture in an infinite plane screen has at- 

tracted attention for many years. Available solu- 
tions have been restricted to a few cases where the 
aperture is circular, elliptical or infinitely long and, for 
rectangular aperture, no adequate theoretical treatment 
has been given as far as the author knows. This paper is 
concerned with the theoretical study on the diffraction 
of plane electromagnetic waves by a rectangular aper- 
ture on an infinite plane screen having perfect conduc- 
tivity. In order to calculate the transmission coefficient 
the integral-equation method, which can be solved by 
method of successive approximations, is used when the 
aperture is long and narrow, on the other hand the vari- 
ational principle is used when this is not the case. 


INTEGRAL EQUATION FORMULATION FOR A RECTANGU® 
LAR APERTURE, AND THE TRANSMISSION COEFFICIENT 


We consider an infinite plane conducting screen of 
infinitesimal thickness, which is perforated by a rec- 
tangular aperture 7. A rectangular co-ordinate system is 
chosen with the origin at the central point of the aper- 
ture, and¥oriented so that the screen lies in x-y plane 
(Figo. 1). 


Fig. 1 


A plane electromagnetic wave is incident upon the 
aperture in the half space Z<0, and it is desired 
to investigate the diffracted field. The plane wave 
E, = Eye" (k is free space propagation constant) is in- 
cident from the left, and sets up a electric field in the 
aperture. The incident field is polarized in x-direction 
only, so that the field in the aperture is also polarized in 
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x-direction only and its amplitude distribution is e(x, y) 
(as yet of undetermined form). It is our first object to 
write the fields in the both regions expressed by Z S0 
and Z=0, and then to equate them on the aperture. 

For this purpose we introduce a Hertzian vector M* 
of the magnetic type, and express the electromagnetic 
field (E, H) as follows: 


E = — jowoV X II*, H=VxXV xX II* (1) 


For our problem II* is taken to have only a single 
component II,*. Hence the field components are 


eolige alive 
z = Jopo BD oy 
Oz Oxdy 
0711 ,* 
E, = 0 ; Hy, = a eae 
oy? 
nee Ollie oll ,* 
E, = — jopo H,= : (2) 
Ox 0ydz 


Subscripts 1 and 2 will be used to distinguish the quanti- 
ties in the regions ZS0 and Z=0 respectively. 
For the region J(Z S0), 


Ls = IIo* | Ges (3) 
where 
Ey , 
ines = (Cae + eike),, (4) 
WO 
1 Ge ike 
yt =-—— | én) atin, 9) 
2r7JwWpo J ; r 
and for the region 2(Z 20), 
1, =e (6) 


£, 7 are the coordinates on the aperture, and r is the dis- 
tance between a point on the aperture and a field point 
x, y, 2; 


CR ess Cee eee 
It will be confirmed that 


_ OIIo* 
lim = 
20 OZ 
. ( aa ioe r 
lim { — = 
20— Oz 0 0, 
oln* x, T 
tim ( “) = " y) (7) 
z-0+ Oz 0 o 


and, that II,,* and II.,* satisfy the wave equation in the 
respective regions. 
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The boundary conditions for the electric field at Z =0 
are that they should be zero on the metal surface, and 
continuous and equal to e(x, y) on the aperture. From 
(2) and (7) these conditions are satisfied at all. The re- 
maining boundary condition to be satisfied is the conti- 
nuity of H, over the aperture. When this has been satis- 
fied, it is easy to see that the boundary conditions for 
the remaining components are necessarily satisfied. 

Equating H, on each side of the aperture, we find 


2 * 


TERA ay 0) = oh lly. = Yop on 7, (8) 


in which 
yo = V€0/po- 


Eq. (8) with (5) is an integrodifferential equation for 
the aperture field distribution e(x, y). 

Now let us introduce the transmission coefficient of 
the aperture. The transmission coefficient T of the aper- 
ture is defined as the ratio of the transmitted energy 
flux W through the aperture to the incident energy flux 
w upon the aperture, and expressed as follows: 
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Il,*(«, y) = (13) 


=i SOE + C(a)e7t#u + D(x)ei*v| 
where C(x) and D(x) are some functions of x. 

If the width a of the slot is small compared with wave- 
length and the primary wave is incident perpendicular 
to the aperture, it will be confirmed that we can put 
C(x) =C(0) (=6), D(X) =DO): (=D) and C=D=B/2 
(B is the unknown constant). Equating the right-hand 
sides of (13) and (5), we find 


{«@ ) se 


This is the integral equation for e(x, y). After finding of 
e(x, y), T will be determined by (12). Referring to Gray’s 
paper [1] let us rewrite the left-hand side of (14) as 
follows: 


f e(E, 


fi {e(E, n) — €(&, y)} cos kr — je(&, ») sin kr 


r 


. 


(14) 


Lediyaet EO pany = [Ho + Bos ky]. 


didn 


W and r = \/(x—£)?+ (y—7)? in the second integral may be 
yi ae (9) written as y—n| because } is much greater than a. 
Using these results, (14) becomes 
a a 7/2 te(E,n) — €(&, y)} cos k| y — | —Je(E, 9) sin k| v — 
f e(&, y) K(x, &; y)d& + fe | | | ae 
—a/2 —a/2¥ —b/2 | yes | 
207 
= a (Eo — B cos ky), (15) 
in which where 
1 a b/2 k 
W= zs Re f a, y)H2y(x, y, O)dady, (10) K(x, &; y) = i mes an 
4 T b/2 r 
1 
w= ai abyoEy?, (11) 


where A», is the conjugate complex of Hz, and a and } 
are the sizes of the aperture as shown in Fig. 1. 

Let us eliminate H», in (10) by means of (8), then (9) 
becomes 


T = (12) 


Re f , y)dxdy. 
aE. e eA: y)dxdy 


THE CASE WHEN THE APERTURE Is 
LONG AND NARROW 


At first we consider the case when the slot is long in 
y-direction and narrow in x-direction. The solution of 
the partial differential (8) will be the form 


(243) + 0-04 (245) 
/(- y) + -8- (+-») 


~cufi(}} ~ c(h 


Cin(%) 


log 


= log x + 0.5772 — C,(x), 


C,(x); integral cosine. 


Since the function K(x, &; y) is slowly varying func- 
tion of y, it may be replaced by its mean value 
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1 b/2 

K(x, ) = — K(«, & y)dy 
b J_o;s 

b sin 28 
=2 (1 = Cx(28) — 
[og — =f lors = Cy(2a) — |,ao 
where 
es wb 
max 


On changing the variables x and é to @ and 6’, respec- 
tively, by putting x =a/2 cos 6 and £=a/2 cos 6’, where 
6 and 6’ vary from 0 to 7, and using the formula 


2 cos mé cos mé’ 1 
pa Se EE See 7 108 2| cos8 — cos 6’ | , 


1 m 
we get 
2 cos mé cos mé’ 
KGS) = M% +4 >> —_——— (17a) 
1 m 
where 


4b sin 26 
Qg = 2| tog — Cin (28) — —— + log 2|. (7b) 
a 268 
Expanding s in e(&, y) in a Fourier series as follows: 


sin 6’e(, y) = >> Cn(y) cos mé’ (18) 
0 

and substituting this into (16) and making some simpli- 

fications, we get 


© Cm 9 
>, | acolo) eo os | 
2 1 mM 


20j Ta 
iF Gos each eee: by 


—b/2 
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where 
Ly) = Sisk d S54 k P 
= 4 (= +»)} sie 4 (=-»)t, 
b 
M(y) = sin ky E {28 a + y)k 
hee ae ( ) 
PHS = 9) | 
— ky| Cin<2k p 
a fen fu 9 
Cin s2k : 
+ Cu 424(5—2)t | 
N(y) = ky| S;42k : t 
N(y) = cos y| Si (= +9 
+ S;42k : 
ees —7)} | 
+sin by €. foe(— +3) 
wn 2 
b 
-exfa(S—}] 
S;(x);  integral—sine. 
12 {Co(n) — Co(y)} cos &| y — 7| — jCo(n) sin | y — | Hees 


ly —a| 


Equating the coefficients of cos m#(m=0, 7, - ++) in Now the boundary condition for e(x, y) at y= +b/2 is 

both sides of (19), we get the following relations: ee; y= +b/2)'=07/thattis, Co(-Eb/2)=0: 

COs CY) “esp; 

f 1 12 }Co(n) — Co(y)} cos k| y — n| — 7Co(n) sin k| y — 

Coly) = <- e+ Beos by) - — f iC of) UI) aa Ue (ey 
OQ a ¥ —/2 | Baars: n| 
where From this condition the coefficient B can be deter- 
_4 mined. Hence the electric field distribution e(x, y) in the 
(Sate a aperture becomes as follows: 


Eq. (20) is also the integral equation for Co(y), and can 
be solved by means of the method of successive approxi- 
mations which HAllen used in his paper concerning the 
theory of electric current antenna. Thus we get the 
following solution of Co(y). 


f 1 
Coy) = =| Eo Tee COSY “Pia {ioL() 
Qe Q¢ 


B _B 
- Su +15 xo |, (21) 


ous Canes ath Gol): aig a Co(y) 
a4 sin@  4/I — cos? ie ey 
Vt es 
4 1 
= vee Soe E — dcos ky 
2 


ce 


+—fi) + Mo) -5-Noh |, 22 
Qe 2 2 
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where 


HU 
14+ — S,(2 
aay (28) 
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(23) 


1 
cos B — ae [sin B5;(48) — cos BCin(48) — j {cos BS:(48) + sin BC in(4) } ] 


G 


The aperture electric field distribution in x-direction 
is given by [1—(2a/a)e]|-"/?, and this distribution seems 
to be correct, since the width of the aperture is so small 
compared with wavelength that the field distribution in 
the same direction may be considered just the same as 
static field. 

Substituting the e(x, y) from (22) in (12) and carrying 
out the integration, the transmission coefficient JT is 
given as follows: 


Ha 1 Im| 1 a sin B 
adg B 
1 i ae ‘a 
a ee EM phi (a+ jad |, (24) 
where 
IP cos*26 
$9 = 25,(28) = ===>) 
B 
gi = 2{sin 26Cin(48) 
— cos 265(48) + S:(68) — S,(28)}, 
qo = 2{sin 26S,(48) 


+ cos (28)Cin(48) + Cin(2B) — Cin(6B)}, 


and Q¢ and 6 are given in (17a) and (23) respectively. 
Now if we rewrite (14) as follows: 


e(E, y ihr — ¢ 
f ED) pn mp e(E, ne (, ») jw 
T r e T 


r 


21] 
= — oy (Eo of B cos ky), 


another solution for 7 is obtained by the similar way. 
This is as follows: 


sin B 


Tv 
r= ——Im| 1 - 


aQy 


1 
a5 ao E28) Boe ea 


ye 6 


YC ey ay 
j26 


F = et+®E(4B) — j2e-#®S;(28), 


where 


E(x) = Cols) 95%), 
4b 
Qg = 2 log — 
a 
il 
D-- ——- Pi2B) 
Qx 
ss (25) 


These formulas are available for the case when the 
aperture is long and narrow. 


‘ 


THE SOLUTIONS BY THE VARIATIONAL METHOD 


In this section we will find the solution by using the 
variational method. Using (5), (8) and (12), the sta- 
tionary expression for the transmission coefficient T of 
the aperture, as derived by Levine and Schwinger [2], is 


(27)? (=) 

as Tras =) 
dab F 

G= fe y)dxdy, 


Ff f domme en (5 s+ He) 


where the notation Im means the imaginary part of. 
F is the quadruple integral and this is very troublesome 
to evaluate, but it can be treated as follows. Now let us 
introduce the integral representation for e~*"/r, that is 


e ike 


dtdndxdy, (26) 


en akN (a—-8)°+ (yn)? 


Va-+ Oo =e 
e-taly—n)—V at BP #| 28] 
== ll ic A/a 
2 * cos a(y — nerve tele _ El 
f atl: na: 
As y and 7 vary from —b/2 to 6/2, z=y—n varies from 


—b to b. Expanding cos aZ in a Fourier series in the 
range from —b to 0d as follows: 


dadB 


dadg. (27) 


cos aZ = >> A, cos q,Z, 


(28) 
n=0 
where 
1 b sin ab 
Aj —| cos aZdZ = ) 
2b 0 a 
1 b 2 sin ab (—1)” 
An=— f cos aZ cos gxZdZ = 2 Sa os ae 
bd b a? — gy? 
nr 
Qn b , 


and poe this into (27), we get 


Cot 2 co 


2 i Styreseotan ean it T,(B)d8, (29) 


r ee bee 


1(n # 0) 


a Sin aBe~Ve+8°—k*|2-€| 


igre ih de, 
nos ante 


Integrating the following functions f,(y) and f_(u) 


(30) 


pet tHe Vu +8? k| a8 | 


Wt — viet 


i= 


along the contours C; and C, respectively as shown in 
Fig. 2, it will be shown that [3] 


e- VB en"| z—€ 
18) = (1) > ee (31) 
where 
Ope k* 1a. 


Inserting (31) into (29), we get 


jkr © 
« : > wXn(| x — £|) cosqa(y — ), (32) 
where 
2% g—VB—an*|2—t| (33) 
( x — |) Me pees 


The integral X, (|a— é|) can be evaluated as follows: 
1) If a,?<0, putting an?= —atn! (An! >0) and chang- 
ing the variable from 6 to ¢ by putting V/6?+-a,/ = ay! t, 
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we get 


| 2} ean’ |2—§|s 
x — £/) -{ 2 
Te Ae eh) 


by using the well-known formula, where Ko(x) denotes 
the modified Bessel function of the order zero. 

2) If a,?>0, the integral (33) can be divided into two 
integrals as follows: 


| | an ei V an®—6?| 2 | 
XG) = =-—7] ft ee d 
a — £|) a es ay Serr: B 


Xnf 


Kolan’ | oo E| ) (34) 


eV Ban? | — an? | x—£| 
8), 
ale: Ve nia 


Changing the variable from 6 in the first integral 
to @ by putting Va,?—6?=a, cos ¢, the first integral 


becomes 
w/2 
— if e 7en| 2—§| cos edd. 
0 


Employing the Jacobi’s expansion in the series of 
Bessel coefficients, 


(35) 


eiZ os = JZ) +2 D (—1)™Jom(Z) cos $s Imdb 


m=1 


eg PA ss (—1)"Jom4i1(Z) cos (2m + 1)¢ 


m=0 


(36) 


and inserting this expansion into (35) and integrating 
term by term, we get 


a —&|) — 2j Danes 


m=0 


vg amt (On | Loe 


PMO ee Nh 


‘De 


ia = Tole 


«—£|) becomes as follows: 


§ . ay ‘ 
AG! x — £| ) = — jy Tole | De E| ) 
= Dom eae E 
ir = yee ee 
m=0 \ 2m + il 
+ Ko(on| x — &|). (37a) 
If w=na/dX is very small compared to unity, 


Dae x —&€|) may be approximated as follows: 


1 J 
X.(|# =£|) = b| 341 cage oe pt 


8 | | tan ei 
x — og ————_— 
; | (E ae Lan | 
where 6,=1, if a», is real, and 6,=0 if Q, is imaginary. 
Now let us assume that e(x, y) will be a product of the 
functions of x and y separately, that is, 


, (37b) 


— &y 


e(x, y) = e1(«)eo(y). (38) 


Using (32) and (38), (26) becomes, after some simpli- 
fications, as follows: 


154 IRE TRA 
T 1 
etait (=), (39) 
2a D 
where 
20 Qn 
D= De (=) ped (40) 
a/2 
fe ex(x)er() Xn(| i | \dadé 
r, = rik ; (41) 


Hieaoed 


b/2 
i i sles icos eave sanane 
es all mae co) 


6/2 2 
| {h aly)dy | 
—b/2 


1) At first let us consider the case when aX. 
In this case we will take as a trial distribution of 
€i(x) the following function: 


wf OT" 


Inserting (43) with (37a) into (40) and carrying out 
the integrations, we get 


(43) 


ol pate ucolpelea) eel 
ie 2-24 (44) 
| an | a/2 


2) If we put e(x) =J, (40) is evaluated as follows: 


Let F(A) bea function of Xand fo\F(u)du =F, (A), then 
we find, by integration by parts, 
Dre 
f f F(u)dudd 
0 0 
D 
= { Fea = DRO)? -f NF(N)AN 
0 


= fo — \)F(A)A. 


Similarly, 


eo = fe — A)F(A)a(a)e. 


Using these results to carry out the x- and £-integra- 
tions in (41), we get 


2 “9 Cs 
1D. — |. a, j— ais: — A)Jo(and)dr 
a? 2 


sop 


mao 2m +1 zt (4 — AJ ants(aud)ant 


+f (a — N)Kalaad)ar |. 
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Using the following formulas: 


ae (x°J (x) | = «J ,1(X), 
Gigs 


— x°K,-1(x), 


I 


d 
AE [xeK p(x) | 


I 


f sade = 2D Iponl®) 
v=0 


and employing the power series solutions for Jp(«) and 
the ascending series solution for Ko(«), we get 


2 

Ds = E {- 4] — = (25 DS J o41(And) — J.(as0)) 

And 2 y=0 

ee) 1 °° 
m7 BO aE (2 bed omp2v4-2( And) + Pants(oat)) 
a,aKy(ana) — 1 
(ana)? 

where 


nd 2m+2v+1 
os ae 
am+1(Qn@) = (and) a7 y!(2m - +y+1)!(2m + 2v + ai 


oni 
CRO Ne D 
2 De (m!)?(2m + 1) 
=f)" =n" 
“ a Pa 2 


mao {ml(2m + 1)}2 seco (m!)2%(2m + 1) 


Qland) = — log ( 


1 1 
Y(m) = — 0.5772 +1 +—+4+---4+—. 
2 m 


3) If we put 


i 4 2y 21/2 
a) =| Gi 
(42) becomes 


= fon fenton 


-sin? 6 sin? 6’d6d0’ 


(46) 
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after changing the variables from y and 7 to 6 and @’ re- 


And 2m 
spectively by putting y=(b/2) cos @ and 7=(b/2) cos 6’ SN (<*) 
and using the Jacobi’s expansion in series of Bessel co- Q(a,a) = — log (=*) > 
efficients. 2 7 mao (m!)*(2m + 1) 
4) If we put é(y) =cos B—cos ky, then (42) becomes Ana?” ana\*™ 
as follows: & Te oo v(m) 2 


+> 


ee | peo AES CP MER NIT MY Soles le <I) 
1 1 
(6/2)? UG) = 05772 1 ee 
Ag,.—k = , 2 m 
(8 cos B — sin B)? 


2 2 : 
—BcosB 4 
T 


Ame te ey (48c) 
Amp = |—————_| | ae. 
B cos 8B + sin 6 2 
1 Ao SS - 
jet ie ces aL? Teen (8/2)? 
Pea Gra Sens une 
WY 2 
7 2 1 — 
ime ( sin B ) ethers = B cos B 
B cos B — sin B Qnum\?) 2 Asati = | 
{1-(=)} B cos 8 — sin B 
28 E 
1 2 
The results so far obtained may be summarized as : 2h + ila\\? 
follows: (2n + 1) 41 — 26 
1 ; 2 1 
ea pil Im (=) Aon = <5) Maye ste ee nine Rate, (48d) 
2a D B cos B — sin B 


ebay 
where 26 


nt 5) Finally, if we take e(y)=cos gy+A cos gsy and 
2 Jp determine the constant A by putting 07/0A =0, then 
D=—AIo+ >, \1 i 7 T,An, (39) becomes as follows, after some simplifications, 
1 
rg ty + t33 — =| 
“| ie = — im 5) (49) 
ae 2a | tiitss3 — tis” 
Bee oe (ay (2) (2)] wie 
ied. 2 Aa? tT 2 ; 
(em 
2:24 1 eS f y 
lo (48b) ihe ar tote To — a 1 es he EAC arene, 
ave | Qn | a/2 2 eat B 
2 2 e with 
1 en = E \- j oy (2 ys J oy4.1(n@) - 1:(a,0)) ar 
Gnd baa Ausi1 ) Aiz1,3 = Asz;1,3 = 9, 
16 
2) il ve) 
—2 2 Jf m-+2y- (Qn@) ms Pom (axa) ue 
eal 2 Jamsints ite A3;3,3 = 9(—], 
16 
Oma: ana K (ana) — -| 1 x 9 
And 5) Az; Fi leah Saag GR eps Pee 
(ana)” es (n? — 1)? ty (n? — 1)(n? — 3) 


2m+2v+1 3 2 
ont Anj3,3 Ee 0( H 
A) 9) n ae 32 
Pon4i(Qnd) = (a,@) le 


“(2m ty + 1)!Qm + 27+3) and Ty is given in (48a) or (48b). 
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Successive approximation due to Gray 
x Successive approximation due to HAllen 
I] — - — - — Variational approximation, based on aperture 
field of the form 


Qxn\ 29-12 ay 2741/2 
on eo) eel ea lll ibe ela 
Variational approximation, based on aperture 
field of the form 


2-1/2 
e(x, y) = [1 — (=) | [cos 8 — cos ky | 


a 


fe) Variational approximation, based on aperture 
field of the form 


Qn\29-12 3 
dey) =[1- (=) ] [ cos 9 + Coos y|. 


NUMERICAL EXAMPLES 


In this section two numerical examples are shown in 
Fig. 3 and Fig. 4. The transmission coefficient is plotted 
in Fig. 3 as a function of 8 in the case a=0.1. The 
curve I is calculated from (24) which was obtained by 
means of the integral equation method. 

The curve II and III are calculated from (48a) with 
(48c), and (48) with (48d) respectively. The points indi- 
cated by crosses are obtained by (25), and the points 
indicated by circles are obtained by (49) with (26). 
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fe} 0.5 1.0 5 2.0 2.5 3: Ome 


Fig. 4—8=1.5. Variational approximation based on aperture field 
of the form e(x, y) =(cos 8—cos ky). 


When the length 0 of the slot is approximately equal 
to \/2, 3/2, 5/2, - - - , the resonance phenomena oc- 
cur. At these points the transmission coefficient takes 
the largest values. 

The transmission coefficient is plotted in Fig. 4 as a 
function of a in the case @=1.5. This curve which is 
obtained from (48b) with (48d) approaches to the value 
unity as it is increased to infinity. 
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The Amplitude Concept of an Electromagnetic Wave and 
Its Application to Junction Problems in 
Waveguides” 

JEAN A. ORTUSIt 


T IS generally known that the concept of complex 

| amplitude of electric and magnetic fields is widely 
used in mathematical physics. 

In this paper, it will be shown that this concept can 

be extended to guided waves and that it enables various 


* Manuscript received by the PGAP, June 10, 1955; revised 
manuscript received December 12, 1955. 
t+ Compagnie Generale de T.S.F., Paris, France. 


functional properties of junctions used in waveguides to 
be computed and represented. 


CALCULATION OF A GUIDED WAVE 


Denoting by U(x, y), the complex amplitude of the 
longitudinal component of the hertzian vector, it can 
be readily shown that the longitudinal and transverse 
electric and magnetic fields are proportional to U(x, y) 
or Whee. 
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A guided wave will be defined as normalized when the 
total flux of the Poynting vector in direction Oz, inte- 
grated during a period of time, is equal to unit power. 
The concept of a normalized guided wave allows for the 
unambiguous definition of the complex amplitude of a 
guided wave. A nonnormalized wave of a given type is 
defined by the value of the function U~ 


U= "GV (x, 3). (1) 


C is a constant and V(x, y) is obtained by solving the 
propagation equation and taking into account the 
boundary conditions: 


V=0 


os 
e;;ds = 9 at the contour. 


The function V(x, y) is characteristic of the type of a 
guided wave. 


COMPLEX AMPLITUDE OF A GUIDED WAVE 


Let the real value of the constant C for the normalized 
guided wave, characterized by the function V(x, y), be 
designated by C\. By definition, the term complex am- 
plitude of the guided wave will be applied to the ratio: 


G 


iL es 
C1 


(2) 


The power transmitted by the guided wave is therefore: 
E = D-D* (3) 


where D* represents the conjugate value of D. 


TRANSMISSION TENSOR 


Let us now suppose that there are m guided waves 
each characterized by functions V;(x, y) which are all 
different. When part of the power £; circulating in the 
corresponding waveguide is distributed among the n—1 
others, the arrangement of 7 guides is called a junction 
or branching circuit. Let A; designate each component 


of the incident vector A with m components, represent- 

ing the complex amplitudes of the 2 incident waves in 

an n-dimensional complex space. Similarly, B; will be 
— 


the component of the transmitted vector B represent- 
ing the whole of the complex amplitudes of the trans- 
mitted waves. Because of the linearity with respect to 
the values of D of the solutions, we can write: 


= 


many. (4) 


where J represents the tensor having for matrix the 
square table lay out 


SPP oe he OBE Cece By 
TET ueweny Ne eae ae (5) 
PB ns hoor WW bey Oo IOESL D Bene 


Tensor T will be called the transmission tensor. 
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Each of its terms T;; represents the coefficient of 
transmission from guide z into guide j7. The term 7%; 
represents the reflection coefficient R; in the direction of 
the incident wave (R;=T;:). The term coupling factor 
fi; (or transmission equivalent between guides 7 and 7) is 
applied to the ratio: 


(6) 


T;;* designating the conjugate quantity of T;;. E; is the 
power incident into guide z, and £,; the power trans- 
mitted into guide j when all other components A; (k 47) 
are zero. 


RECIPROCITY THEOREM 


Let us suppose that the media contained within each 
guide are all isotropic (€; and uw; are then scalar values 
differing in each medium). It can be shown! that the 
transmission tensor is then a symmetrical tensor. 


(7) 


This property is independent of the choice of the 
planes of origin. Formula (7) involves two relations: the 
transmission equivalent (or coupling factor) is the same 
for the two opposite directions of energy transfer, fi; =fji, 
and the transit time is also the same in both directions. 

Eq. (7) constitutes the reciprocity theorem. It is im- 
portant to note that this theorem is invalid when at least 
one anisotropic medium exists (ferrite or electronic 
plasma). 


Tay = Ty. 


STUDY OF THE TRANSMISSION TENSOR 
IN LossLess MEDIA 


Let us suppose that the media under consideration 
(electronic plasmas, isotropic, or anisotropic material 
media) are without electromagnetic losses of any kind. 
Also, since the media may be anisotropic, the transmis- 
sion tensor is not necessarily symmetrical and the rela- 
tion 7;;=7;; between the tensor components will, in 
general, not be observed. However, because of the ab- 
sence of electromagnetic losses, a certain number of rela- 
tions must exist between the transmission vector com- 
ponents. This new form for the tensor will be sought by 
supposing, therefore, that all media carrying waves 
introduce no electromagnetic losses. 

This case results when the symmetrical components 
of the electric and magnetic induction tensor are real 
and when their antisymmetrical components are purely 
imaginary. In isotropic media, € and yp are real numbers. 


1 This demonstration is given in: C.R.A.S. Part 217, p. 677, 1943. 
Note by Messrs. Gutton and Ortusi, “On the form of the reciprocity 
theorem on hertzian waves.” It is important to note that (7) applies 
to amplitudes D of (2) and not to the amplitude of the electric vector 
or even the amplitude C of (1) which can be different for the two 
senses of propagation if the waves of guides 7 and 7 are not similar. 
This point is an important subject of error in the application of the 
reciprocity theorem, in particular when the media are not the same 
even with similar waves. 
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Determination of the Necessary and Sufficient Relations 
Between the Transmission Tensor Components in Loss- 
less Media 


Let us designate by L, the total incident power 


>> 
Le = >, Ay: A* = A-A* 


and by Ly the total power transmitted 


> 
Le = >) BB; = B-B*, 
— — —s — 
where A* and B* are the conjugate vectors of A and B. 


—_-? 
L, and Ly then represent the lengths of the vectors A 


ay 
and B. The theorem of the conservation of energy, in 
the case of lossless media, shows that L,=Ly. The 
transmission tensor is therefore an orthogonal tensor (pre- 
serving angles and lengths). 

This condition is necessary and sufficient. A junction 
of m guides, without electromagnetic loss in the media, 
is therefore characterized by an orthogonal transmission 
tensor. The components of this tensor are thus con- 
nected by the classical n? relations which fix the changes 
of coordinates on rectangular axes. 


Let T* and T, be the conjugate and transposed ten- 
sors of the tensor T defined by their respective com- 
ponents 7;;* and Tj; (tensor T; is therefore identical to 


tensor T when the media are isotropic, but differs from 
it in the contrary case). 


Writing that tensor Ts orthogonal, it will be seen 
that we have: 
Tals esov Tee a) 
T;—! is the inverse tensor of the transposed tensor of iis 
The components of Tas are given by Krammer’s rule. 
We have: 


ie aes 


where M,; represents the minor of component 7;; of 
tensor T and A the determinant of (5) formed by the 


components of ie 
We thus obtain the 2? equations which give the com- 


ponents of the conjugate tensor T*. 


Mi; 
yt ha 
A 
which can be written: 
Mi (1)? AG 
a Det i 10% 8) 


A;; being the determinant obtained by suppressing 
column 7 and row 7 in (5). 

It should be noted that (8) is valid for any type of 
wave and for any media, isotropic or not. Eq. (8) thus 
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constitutes the necessary and sufficient conditions to 
produce lossless junctions. In particular, it comprises 
the following relations (obtained, for instance, by de- 
veloping A along each row and column): 


1 whatever the value of 7. 


j=n 

* 
Dean 
j=l 


j=n 
>> Tis: T* = 1 whatever the value of i. 
j=1 
A has a modulus equal to 1. We can therefore write 
A=e?; ¢ will be called the junction phase angle. 
Eq. (8) is usually written in the equivalent form 


= e'7@ whatever thevalues ofiandj. (9) 


(ag) 


+A, 
thee 
It will now be shown how (9) can very simply define 
the principal categories of junctions and the various 
functional properties of these junctions of guides. 
Junctions oF Two GUIDES 


Isotropic media is given in (5), in this case, by making 


Shee al 
T 12 a T4 = fe T Ry ° 
Eq. (9) becomes: 
Rs Ri I 
RR, — T? = — = — = — — = @* 
Ry* Rot T* 


R, and R, having the same modulus, the planes of origin 
may be chosen to give the relations: 


Ry = Rp = cos ue. 
With respect to those planes, we have: 
T =isin u-e*, 
u is the reflection angle (—1/2<u<7/2): 


0<u<7/2 corresponds to a capacitive obstacle; 
—7/2<u<0 corresponds to an inductive obstacle. 


It is easy to see, when the longitudinal dimension of 
the obstacle is small, that u is related to the normalized 
susceptance Z by: 


Zz in 
= — — tgu 
7 


The measurement of u (and consequently of Z) is 
made, in the laboratory of the Compagnie Generale de 
TSF, according to the method of Fig. 1. 

The obstacle is put at an approximately equal dis- 
tance of \/2 between a detector and a conducting wall. 
The frequency of the source is modified in order to have 
two zeros at the detector which gives the u angle by 
the two corresponding values of the electric distance 6, 
6,=0 and #.=u. This is a very accurate method for 
measuring small values of the u angle of obstacles used 
for uhf filters of large surtension. 


1956 


Nonisotropic Media-Gyrator 


The tensor T is no longer symmetrical when the media 
are not all isotropic. By definition, when Ty,= — T.; the 
two-way junction forms a gyrator. Eq. (5) becomes a 
gyrator with lossless media, by making 

Kop Tr 
we Ks 


Ty = —Tn= 7 


Eq. (9) becomes: 
RAVER ; 
RyR, + T? = — = — = — = ot, 
ee Rs thy 
It is thus seen that, for planes of origin so chosen that 
Ry = Ro 
R = cos u-e# 
T = sin u-e*, 
A tuned gyrator (Ri=R2,=0) corresponds to the case 
where u =7/2; we then have, making ¢ =0, the relations: 
Li ll and 1b ae 
The incident wave, in one direction, suffers no change 
of phase and no reflection. In the other direction, it suf- 
fers a phase change of 180°, without reflection; this con- 


dition is equivalent to a transit time, modified by half 
a period. 


DETECTOR 


Fig. 1 


JuNcTION AMONG FouR GUIDES 


As an example of application of (9), the necessary and 
sufficient conditions will be sought which will give, in 
- lossless media assumed to be isotropic, a matched di- 
rectional coupling. The latter is characterized by the 
matching requirements: R; = R2= R; = R,=0 and that of 
' decoupling 713;=0 (guides 1 and 3 being the decoupled 


guides). The matrix of tensor T is: 


0 Ti. 0 Tus 
Pix Over Poq OL o4 
0 T 23 0 T 34 
Die Les Fae 0 
Eq. (9) then becomes: 

A =-Al/2 nee Alle. Ti mea? Tos 

2 T33* ites 
E 0 
SSO vi 2 eit aa = ¢tid (10) 
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Designated by fi;=7;;:T;;* and by $.; the coupling 
factors and the arguments of the transmission coeffi- 
cients 7,;, (10) becomes: 


fis = fos = cos? w 
fis — S23 = sin? w 


(11) 


2 = 0 
giz + bas = x2 
o14 + G23 = — t+ 26. 


Eq. (11), which is quite general, is valid whatever the 
form of the matched hybrid junction. In particular, it 
allows the following important theorem to be stated: 
matching of a hybrid junction, in which guides 1 and 3 
are decoupled, necessarily involves complete decoupling 
between the other two guides 2 and 4. 

The hybrid junction is necessarily formed by a series 
T and a shunt J of which the coupling factors are 
complementary. In the directional couplers, made with 
two distinct coupled guides, the phase angles of the di- 
rectly transmitted wave and of the wave transmitted by 
coupling have a difference of 7/2. 


THE FUNCTIONAL PROPERTIES OF THE JUNCTIONS WITH 
ANISTROPIC MEDIA 


Anisotropic media provide certain functional proper- 
ties in junctions which, from fundamental principles, are 
impossible of realization in isotropic media. In this sec- 
tion, two special cases will be examined. 


Circulators 


By definition a junction will be denominated a circu- 
lator when it consists of 2 guides defined by the follow- 
ing conditions regarding the components of its trans- 
mission tensor: 


T;; = 0 except for 7 =i orj=i+1 


with, of course, T7140 (n+1=1). 
A circulator is shown in Fig. 2. For example, a 4-guide 
circulator will have a transmission tensor represented by 


the following: 
Rie 0 Onties 


Ti2 Rp O 0 
O T2, R3; O 
0 Ousitisat ia 


It will be noted that the circulator must necessarily con- 
tain at least one anisotropic medium, since the trans- 
mission tensor is not symmetrical. 

Matching Conditions for Circulators: It can be shown 
that, in the case of lossless media, a circulator must 
necessarily be matched. 

Eq. (9) then becomes: 


7#i+1 and j¥i 


Ilo; RJ] 0; 


ee et a2) 
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making 7;, :41=6; where the symbol II represents the 
product of the terms excluding terms whose numbers 
are specially indicated. These equations show that all 
R; coefficients are zero. Some mistuning introduced in 
one guide of the junction necessarily introduces an un- 
desired coupling with one of the guides numbered 
j#i+1. When the circulator is matched, all 0; coeffi- 
cients have a modulus equal to 1 and their phases satis- 
fy the relation 


od; = nd. 


These conclusions do not apply in the case of a lossy 
medium. 


Use of Circulators 


Associated with band-pass filters, circulators are a 
very elegant means for effecting the junction of radio 
links. The diagram of Fig. 3 shows an arrangement with 
a single circulator with m guides. 

Guide 1 contains the initial wave frequency modu- 
lated from the value F=F, up to the value F=F,_4. 
Each guide 7 (with 2$7<n—1) terminates in an anten- 
na or in a receiver tuned to frequency F;. It contains a 
band-pass filter which passes, with a small reflection co- 
efficient, the wave of frequency F;+AF; and reflects 
with a reflection coefficient close to 1 waves of frequen- 
cies F; (7 #1). 

Guide m terminates in an attenuation impedance and 
does not contain a filter. The mode of operation of the 
arrangement is immediately deduced from an examina- 
tion of Fig. 3. No reflected wave can return to the input 
guide 1 whatever the quality of the filters F;. Inade- 
quate attenuation outside the passband produces linear 
cross talk in the transmitted wave (or the received wave, 
as the case may be). Too great a reflection within the 
passband produces an absorption of energy in the 
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terminal impedance of guide x. No nonlinear cross talk 
need be expected if the modulation of the transmitter is 
satisfactory, since guide 1 is always matched. It is not 
essential to produce an -section circulator in order to 
obtain the separation of m channels. A series of 4-section 
circulators can also be used. 


Ae(i) 


Fig. 3 


Fig. 4 shows an example of an arrangement which is of 
special interest when the channels are grouped in pairs. 


Fig. 4 


Isolators : 


An isolator is a junction of two guides satisfying the 
following condition: 


An isolator is represented by the diagram of Fig. 5. 


+ 
© © 


It necessarily contains a nonisotropic medium except 
where 7» is also zero and where no coupling exists be- 
tween guides 1 and 2 in either direction. 


Fig. 5 
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Necessary Presence of an Absorbing Body: It will now 
be shown that it must necessarily also contain an ab- 
sorbing medium. If the medium is assumed to be loss- 
less, (9) becomes in the general case of a junction of 2 
guides of any sort: 


Ti2 ful To 
To* T 12* 


wen 


From these equations, we deduce: 
R;-Ri* = Ry: Ro* 


Py: T1* = To1-T2* or fie = fai. (13) 


Eqs. (13) provide the following theorem: In lossless 
media, even when anisotropic, the transmission equiva- 
lent and the modulus of the reflection coefficient are 
independent of the direction of energy transfer. No con- 
dition exists for transit times which may be different in 
the case of anisotropic media. 

In particular, in the case of isolators, condition 
T>,=0 can only be verified if 7;,=0, that is to say if no 
coupling exists between the two guides. 

A lossless isolator cannot therefore be devised. This 
remark becomes generalized in the case of several guides. 
This theorem explains certain phenomena in guide junc- 
tions which appear to be paradoxical. 

Consider, for instance, Fig. 6, which is a guide carry- 
ing wave Ho and a ferrite insert of thickness e and 
length / placed in the zone of the guide where the mag- 
netic field is circularly polarized. It is well known? that, 
in this zone, the permeability tensor becomes a spherical 
tensor whose single component w+yp’ depends on the di- 
rection of propagation. The applied continuous mag- 
netic field and the ferrite insert may be so chosen that 
the effective magnetic permeability u, is equal to € in 
the direction! —2 and zero in the direction 2—1. 


©. | 
{ 

Vee eeeaaa ages ereseaeees | 

1 
a= 


Fig. 6 


In these conditions, if the length / is sufficient, we 
shall obtain the isolator conditions: 


T2 = 0 with T 12 7 it, 


However, these conditions can only be attained with a 
suitable thickness e. 

When thickness e tends to zero, no decoupling occurs 
between the alternating magnetic fields in the regions 


2J. R. Roven, “Ferrite in microwaves applications,” Bell Syst. 
Tech. Jour., vol. 32, p. 1333; November, 1953. 
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situated above the ferrite block and below it in the 
guide of Fig. 6. The coefficient 72 remains equal to 1. 
When the thickness e increases, these alternative mag- 
netic fields are soon decoupled and, in the direction 2 —1 
for which u.=0, these fields then tend towards the 
value corresponding to the guide limited by the ferrite 
insert and the lower metal part. 

When the frequency is sufficiently low, it can be seen 
that 7, then decreases very rapidly while Ty, remains 
equal to 1. But the general conditions of the foregoing 
theorem impose that the decrease of 7» shall necessarily 
be attended by an absorption in the ferrite block pro- 
portional to the thickness e. 

It is in fact this absorption which makes it possible to 
decouple the alternating magnetic fields in the two parts 
of the guide and, hence, to obtain a value of JT»; less than 
that of Zi. and which can even be zero. 

With certain ferrite elements, it is possible to obtain 
the conditions of an isolator (Tj,=1 and 7»=0) with 
the necessary losses remaining small. In that case, one 
can have | R,| =0 and | R,|1. 

With that isolator and a conducting wall farther into 
the guide, one can realize a cavity whose surtension is 
nearly independent of the frequency passband. That 
device constitutes an energy accumulator which can 
have several applications. For instance, in Fig. 7, in- 


short pulse 


(1/100 pts) 


long pulse (1 fis) 


4 


FERRITE 


MAGNETRON 


Fig. 7 


stead of the conducting wall, it is possible to put an 
impedance Z formed by an ionic or electronic commuta- 
tor; 1.e€., a device whose impedance is modified by the 
voltage V. For instance, we have 


VK<Vo-Z=0 
ns VoroZ = oO, 


When V is less than Vo, the cavity formed between 
the isolator and the impedance Z can be excited by a 
magnetron acting with relatively long pulses (1 us for 
instance) .® 

When the voltage V is greater than Vo, the cavity is 
open and one obtains, in the characteristic impedance 
Z, of the guide, pulses of large amplitude and very short 
time (107% us) when the isolator losses remain small. 

Use of Isolators: Isolators can be used in two different 
ways: 


3R. R. Warnecke and J. Ortusi, U. S. Patent No. 2,543,213. 
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Associated to an oscillator, isolators prevent fre- 
quency drift caused by waves reflected from the an- 
tenna. This advantage is specially valuable for a fre- 
quency modulated klystron. 

Similarly, in the case of its use on measuring benches 
where the isolation factor of the isolator (ratio Ti2/T») 
expressed in decibels, does not need to be particularly 
high. From the point of view of reaction on the trans- 
mitter an isolator of 20 decibels is equivalent to an iso- 
tropic attenuator of 10 decibels but has the advantage 
of almost fully maintaining the power of the generator. 
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Introduced in a delay line, isolators avoid, in traveling 
wave tubes, coupling between the input and output of 
the tube. Such coupling may arise either through a mis- 
match in the output guide,.or through the presence of a 
wave whose direction of propagation is opposite to that 
of the circulation of energy. The use of an attenuating 
impedance with isotropic medium, uniformly distributed 
or not, can suppress spurious oscillations only at the ex- 
pense of the tube efficiency and of a deterioration of the 
gain characteristic with frequency. The use of an isola- 
tor, comprising an anisotropic body in the delay line, 
avoids these disadvantages. 


Phenomenological Vector Model of Microwave Reflection 
from the Ocean* 
C. I. BEARD}, I. KATZ}, anp L. M. SPETNER{ 


Summary—A model of one-way transmission of microwave elec- 
tromagnetic signals over the ocean surface is developed from experi- 
ment. The received signal is described as a vector sum of a constant 
direct signal, a coherent reflected signal, whose amplitude and phase 
are fixed by geometry and sea state, and a fluctuating reflected com- 
ponent of random amplitude and phase. By interpreting experimental 
data in the light of this phenomenological model it has been possible 
to relate, quantitatively, the coherent and incoherent reflected signal 
and total signal to geometry and sea state. The results give support 
to the theoretical expression previously derived by Ament and others 
relating the coherent reflected signal to “apparent ocean roughness.” 
In addition, the general shape of the curve relating the incoherent 
scattering to “apparent ocean roughness” has been established and 
its asymptotic value found. 


INTRODUCTION 


N ONE-WAY transmission of microwave electro- 
if magnetic signals over the ocean between transmit- 

ting and receiving points within the line of sight, 
the received electric-field strength can be represented as 
the vector sum of a direct signal and a water-reflected 
signal. The resultant is a random, time-varying signal 
which is a function of geometry and ocean-surface pa- 
rameters, or “sea state.” Theoretical treatments of this 
forward scattering have been limited to special cases! or 
are formal solutions which have not been reduced to 
specific functional relations.2-4 Two recent reports 


* Manuscript received by the PGAP, August 11, 1955; revised 
“copy received, February 8, 1956. 

t+ Applied Physics Lab., Silver Spring, Md. 

1S. O. Rice, “Reflection of electromagnetic waves from slightly 
rough surfaces,” Comm. on Pure and Applied Math., vol. 4, pp. 351- 
378; August, 1951. 

2Lord Rayleigh, “On the dynamical theory of gratings,” Proc. 
Roy. Soc. A, vol. 79, p. 399; 1907. 

3 C, T. Tai, “Reflection and refraction of a plane electromagnetic 
wave at a periodical surface,” Harvard Tech. Report No. 28, Cruft 
Laboratory; January 15, 1948. ; 

4V. Twersky, “Multiple scattering of radiation by an arbitrary 
configuration of parallel cylinders,” J. Acoust. Soc. Am., vol. 24, pp. 
42-46; January, 1952. (More complete bibliographies are given in 
the papers of Rice and Twersky.) 


more closely related to the present work and which de- 
rive explicit relations for the coherent reflected term are 
those by Ament’ and Goldstein and Goldmuntz.® In 
addition, extensive unpublished theoretical work has 
been done at the Applied Physics Laboratory by Spetner. 

On the other hand, past experimental work has been 
regrettably deficient in that the ocean surface fluc- 
tuations have not in general been measured. It is only 
in the past few years that sea surface measurements 
have been included in propagation experiments to 
provide the data necessary for formulating an experi- 
mental description of forward scattering. The several 
experiments providing the data that made possible the 
analysis presented in this paper were performed by the 
Electrical Engineering Research Laboratory (EERL) of 
the University of Texas. EERL and the Applied Physics 
Laboratory of the Johns Hopkins University have been 


. associated since the beginning of 1951 in studies of 


microwave reflection from the water. 

It is the purpose of this paper to report a phenomeno- 
logical vector model of reflection which has been induced 
from these experimental data and which quantitatively 
relates the coherent and incoherent reflected signals and 
the total signal to geometry and sea state. This repre- 
sentation ties together many diverse experiments into a 
consistent pattern which has heretofore been lacking. 

This paper will describe the phenomenological model 
by first giving its vector representation and postulates, 
followed by the experimental supporting evidence for its 
validity. A homogeneous atmosphere will be assumed 
throughout. 


5W. S. Ament, “Toward a theory of reflection by a rough sur- 
face,” Proc. IRE, vol. 41, pp. 142-146; January, 1953. 

° H. Goldstein and L. A. Goldmuntz, “Mean square signal com- 
putation,” Nuclear Development Associates Report NDA-18-2; 
August 18, 1952. ah 
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VECTOR MODEL AND POSTULATES 


If electromagnetic radiation is transmitted over a 
calm ocean and the resulting field is probed vertically, 
a well-defined interference pattern will be found, indi- 
cating a strong coherent reflection from the surface. If 
this experiment is repeated as the ocean becomes 
rougher, the interference pattern will be less pronounced. 
With very rough surfaces there will be no identifiable 
interference pattern, indicating incoherent scattering is 
predominant. Attempts to describe this process analyti- 
cally led to the following assumption: 


Postulate 1 
The total signal (7) is the vector sum of (Fig. 1): 


1) the direct ray (D) of constant magnitude and 
phase, 

2) a coherent reflected ray (C) whose amplitude 
and phase are fixed by a given geometry and sea 
state, and 

3) a fluctuating or incoherent reflected component 
(IZ) of random amplitude and phase. 

(The quantities 7, D, C, and I represent electric 

field strengths.) 


Fig. 1—Vector model of reflection. 


If the vector J is assumed to be the resultant of indi- 
vidual vectors from a large number of independent ran- 
dom scatterers, following Goldstein’s treatment of back- 
scattering from the sea,’ then it follows: 


Postulate 2 


The incoherent vector J is the resultant of two orthog- 
onal independent vectors each representing a Gaus- 


7D. E. Kerr, “Propagation of Short Radio Waves,” McGraw- 
Hill Book Co., Inc., New York, N. Y., pp. 553-554; 1951. 
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sian random process with zero mean. (For conven- 
ience, the directions of the two orthogonal vectors are 
taken along and normal to C and their standard devi- 
ations and spectra are assumed to be equal.) 


DESCRIPTION OF EXPERIMENTS 


This paper draws on experiments performed by as- 
sociated laboratories in the course of investigations into 
propagation problems. These experiments were quite 
varied in nature using several frequencies in the micro- 
wave region; all involved only one-way propagation. 
Certain facts were gleaned from each of the experiments 
which, viewed together in the light of the model, form a 
consistent pattern. Since these investigations form the 
basis for the data a necessarily short résumé of each will 
now be presented. More details may be found in original 
reports. 


Lake Austin 


In February, 1952, the Electrical Engineering Re- 
search Laboratory (EERL) of the University of Texas 
performed an experimental study of the interference 
patterns existing over a relatively smooth body of water, 
Lake Austin, Texas.* The transmitter was placed near 
the water’s edge on one side of the lake while continuous 
measurements of field strength were made by a receiver 
mounted on a 50-foot elevator across the lake, a distance 
of 1,430 feet. (Such experiments are generally referred 
to as “height-gain measurements.”) Three wavelengths 
were used: 9, 3.2, and 0.86 centimeters. 


Gulf of Mexico 


In August, 1952, EERL carried out a series of propa- 
gation studies over two paths in the Gulf of Mexico to 
measure the reflection characteristics of ocean water at 
radio wavelengths of 9.0, 5.3, 3.2, and 0.86 centi- 
meters.°-1+, One path was across Barataria Pass at 
Grand Isle, La., a distance of about 2,800 feet, and the 
other was between two oil drilling platforms 5,000 feet 
apart, eight miles offshore. Height-gain measurements 
were made at Barataria Pass whereas fixed-point meas- 
urements over a time interval of a few minutes near 
interference maxima were made on the platform path. 
Some measurements were also made at 3.2 cm with dual 
antennas at the receiver end. 


San Clemente Island 


A continuous interference pattern was measured by 
EERL in October, 1953 using an airplane-mounted 
transmitter and a stationary 3-cm vertically-polarized 


8 A, W. Straiton, J. R. Gerhardt, A. H. LaGrone, and C. W. 
Tolbert, “Reflection of Centimeter and Millimeter Radio Waves from 
the Surface of a Small Lake,” EERL, Univ. of Texas, Report No. 63; 
May 15, 1952. 

9 A, H. LaGrone and C. W. Tolbert, “Reflection Studies of Milli- 
meter and Centimeter Radio Waves for Gulf of Mexico Paths,” 
EERL, Univ. of Texas, Report CM-753; October 31, 1952. 

10W. J. McKune and H. W. Smith, “Comparison of Power 
Spectrum Estimates of Overwater Microwave Radio Signal and As- 
sociated Water Waves,” EERL, Univ. of Texas, Report No. 68; May 
6, 1953. 

11 A, H. LaGrone, A. W. Straiton, and H. W. Smith, “Synthesis 
of Radio Signals on Overwater Paths,” EERL, Univ. of Texas, Re- 
port No. 71; April 30, 1954. 
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receiver 45 feet above the shoreline on San Clemente 
Island, California. The airplane flew toward the re- 
ceiver at a constant altitude of 125 feet starting at a 
range of 10,000 yards.” 


Golden Gate 


Further over-water propagation measurements were 
made by EERL in the spring of 1953 over the Golden 
Gate at San Francisco, California.!! The geometries of 
the two over-water paths used are summarized in 
Table I. Heights are above mean lower low water. 


TABLE I 
Receiver Transmitter 
Path Height Height ve 
(feet) (feet) 
I 186 85 16,000 
II 186 260 16,500 


The paths extended across a water surface approxi- 
mating open-ocean conditions. This experiment was 
quite varied in nature; it included measurements of total 
signal and signal with a null placed on the specular re- 
flection to obtain the direct signal magnitude for nor- 
malization; simultaneous reception was accomplished on 
two antennas separated vertically by 50 and 100 wave- 
lengths alternately. Extensive ocean-wave and meteoro- 
logical measurements were made." 


San Clemente, Convair 


An experiment similar to the one previously described 
under the EERL San Clemente Island Tests was car- 
ried out by Consolidated Vultee Aircraft Corporation 
at San Clemente Island, California, during May, 1954.1* 
The airplane runs were made at several constant al- 
titudes on incoming courses starting at ranges of about 
7,500 yards. A 3-cm transmitter using a wide antenna 
pattern was carried in the airplane; reception was ac- 
complished at a 200-foot-altitude stationary point at 
the shoreline. 


EXPERIMENTAL EVIDENCE SUPPORTING 
PHENOMENOLOGICAL MODEL 


From the experiments just described data were ob- 
tained which are now offered in evidence of the model’s 
validity. The initial clue was in the form of numerical 
values of the Rayleigh criterion of roughness. When it 
was clear that a transition zone existed between smooth 
and rough surfaces as measured by the parameter 
hy/X, a natural question was raised whether this pa- 


12 Since the pattern referred to above was obtained at APL’s re- 
quest only as an incidental sideline to other activity and was worked 
up at APL, no published reference is available. 

13K. H. Jehn, J. R. Gerhardt, D. F. Metcalf, and S. J. Prosser, 
“Some Meteorological and Oceanographic Characteristics of the 
Golden Gate, California, Area,” EERL, Univ. of Texas, Report CM- 
760; February 26, 1954. ; 

144 Records from these runs were made available to APL for re- 
duction and analysis; no published reference is as yet available for 
this experiment. 
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rameter was indeed a quantitative measure of surface 
roughness; the present section will show data supporting 
this. (The quantity h is here defined as the rms water- 
surface fluctuation measured from the mean water 
level; y is the grazing angle; and J is the electromagnetic 
wavelength.) From the form of the height-gain runs 
came a suggestion that incoherent noise increased with 
the same parameter; the experimental data led to the 
quantitative relation between J and hy/X. 


Transition Region ‘ 


Until recently the best estimate of the Rayleigh 
smooth-to-rough transition in the microwave region was 
given by Kerr as hW/X=25," where y is in milliradians. 
Evidence from height-gain and airplane runs at Lake 
Austin, Barataria Pass and San Clemente Island indi- 
cated transition values of: the roughness parameter 
hy / between 8 and 60. Although at first it would seem 
that this is a large spread it should be noted that the 
values bracket the critical value given by Kerr. Further- 
more, the random nature of the incoherent component 
makes the selection of a transition point quite difficult. 
Rather, the change from a smooth to a rough surface 
should be thought of as a zone of transition, not a point. 

Of the three experiments giving evidence on the tran- 
sition region, the Gulf of Mexico set was of most inter- 
est. In this set there were three different frequencies 
involved with both horizontal and vertical polarization. 
Yet a roughness parameter value of about 10+2.5 was 
selected for all runs. The variation was indeed small con- 
sidering the qualitative nature of what is called a tran- 
sition. 


Curve of Coherent Term 


The term “reflection coefficient” is generally taken to 
mean the ratio of coherent reflected ray to that of the 
direct ray, C/D. Since C/D decreases with increasing 
even for a smooth sea, in the following discussion this 
ratio will be divided by the smooth-sea reflection coeffi- 
cient r. One point on the curve of the coherent term vs 
effective roughness has already been established, !¢8 
namely, C/Dr=1 for hY/X=0. Values of the coherent 
term for hY/X>0 were obtained from the San Clemente 
(EERL) airplane run and the Golden Gate experiments. 

The continuous interference pattern obtained at San 
Clemente Island (EERL) was over a sufficient range of 
roughness to yield several values for the coherent term. 
They were obtained by drawing smooth curves through 
the maxima and minima of the pattern in the region of 
interest. Apparent reflection coefficients were computed 
from the upper and lower smoothed curves. These are 
plotted in Fig. 2. No dependable points were obtained 
for higher values of )/d because the smoothing process 


** This number was obtained by taking Kerr’s obstacle height as 
F9 =5.1 h. For this transformation see M. S. Longuet-Higgins, “On 
the statistical distribution of the height of sea waves,” Journal of 
Marine Research, vol. 11, p. 245; December, 1952. 

16 1). E. Kerr, op. cit., pp. 402, 423. 
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is no longer valid when the incoherent term becomes too 
large. The ordinates of the two points are estimated to 
have an accuracy of +20 per cent. Because of the un- 
certainty in wave height, the abscissa position is known 
to only +50 per cent. The solid curve is one previously 
derived from theoretical considerations by others °:* and 
has the analytic form 


C 
Dr OX) ae 8(rhp/n)?. 


In this equation y is in radians; everywhere else in this 
paper yw is in milliradians. 


©O=EERL April 1953 Golden Gate 
O= October 1953 San Clemente 


Vertical .Polarization 


fe en 8 (THY /A)> 


C/Dr 


(0) 100 
hy / d (mils) 


Fig. 2—Coherent term, theoretical and experimental. 


Additional points for higher hY/d on this graph were 
obtained from the Golden Gate experiment. These are 
the three points plotted between values of 100 and 250 
in b/d. These were obtained from the records of total 
signal at three antennas, recorded two at a time. 

The following analysis describes the method used for 
extracting the coherent term from simultaneous records 
of total signal without knowing one’s position in the 
interference pattern. According to the phenomenological 
model (Postulate 1) the total signal may be written 


T? = (D+ Ccos¢+1, cos ¢ — Iz sin ¢)? 
+ (C sing + I, sin ¢ + I, cos ¢)?, 


where ¢ is the angle the coherent vector makes with the 
direct ray vector and J, and J, are components of the 
incoherent vector along and normal to C. Assuming the 
mean values of J, and J, are zero (Postulate 2) and that 
¢ does not change appreciably from atmospheric fluctu- 
ations or water level changes during the 10-minute re- 
cording, T?=D?+C?+2DCcos ¢+/?. The phases at 
the three antennas, spaced vertically may be given as 
Ah; 


f= he = Qhe 
AR i 


gy = gqhy and $a = Gha 


where h; is transmitter height and ha, hy, and h, are re- 
ceiver antenna heights, respectively. The mean total 
signal at each antenna is given by: 
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Te? = D? + C? + 2DC cos (¢, + 2A) + TA 
Te = D?+ C2 + 2DC cos (6 + A) + Ty? 
T2 = D?+C?2+ 2DC cos¢. + T,2 


where 
A = g(hy — he) = q(ha — hs). 


The assumption is made that J? is identical at all three 
antennas, a reasonable assumption for the geometry in- 
volved in this experiment. Now 


T.? — T2 = 2DC[cos (¢. + 2A) — cos ¢] (1) 
T, — T2 = 2DC{[cos (¢. + A) — cos ¢.). (2) 


The value of the direct ray was measured in the 
Golden Gate experiment by using a binomial-array 
transmitter to place a broad null in the specular reflec- 
tion direction. The value of A is computed from geome- 
try. Thus there remain two unknowns C and ¢, which 
are now determinable from the two equations. This 
method of obtaining the coherent term is quite general 
and is independent of the statistics of the incoherent 
term. Its reliability is dependent on how accurately 
D, T? and the geometry are determined. 7,2 and T,2 
were independently assumed to be different from their 
measured values by up to +20 per cent (corresponding 
to the 10 per cent rms error found in D) and the result- 
ing errors in C were calculated. This represents the most 
pessimistic range of error since some correlation was ob- 
served between D, and D,. The greatest error in A is 
caused by the possible error of a transit sighting on the 
transmitter position. From surveying texts this was esti- 
mated as +1 minute and the rms assumed to be 3 min- 
ute. The error in C/Dr from this source is independent 
and was added to the first by the square root of the sum 
of the squares to form a resultant worst possible rms 
error of +100 per cent in the values of C/Dr. 


Incoherent Component 


The fact that, as W/d increases, the incoherent noise 
term J increases was based on examination of the data 
from the various experiments already mentioned. The 
point at hY/X=0 is immediately established by the ab- 
sence of “noise” on height-gain runs over a smooth lake 
surface. Before presenting further experimental data 
for hY/X>0, a derivation is necessary to indicate the 
method of obtaining J from total signal-strength meas- 
urements. 

First, for simplicity, let 6=0 so that C is in phase 
with D, and let D;}=D-+C. Then from Postulate 1: 


C= (Ort iy? ery : 
Assuming 
Te 


nee aol 
(Di + Iy)? 
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a series expansion can be made for 7. Expanding the 
individual terms in series and neglecting terms of higher 
order than J;,*, one obtains: 


r= (+1) += (1 By al ) 
are 2D, De a 
ioe 
 -8p3 


Making the following three assumptions of Postulate 2, 
that (a) I,=J,=0, (b) T2=I,, and (c) Jy and J, are 
Gaussian (then J,4= 304, where J,2=0°), one obtains: 


4 


2D/ 


Now T= Deg? under assumptions (a) and (b) 
above. Let 7,4 be the standard deviation of 7. Then, 


T= P= T Se 


Since, as will be shown later, ¢/Dr does not exceed 0.33, 
and since, at a maximum in the interference pattern 
D, cannot drop below D, the largest error in 7,@* arising 
from the approximation 7,40 is 3 per cent. 

In any experiment, therefore, the standard deviation 
of the total signal at or near an interference maximum 
is taken as the standard deviation of one of the two 
orthogonal components (J, or J,) of the incoherent noise 
vector I. 


0.6 


EERL 1952 Gulf of Mexico 


Vertical Polarization 
A tS a =i 
@ =X v=Q 


hw/d (mils) 


Fig. 3—Experimental values of incoherent term 
showing trend with hy/). 


Let us first examine data from the earlier Gulf experi- 
ments (Fig. 3) that led up to the final results, primarily 
to show the trend with ocean roughness. Fig. 3 gives the 
complete plot of the vertical-polarization data for all 
values of hw and X but plotted as a function of hW/d for 
the 1952 Gulf of Mexico experiments.!7 An estimate of 


17 Only four of the points on Fig. 3 are given in reference 9; the 
other points have been worked up from EERL’s original data. On 
August 7, 1952, the wave-gage accuracy was open to question and 
a best estimate from all data was made that h ~0.5 +0.25-feet rms. 
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the uncertainty bounds of the data, taking into account 
the wave uncertainty and the spread in the radio data, 
is given by the dashed lines. The solid line is a smooth 
curve drawn through the points. The trend shown in 
this curve is clear. For low values of hW/) the fluctuating 
component is near zero, and at high values of hy/X the 
curve levels off. Since the direct ray was not measured, 
however, no normalization to D is possible for these 
points. It may be noted, however, that normalizing to 
D and dividing by r would result in a steeper slope at 
low values of hw/d. Analysis of horizontal-polarization 
data obtained in the Gulf of Mexico® showed the same 
trend with #/d as shown in Fig. 3. 


Vertical Polarization 

©,®=EERL April 1953 Golden Gate ~ 
O=EERL August 1952 Gulf of Mexico (scaled 0.86cm) 
4= Convair May 1954 San Clemente 


o/Dr 


200 
hw/) (mils) 


(0) 100 300 400 


Fig. 4—Normalized values of incoherent term from experiment. 


The latest and most reliable results concerning inco- 
herent noise and for which normalization was possible 
are given in Fig. 4. At the lowest values of hy/\ reached 
in these data the value of o has already leveled off. One 
may only conclude that the transition region for 
“smooth” to “rough” for the incoherent energy lies be- 
low hy /\ =50. For the points denoted by open circles the 
direct-ray strength was measured by using a binomial- 
array transmitter which placed a null in the specular- 
reflection direction. For the other points, D was not 
measured separately, but at these high values of hW/n, 
C0, and T&D, so that these points are approximately 
normalized. The circled points represent data samples 
of ten minutes in length and the other points are from 
samples of variable lengths less than this. 

Another source of incoherent noise information is the 
measurement of the signal received from an airplane 
transmitter at large values of hy. The normalization 
approximation is used that T&D. In this case, however, 
T is varying as 1/R and hence in computing the stand- 
ard deviation each departure from the mean has to be 
divided by the mean at that point. In order to get 
enough points to compute Ta, it is preferable to use a 
very short section of each airplane run (to limit the 
variation in grazing angle) and to use many such runs 
to obtain a reliable value of 7,4. Four airplane flights 
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made by Convair off San Clemente Island on May 4, 
1954 were analyzed in this manner. The assumption in 
using these runs from a moving point is that the magni- 
tude of the incoherent noise is the same as in the static 
case and that only the spectrum may change. A sum- 
mary of the data is given in Table II below and the 
average is plotted on Fig. 4. The solid line is drawn at 
the level of the average of all the points. It is extended 
arbitrarily to zero by a broken line. The horizontal 
dashed lines indicate + twice the rms error.or the +20 
per cent bounds about the solid line. 


TABLE II 
Air- 
nue plane T sa hy / a T sa 
No. Height (mils) 
(feet) r Dr 
7 50 0.15 300 to 400 | 0.45 to 0.33 — 
11 50 0.13 300 to 400 | 0.45 to 0.33 a 
6 100 0.07 300 to 400 | 0.45 to 0.33 — 
9 150 0.10 300 to 400 | 0.45 to 0.33 == 
Ave = 0.11 +0.04| 300 to 400 0.39 OFSEEORL 


CONCLUSION 


By using a vector representation of the total, direct, 
coherent-reflected and incoherent-reflected signals it 
has been possible to arrive at quantitative relationships 
. between these signals and sea state: 


a) Experimental evidence is given in support of a 
theoretical expression previously derived by 
other workers of the relation of the coherent re- 
flected signal to “apparent ocean roughness.” 


b) The general shape of the curve of the incoherent- 
reflected term to “apparent ocean roughness” 
has been established and quantitative values 
have been assigned. 


Although this vector representation may be useful in 
predicting certain phenomena related to forward scat- 
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tering over smooth- and rough-water surfaces, it is 
realized that it is only a beginning. All the experiments 
mentioned in this paper were performed before the 
vector model was formulated and are only partially 
suitable. New experiments are needed to explore con- 
cepts suggested by the model. (Some experiments 
pointed in this direction were accordingly planned and 
were performed in April and May of 1955 by EERL 
under the sponsorship of APL.) For example, it is highly 
desirable to obtain information regarding the statistics 
of J. One method for obtaining such data is to measure 
rf phase as well as amplitude. In addition, in the plots of 
C and I vs b/d insufficient data are available to estab- 
lish any dependence upon water-wave direction. There 
is reason to believe that a wave-direction effect would 
exist since in the special case of short paths with the 
propagation direction perpendicular to the wave direc- 
tion as at Barataria Pass,®? long swells may lift the 
entire Fresnel zones up and down. This effect should de- 
crease as the radio-propagation direction approaches 
parallelism with the wave-propagation direction. Also 
for longer paths and short-crested waves the Fresnel 
zones are more likely to include a number of waves 
rather than just one, even when the radio- and water- 
wave-propagation directions are perpendicular. 

Likewise, all types of “seas” have been lumped to- 
gether in the parameter /, the rms water-surface fluctua- 
tion. For short ranges and long-period swells there may 
be only a few swells in the illuminated region. Since the 
statistics assumed for the incoherent term in Postulate 2 
depend upon the existence of at least 5 or 10 scatterers 
in the reflecting area, one would expect a departure from 
the simple model under these conditions. 

Refinements such as the effect of wave direction and 
the number of scatterers are suggested for future in- 
vestigations. What has been attempted in this paper is 
the development of the simplest general model that 
would tie together presently available experimental 
facts. 
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Report on Comparative 100 MC Measurements for Three 
Transmitting Antenna Heights” 
A. P. BARSIST anp R. E. McGAVINT 


Summary—This report evaluates measurements taken during 
August, 1952, at a frequency of 100 mc as transmitted from three 
transmitting sites at elevations ranging from 6,220 feet to 14,110 feet 
above mean sea level. Six receiving sites were used ranging in dis- 
tance from approximately 50 miles to 620 miles from the transmitter 
sites. Results are presented in terms of hourly medians of recorded 
field intensity and their distributions, as well as the over-all median 
values and deviations derived from these distributions. 


INTRODUCTION 
ee JULY and August of 1952 mobile re- 


cording equipment was set up near Anthony, 

Kan., and Fayetteville, Ark., in order to investi- 
gate radio wave propagation over the Cheyenne Moun- 
tain path during a typical summer month. This pro- 
cedure served to supplement the regular Cheyenne 
Mountain Field Station recording program! which pro- 
vides for continuous recordings at four receiving sites at 
frequencies ranging from 92 to 1,046 mc. The acquisition 
of a mobile 1 kw transmitter operating at a frequency of 
100 mc provided a possibility of comparing propagation 
from two 100 mc transmitting sources located at differ- 
ent ground elevations, and operating during alternate 
hourly periods. Accordingly, such alternate transmis- 
sions were made during the period August 6 to August 
13, using the mobile transmitter at Camp Carson which 
is lower than the regular fixed 100 mc installation, and 
during August 14 to August 20, using the mobile trans- 
mitter at a location on top of Pikes Peak which is much 
higher than the fixed installation. It is the purpose of 
this paper to evaluate the results of these measure- 
ments. 


TRANSMITTING AND RECEIVING STATIONS 
AND TRANSMISSION PATHS 


Tables I and II opposite give pertinent data on the 
antennas and paths used for these measurements. Spe- 
cial rhombic antennas were erected for transmitting at 
Camp Carson, and for receiving at Anthony and Fay- 
etteville, and a five-element horizontal Yagi array was 
built for the Pikes Peak site. The Pikes Peak installation 
including the transmitter truck is shown in Fig. 1. The 
90-degree corner reflector used on Cheyenne Mountain 
is fed by a folded dipole located 0.28 wavelength from 
the corner of the reflector. 

Pattern and gain measurements on the Camp Carson 
antenna were made by comparing the field produced by 
it with the field produced by a dipole mounted on the 


* Manuscript received by the PGAP, February 8, 1956. 

+ National Bureau of Standards, Boulder, Colo. 

1 “Tropospheric propagation research,” NBS Tech. News Buil., 
vol. 36, pp. 120-124; August, 1952. 
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rear pole of the rhombic, the measurements being made 
with mobile receiving equipment 15 to 20 miles distant. 

A similar procedure was used to determine the gains 
and patterns of the Anthony and Fayetteville receiving 
rhombics, except that the mobile unit was used for 
transmitting, and the signals were received on the 
rhombic, and on the auxiliary dipole mounted on the 
rear pole of the rhombic (see Fig. 2). In all three cases, 


DIPOLE RHOMBIC 


PATH OF 
TRANSMITTER 


qT 
[e} CHE YE, NNE s 
Ie 


TERMINATING 
RESISTOR 


FOR ANTHONY, KANSAS - A= | MILE 
FOR -FAYETTEVILLE, ARK. — A= V> MILE 


Fig. 2 


tests were also made on the dipole antenna with the 
rhombic antenna lowered, to ascertain the effect of the 
rhombic on the gain of the dipole, so as to obtain the 
corrections therefor. For the two receiving rhombics, 
these tests were made only along the main axis of the 
rhombic, and the correction so derived was assumed 
constant over the range of azimuth angles used in the 
experiments. Fig. 3 shows a diagram of the antenna 
matching section used for the receiving rhombics. 
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TABLE I 
PERTINENT Data—TRANSMITTING AND RECEIVING SITES 


a) Transmitting Sites 
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s Elevation 
Location ne Antenna Type of 
NI ove Mean Antenna 
N. Lat. W. Long. Sea Leyel 
Cheyenne Mt. 38.764° 104.864° 8,805 feet Corner 
Summit Reflector 
Camp Carson 38.696° 104.830° 6,258 feet © Rhombic 
Pikes Peak 38.839° 105.042° 14,115 feet Yagi 
b) Receiving Sites 
; Elevation 
Location of Antenna Type of 
5 _ Above Mean Antenna 
Nea cata Via leone. Saauieceal 
Kendrick 38.569° 103.984° 5,279 feet Dipole 
Karval RISA NO 5,079 feet Dipole 
Haswell Soo Oss, za 4,334 feet Dipole 
Garden City 37.833° 100.858° 2,874 feet Dipole 
Anthony | 37.240° 97.898° 1,375 feet |Rhombic 
Fayetteville 36.107°  94.107° 1,385 feet | Rhombic 


c) Receiving Sites—Distances from Transmitters (Miles) 


Cheyenne Mt. Camp Carson Pikes Peak 


Kendrick 49 .3 46.6 60.0 
Karval One 68.0 80.5 
Haswell 96.6 93.8 107.4 
Garden City 22625 223.6 DEI A 
Anthony 393.5 390.7 404.1 
Fayetteville 617.7 614.0 628.1 


In the case of the Yagi antenna, the pattern measure- 
ment was made at Camp Carson, by rotating the an- 
tenna and measuring the relative signal strengths re- 
ceived from Cheyenne Mountain. The gain was meas- 
ured at the maximum of the main lobe by using it as a 
transmitting antenna and comparing its field at 25 miles 
with that of a dipole set up in place of the Yagi. 


620 Ohms. 


50 Ohms 


Fig. 3—Impedance matching section, 100 mc rhombic antennas, 
Anthony and Fayetteville. 


All antenna gain data are given in decibels relative to 
an isotropic radiator under the assumption that all di- 
poles used either for comparative gain measurements, or 
as receiving antennas, compared closely to the ideal 
case as far as gain and pattern data were concerned. 
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TABLE II 
ANTENNA GAINS AND LINE LOSSES 
(ALL GAIN VALUES IN DECIBELS RELATIVE TO AN 
IsoTRoprc ANTENNA) 
a) Transmitting Antennas 
Toward Receiving Site 
Corner 
Reflector Rhombic Yagi 
Cheyenne Cp. Carson Pikes Peak 
Mt. Summit 
Kendrick 10.0 12.5 8.6 
Karval 9.9 as) 8.1 
Haswell 10.0 16.8 8.6 
Garden City 10.0 18.2 8.6 
Anthony 10.0 17.9 8.5 
Fayetteville 10.0 LEG 8.5 
b) Receiving Antennas other than Dipoles 
Toward Transmitting Site 
Rhombic Rhombic 
Anthony Fayetteville 
Cheyenne Mt. Summit 16.4 19.6 
Camp Carson 16.2 20.0 
Pikes Peak 16.4 19.4 


c) Line Losses Line Loss, decibels 


Transmitting Installations 
Cheyenne Mountain 
Camp Carson 
Pikes Peak 

Receiving Installations 
Kendrick 
Karval 
Haswell 
Garden City 
Anthony 
Fayetteville 


i) 
ee 


NwNoodo 


Dipole receiving antennas were used at Kendrick, Karval, Haswell, and Garden 
City. The theoretical gain value used for these antennas was 2.15 db relative to an 
isotropic radiator. 

Line losses were determined either by direct measure- 
ments, or by assumption of a loss figure per unit length 
in accordance with the manufacturer’s specifications. 


CALIBRATION AND DATA RECORDING 


Transmitter output power was determined by direct 
measurements with radio frequency wattmeters. Com- 
plete logs for all transmitting operations were kept. The 
transmitter output power was held closely to the nomi- 
nal values; any departure amounting to more than 
+0.1 db was taken into account in the data analysis. 
Operating frequency was kept within +200 cycles, any 
greater variations were corrected when detected. 

All receiver calibrations were made using signal gen- 
erators which were checked against a similar unit kept 
in the laboratory and used as a comparison standard 
only. At all times, suitable 50 ohm pads were inserted 
between the signal generators and the receivers. 

At Kendrick, Karval, and Haswell the output of the 
receivers was fed into Esterline-Angus graphic recording 
ammeters which were operated at a chart speed of three 
inches per hour. Calibrations were accomplished by con- 
necting the signal generator to the receiver input in 
place of the’antenna lead, and marking a series of signal 
levels on the charts. 
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At Garden City, Anthony, and Fayetteville, narrow 
band receivers were used in conjunction with Esterline- 
Angus recorders as well as with time totalizers similar to 
the apparatus described in detail by R. P. Decker.” The 
National Bureau of Standards time totalizer essentially 
contains 10 channels, each of which includes a de ampli- 
fier, and a bistable multivibrator actuating a fast-acting 
relay. If the receiver output voltage exceeds a preset 
level on any of the ten channels, the relay is actuated 
and a motor started which drives a revolution counter 
unvil the receiver output voltage drops again below the 
preset level. All ten channels are preset to the calibra- 
tion levels which are also marked on the recording 
charts, and the face of the counter panel is photographed 
automatically at hourly intervals by a relay-operated 
camera. Each film frame obtained by this method indi- 
cates also the date, the time of the day, and the preset 
calibration levels. 

At Anthony and Fayetteville the recording charts 
were operated at speeds ranging up to 1.5 inches per 
minute in order to obtain data suitable for analysis in 
case of failure of the totalizing or photographic equip- 
ment, and in order to provide records for subsequent in- 
vestigations of fading range and fading rate. 


PROCEDURE USED IN ANALYZING 
AND PRESENTING DATA 


In order to utilize the concept of transmission loss 
it was found convenient to calibrate the signal generator 
output in terms of decibels below 2 kw transmitter out- 
put. This permits direct tabulation of the recorded data 
in decibels of transmission loss, requiring only correc- 
tions for variations in the signal generator calibrations, 
for transmission line losses, and for transmitter power 
values different from 2 kw (applicable to the Camp Car- 
son and Pikes Peak Data). 

Where time totalizers were used, the differences be- 
tween hourly counter readings were expressed as per- 
centages of total times between readings, for each level 
and these were then graphed on probability paper and 
the median values obtained. Where time totalizers were 
not available, the Esterline-Angus charts were scaled 
for the same type of hourly data. 

After all hourly medians were determined, they were 
corrected as noted above, and the resulting hourly 
medians of transmission loss in db were classified into 
four groups as follows (for each of the six receiving 
sites): 


Cheyenne Mountain Transmissions Aug. 6-13 
Camp Carson Transmissions Aug. 6-13 
Cheyenne Mountain Transmissions Aug. 14-20 
Pikes Peak Transmissions Aug. 14-20 


2R. P. Decker, “Notes on the analysis of radio propagation 
data,” Proc. IRE, vol. 39, pp. 1382-1388; November, 1951. 
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For each of the final twenty-four groups, a cumulative 
distribution of all hourly medians was obtained, plotted 
on probability graph paper, and from each of these the 
over-all median and the interdecile range of hourly me- 
dians (as a measure of the standard deviation of hourly 
medians) was obtained (see Figs. 4 to 15, pp. 000-000). 

Table III shows the over-all median values, expressed 
as transmission loss in db, obtained for the four trans- 
mission periods at all receiving sites. In order to com- 
pensate for the difference in the number of recording 
hours available for each transmission path and the dif- 
ferent periods involved, the data were adjusted to the 
Cheyenne Mountain weighted average in the following 
manner, considering the fact that the received hourly 
median fields when expressed in db, are approximately 
normally distributed. 


TABLE III 


TRANSMISSION Loss VALUES 
OVER-ALL MEDIANS FOR 100 Mc MEASUREMENTS 
Avucust, 1952 


August 6-13 August 14-20 
Receiving 

Site Cheyenne Camp Cheyenne Pikes 
Mountain Carson Mountain Peak 

Kendrick 122.4 Ses WALES) 116.4 
Karval 123.9 151.8 12325 120.2 
Haswell 138.2 5250 138.3 126.5 
Garden City 185.0 180.6 185.5 182.0 
Anthony 177.8 7h as 181.2 178.6 
Fayetteville 194.7 191.3 195.0 193.0 


Above data adjusted to Cheyenne Mountain weighted average. 


Cheyenne Camp Pikes 

Mountain Carson Peak 
Kendrick 121.8 136.7 116.9 
Karval 123.8 USA 120.5 
Haswell 138.2 152.0 » 126.6 
Garden City 185.1 180.7 181.6 
Anthony 179.4 178.7 180.4 
Fayetteville 194.8 191.2 192.8 


For each receiving location the weighted average of 
the Cheyenne Mountain transmission loss was first de- 
termined by multiplying the over-all median for the 
August 6-13 period by the total number of hours con- 
tained in the distribution for that period, adding to that 
quantity a number obtained similarly for the August 
14-20 period, and dividing the total by the total num- 
ber of hours for both periods. This gives the weighted 
average of the Cheyenne Mountain transmissions. Sub- 
sequently the differences between this average and the 
over-all median for each of the two periods were deter- 
mined, and the Camp Carson and the Pikes Peak 
median corrected by algebraic addition of the difference 
for the appropriate period. The adjusted data are also 
shown in Table III. 
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Fig. 6—Site 2, Karval, August 6-13, 1952. 
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TRANSMISSION LOSS IN DECIBELS 
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Fig. 9—Site 3, Haswell, August 14-20, 1952. 
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Fig. 12—Site 5, Anthony, August 6-13, 1952. 
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Fig. 15—Site 6, Fayetteville, August 14-20, 1952. 
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The interdecile range of hourly medians in Table IV 
was adjusted in a similar way; but in this case the 
squares of the db values (a measure of the variances), 
were multiplied by appropriate number of hours instead 
of first-power values used in adjusting the medians. 


ABIES IV: 


INTERDECILE RANGE 
(DIFFERENCE IN DECIBELS OF MEDIANS EXCEEDED 
10% And 90% or THE ToTAL RECORDED Hours) 
Aucust, 1952 


August 6-13 August 14-20 

gp coovvang Site Cheyenne Camp Cheyenne Pikes 
Mountain Carson Mountain Peak 

Kendrick 5.4 8.6 3.4 ons 
Karval 7.0 1De3 5.9 4.1 
Haswell 141 Life Sloe) 5.8 
Garden City ies Zs 9.0 ORS 
Anthony WEG 7.8 iG) WALES 
Fayetteville 1S 10.7 a6) Teeth 


Above data adjusted to Cheyenne Mountain weighted average. 


Cheyenne Camp Pikes 

Mountain Carson Peak 
Kendrick 4.4 7.6 4.3 
Karval 6.6 11.9 4.8 
Haswell 1320) 14.9 ihe 
Garden City 155 1527) 170 
Anthony 13.8 9.0 13 
Fayetteville 9.2 8.4 eS 


It was felt that the number of recorded hours avail- 
able was not sufficient to permit a further breakdown of 
the results in order to determine diurnal variations. The 
total number of hours for which useful recorded data 
were obtained is shown in Table V. 


Ane} BI Ds 


NUMBER OF USEFUL RECORDING HOURS FOR 
COMPARISON RUNS 
Auecust, 1952 


August 6-13 August 14-20 
Receiving Site Cheyenne Camp Cheyenne Pikes 
Mountain Carson Mountain Peak 
Kendrick 53 53 Lae 65 
Karval (P i 40 33 
Haswell 83 87 UD: 65 
Garden City 67 71 25 25 
Anthony 75 70 63 54 
Fayetteville o2 52 47 42 


Finally, all transmission loss values for the adjusted 
medians were converted into decibels of attenuation 
relative to the free space field. This was done in accord- 
ance with methods described previously.? It should be 
noted that this conversion makes use of the antenna 
gains determined by assumptions or close-in measure- 
ments, and for the more distant receiving sites does not 
take into account the “loss-in-gain” concept which ap- 


3k. A. Norton, “Transmission loss in radio propagation,” PRoc. 
IRE, vol. 41, pp. 146-152; January, 1953. 
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pears as one of the conclusions of the Scattering Theory, 
as applied to long distance tropospheric propagation.* 

Table VI shows the adjusted results in db relative to 
the free space field compared to computed values for the 
sites within the radio horizon and in the diffraction zone 
(Kendrick, Karval, and Haswell). The computed values 
are based on the assumption of an effective radius of the 
earth corresponding to normal refraction, which is ap- 
plicable to the paths concerned in accordance with 
available meteorological data,> and the approximation 
of the actual path profiles by second degree curves fitted 
to appropriate portions of the paths.*7 


TABLE VI 


COMPUTED AND MEASURED FIELDS IN DECIBELS 
BELOW FREE SPACE VALUES 
AuGust, 1952 


Measured values adjusted to average of Cheyenne Mountain 
runs. 


Cheyenne Camp Pikes 
Mountain Carson Peak 
Com- Meas- Com- Meas- Com- Meas- 
puted ured puted ured puted ured 
Kendrick 7 Ome 23.5 31.5 41.4 Ost eeel SiR5: 
Karval OE DES 37.5 46.1 =O Ol Om 
Haswell 28.4 34.0 5725) *5479 210200 
Garden City 138 — dilie — 68.2 
Anthony oe Uthns — 84.4 — 76.5 
Fayetteville — 91.9 = 96.4 — 88.1 


The calculation of expected fields for receiving loca- 
tions far beyond the radio horizon is treated in a paper 
delivered at the 1953 IRE National Convention.‘ 

The adjusted median values from Table 6 are shown 
on Fig. 16, plotted vs distance. For greater clarity in 
presentation a logarithmic distance scale has been 
chosen. The lines connecting various points denoting re- 
sults for the receiving sites should not be construed as 
representing conditions at distances between the actual 
locations. Similarly, the line indicating the radio horizon 
is only approximate, and serves to illustrate the location 
of the sites relative to the radio horizon. 

In a similar way the adjusted interdecile range is 
plotted vs distance and shown in Fig. 17 illustrating the 
expected signal variations as a function of distance. 


DISCUSSION OF RESULTS 


Before attempting discussion or evaluation of these 
measurements, the following should be kept in mind: 


4J. W. Herbstreit, K. A. Norton, P. L. Rice, and G. E. Schafer, 
“Radio wave scattering in tropospheric propagation,” 1953 IRE 
CONVENTION RECORD, Pt. 2, Antennas and Communication, pp. 85-93. 

5B. R. Bean, “The geographical and height distribution of the 
gradient of refractive index,” Proc. IRE, vol. 41, pp. 549-550; April, 
1953. 

6K. A. Norton, “Transmission loss of space waves propagated 
over irregular terrain,” TRANS. IRE, vol. AP-3, pp. 152-166; 
August, 1952. 

7K. A. Norton, “Calculation of ground wave field intensity over 
a finitely conducting spherical earth,” Proc. IRE, vol. 29, pp. 623- 
639; December, 1941. 
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a) The calibration check of the signal generators 
could be made only after transportation over con- 
siderable portions of rough roads. It is unknown in 
what way or to what extent the signal generators 
may have been affected by transportation. 

b) Although considerable effort was made to match 
antennas and receiver inputs to the nominal 
50-ohm impedance, the stability of these values is 
not certain. Any mismatch in the antenna-receiver 
chain, or in the generator-receiver chain would atf- 
fect all calibrations. 

c) The experiments were conducted over a limited 
number of hours, and, due to some equipment 
failures, not all hours of transmission could be 
utilized at all receiving sites. The method of ad- 
justing all data to the weighted average of the 
Cheyenne Mountain transmissions is not neces- 
sarily ideal, but constitutes a reasonable approach 
to the correct values. 


Altogether it is felt that the relative values of the 
interdecile range could be determined more accurately 
than the absolute values of transmission loss. Apart from 
calibration difficulties, the latter are also subject to 
variation in antenna gains from the values assumed. It 
is very apparent that there is not much correlation be- 
tween computed within line-of-sight and diffraction 
fields, and the actual measurements. This may be chiefly 
attributed to the oversimplification one has to resort to 
in order to represent the rough terrain by a second- 
degree surface. The chief factor contributing to the 
computed field seems to be the character of the terrain 
in the vicinity of the receiving antenna, where, within 
the radio horizon, reflections are thought to take place. 
Studies are now under way to obtain sufficient data on 
transmissions within line-of-sight over various paths at 
various frequencies to correlate departures from com- 
puted values with relative terrain roughness and other 
features. 

However, the results obtained from the August meas- 
urements definitely confirm the existence of relatively 
high fields far beyond the radio horizon. Beyond a dis- 
tance of approximately 100 to 200 miles from the radio 
horizon the median field does not decrease at the high 
rate of the diffraction region, and the deviations from 
the medians tend to decrease from their maximum 
values at approximately 100 miles beyond the radio 
horizon (see Figs. 16 and 17). In interpreting Figs. 16 
and 17, one should also bear in mind that the lines con- 
necting the points which represent measurements at a 
particular site are drawn only to show the trend for the 
three transmitting heights compared. Especially a meas- 
uring location between Haswell and Garden City (at 
about 150 miles from the transmitters) would have been 
helpful in establishing the distance at which maximum 
long-term signal variations occur. 

The graph shows clearly that the influence of the an- 
tenna height on the over-all medians as well as on ex- 
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Fig. 16—Transmission loss relative to free space, 100 mc, 
August, 1952 (assuming antenna gains are realized). 
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Fig. 17—Variations of hourly medians, 100 mc, August, 1952, 
for three transmitting antennas. 


pected signal variations diminishes with distance from 
the transmitter. 


CONCLUSION AND RECOMMENDATIONS 


It is concluded that the measurements presented in — 
this report are useful for the study of long distance 
tropospheric propagation of very high frequency signals 
emitted from various types of antennas located at vari- 
ous elevations. They should be supplemented by: 


a) Identical experiments at other times of the year. 
b) Similar experiments at other frequencies. 
c) Study of short-term variations and the fading rate. 


It is also desirable to increase the reliability of the data 
by better methods of calibration of the equipment, and 
comparison with even more reliable standards of refer- 
ence. The procedure for antenna gain and pattern meas- 
urements may also be improved. It is expected that 
work now under way on application of Booker-Gordon 
scattering theory to propagation studies will offer a 
more realistic approach to effective gain figures of an- 
tennas operating on long path transmissions. 
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The Effect of Superrefractive Layers on 50—5,000 MC 
Nonoptical Fields* 
E. E: GOSSARD} anv L. J. ANDERSONt 


Summary—A body of radio data taken over an 80 nautical mile 
overwater link on four frequencies is analyzed, and the relationships 
to the vertical structure of the atmosphere determined statistically. 
The results are compared with the attenuation coefficients yielded 
by the theory for the first mode of a bilinear model as given by Furry. 

It is found that observation agrees reasonably well with theory 
at about 100 mc, but departs considerably for lower or higher fre- 
quencies. It is found that for a given height of the top of the super- 
refractive layer, an optimum frequency generally exists for which 
attenuation is a minimum. 


T HAS BEEN recognized for some time that no 
| existing method based solely on theory is adequate 

for predicting radio fields at all tropospheric fre- 
quencies in the diffraction region beyond the optical 
horizon. In order to predict such fields from meteoro- 
logical data, it is thus necessary to rely upon empirical 
methods. 

This paper describes a statistical investigation of the 
relationship of observed radio fields to certain critical 
parameters describing the meteorological profiles. The 
empirical results thus obtained are compared with the 
results of mode theory. An attempt has been made to 
keep the method as objective as possible and reduce 
assumptions to a minimum. 


THE DATA 


The radio data used were taken over an 82 nautical 
mile link between San Diego and San Pedro (Fig. 1). 
The transmitters and receivers were located 100 feet 
above sea level and the path was entirely over water.! 
Four frequencies were used: 52, 100, 547, and 5,000 mc. 
The fields were recorded in db above 1 pv and the re- 
spective free space fields were 64, 56, 38, and 68 db 
above 1 wv. Although detailed wiresonde data were 
available for some of the time during which the link was 
in operation, it was decided to use radiosonde data in 
this analysis since the results would then be directly 
applicable to other links for which only standard radio- 
sonde data were available. In the Pacific Standard 
Time zone radiosonde ascents are made starting at 
7 A.M. and at 7 P.M., hence the radio fields tabulated 
were averages of the hourly maximum fields recorded 
between 6 and 9 o’clock, both A.M. and P.M. 


ANALYSIS OF THE DATA 


All available San Diego radiosonde observations 
taken while the links were in operation were plotted in 


* Manuscript received by the PGAP, July 13, 1955. 

+ U. S. Navy Electronics Lab., San Diego 52, Calif. 

+t Smyth Research Associates, San Diego. 

1J. B. Smyth and L. G. Trolese, “Propagation of radio waves in 
the lower troposphere,” Proc. IRE, vol. 35, pp. 1198-1202; 
November, 1947. 
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Fig. 1—Map showing location of links. 


terms of a modified refractive index, B, where B is de- 
fined by 


B= (n — 1) X 10° + .012h. 


n is refractive index and h is height in feet. Those cases 
for which the refractive index profile was relatively sim- 
ple in form were then selected for the analysis. Specifi- 
cally, those cases were rejected which required more 
than a trilinear model to adequately represent the 
sounding. Thus, only one superrefractive layer—either 
surface based or elevated—was examined (Fig. 2). 

To determine the form of regression equation to be 
used in the analysis, scatter diagrams of all data for the 
various frequencies were plotted to determine the gen- 
eral form of the dependence of field strength on height 
and intensity of the layer. A previous study by Hoff- 
mann (unpublished) had shown that layer thickness Ah 
was not statistically significant. The dependence of the 
attenuation coefficient on the character of the layer ac- 
cording to mode theory was also taken into account in 
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TABLE I 
: | ; : 
Bred @lass*) |i ry2 r13 112.3 113.2 11(28) S12 $1.3 51.23 b c 
mcs 
Se E 104 0.522 —0.056 0.523 —0.075 0.521 6.10 Tol: 6.10 2520 — 3.66 
52 Susi 58 0.622 —0.035 0.685 —0.355 0.686 Sai) 742 5.41 D5 re/i — 8.1 
100 E eel 0.543 —0.445 0.644 —0.578 0.729 6.07 6.47 4.95 30.7 —26.8 
100 S.B. 72 0.743 —(0.419 0.742 —0.416 0.793 Ui 9.71 6.51 4122 —21.1 
547 E 114 0.416 —0.524 0.502 —0.593 0.679 10.9 10.3 8.77 40.1 —48.2 
547 93 0.210 —0.698 0.713 —0.858 0.865 14.1 10.3 7.24 45.9 —50.4 
5,000 ae 78 0.282 —0.543 0.289 —0.545 0.596 1525 {3h0 13.0 16.3 —45.2 


*E =Elevated cases, S.B.=Surface Based cases, T =Total cases. 


selecting the regression equation. The final equation 
used was 


F=a+bdlogAB+clogh 


where F is the radio field strength in db below the free 
space field. AB is the refractive index at the surface 
minus that at the height h (the top of the superrefrac- 
tive layer) (Fig. 2). 


Surface 


Fig. 2—Simplified trilinear and bilinear B curves 
illustrating h, AB and Ah. 


For the three lower frequencies sufficient data were 
available to permit the sample to be divided into two 
classes: 


a) superrefractive layer surface based. 
b) superrefractive layer elevated. 


RESULTS OF ANALYSIS 


The results of the statistical analysis are listed in 
Table I, above where: 


. 


N =number of observations in sample 

r12=total correlation between signal and log AB 

713 = total correlation between signal and log h 
719.3 = partial correlation between signal and log AB 
713.2 = partial correlation between signal and log h 
13) =multiple correlation of signal with log AB and 


log h 

S1.2=residual standard deviation (effect of log AB 
removed) 

51.3 =residual standard deviation (effect of log h re- 
moved) 


$1.03 =residual standard deviation (effect of log AB 
and log # removed) 
b=regression coefficient for log AB 
c=regression coefficient for log h 


A covariance analysis of the data indicated the differ- 
ence in the regression coefficient between the surface 
based and elevated layers was not statistically signifi- 
cant at the 0.05? significance level, although the differ- 
ence in the 100 mc “b” coefficients was almost signifi- 
cant. Therefore, for prediction purposes it is not worth- 
while to separate the two classes and the weighted 
average of the coefficients was determined. 

The multiple and partial correlation coefficients were 
in all cases significant at the 0.05 level, often at the 0.01 
level and sometimes at the 0.001 level. 

The plot of the regression coefficients vs frequency is 
shown in Fig. 3. 


DISCUSSION 


The fact that there is no significant difference be- 
tween the surface based and elevated cases makes it 
reasonable to compare the results of the analysis with 
those predicted by bilinear mode theory. In this model 
the critical meteorological parameters may be expressed 
in terms of h and AB, and the solution for the first mode 
is presented as a function of height (the height gain 
function) and a function of distance involving an at- 
tenuation coefficient. An effort was therefore made to 
obtain attenuation coefficients from the statistical anal- 


* If correlation is significant at the 0.05 level one chance in 20 
exists that the observed relation is due to chance. 

°W.H. Furry, “Theory of characteristic functions in problems of 
pine propagation,” M.I.T. Rad. Lab. Report 680; February, 
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ysis in order to compare them with those of theory. In 
order to keep the results as empirical and objective as 
possible, it was decided to approximate the effect of the 
height gain function by simply determining the field of 
the radio horizon by some standard method such as that 
outlined by Domb and Pryce.‘ The attenuation coeffi- 
cients then represent simply the linear rate at which the 
field (in db) decreases with distance beyond the horizon 
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Fig. 3—Plot of regression coefficients vs frequency 
for San Pedro-San Diego links. 


point (Fig. 4). The effect of superrefraction on the 
horizon field can be accounted for approximately by 
altering the effective radius of curvature of the earth, 
thus moving the horizon point farther from the trans- 
mitter and lowering the field. This method breaks down 
as trapping is approached, yielding negative attenuation 
coefficients. For large AB’s and thin layers, a further 
difficulty is encountered in assigning an effective earth’s 
radius. If we assume, however, that the height gain 
function corresponding to the first mode of the bilinear 
model is valid for such layers, we can solve for an effec- 
tive attenuation coefficient needed to yield the observed 
fields. Actually it developed that the attenuation co- 
efficients obtained by this approach were similar even 
for these extreme conditions to those obtained assuming 
a horizon field that was independent of refraction. This 
must be considered accidental, since the horizon field 
under large AB and thin layer conditions must depart 
considerably. from the standard case in a way that can- 


4C. Domb and M. H. L. Pryce, “The calculation of field strengths 
over a spherical earth,” JIEE, vol. 99, pp. 325-339; September. 
1947. 


Gossard and Anderson: Effect of Superrefractive Layers on Nonoptical Fields 


LP | 


not be taken into account by any simple procedure, but 
the similarity of this result with that obtained by using 
the height gain function adds confidence in the fixed 
horizon point technique—particularly when the super- 
refraction is not extreme. 
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Fig. 4—Curved earth representation of horizon ray, showing diffrac- 
tion region where attenuation coefficients are determined. Solid 
lines show field strength vs range for various attenuations. 


The empirical attenuation coefficients so obtained are 
compared with the theoretical in Figs. 5-8. 

Obviously, the regression analysis only yields values 
above the kvee in the attenuation curves, and only data 
above this critical # were examined. Furthermore, the 
noise level of the receiver was above the standard dif- 
fracted field, so the data on which this analysis is based 
were all taken under somewhat superrefractive condi- 
tions. To properly weight the near standard area of the 
graphs, the fields were assumed to go to the standard 
diffracted field as AB—0 and the curves in this region 
obtained by interpolation. 

Interesting features of the analysis are the following: 


1) In the present analysis the difference in the regres- 
sion coefficients between surface based and elevated 
cases was not significant. 


2) The departure of the observations from mode 
theory is such that fairly good agreement is noted at 
about 100 mc, becoming worse as frequency is either 
increased or decreased. 


3) The effect of elevating the layers does not diminish 
the field nearly as rapidly as would be anticipated by 
theory. 


4) For a given layer height, the variation of the re- 
gression coefficients with frequency as shown in Fig. 3 
leads to an optimum frequency for which the field is 
maximum (minimum attenuation coefficient). As the 
layer gets higher, this maximum field occurs at a lower 
and lower frequency. This is principally the result of the 
change in slope of the attenuation lines with increasing 
frequency. Physically it could be explained in terms of 
two opposing trends resulting from scattering within 
the superrefractive layer. Such scattering (perhaps by 
gravity waves) has been suggested by Anderson and 
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Fig. 5—Attenuation coefficients 52 mc. Dashed 
curves bilinear mode theory. 
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Fig. 6—Attenuation coefficients 100 mc. Dashed 
curves bilinear mode theory. 


Smyth.® Thus, as the frequency increases scattering be- 
comes more effective, but for a given height of the layer 
increasing the frequency will lessen the over-all effect of 
the layer. Therefore if the layer is very high the field 
would decrease with increasing frequency even though 
the scatterers are more effective. For h between 2,000 
and 3,000 feet the optimum frequency is about 400 mc. 


APPLICATION 


Although the purpose of this paper is to present the 
results of a statistical analysis of experimental data, it 
should be pointed out that the basic purpose in doing 
the work was to develop a method of predicting radio 
field strength beyond the horizon. The purpose in using 


5 L. J. Anderson and J. B. Smyth, “The effect of uniform layers 
on the propagation of radio waves,” TRANS. IRE, vol. AP-2, pp. 28- 
34; March, 1952. 
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Fig. 8—Attenuation coefficients 5,000 mc. Dashed 
curves bilinear mode theory. 


only radiosonde data for refractive index profiles and in 
presenting the results as attenuation coefficients beyond 
the horizon point, was aimed at the development of an 
objective method which could be applied to links of 
differing geometry as well as differing climatic condi- 
tions. Work is proceeding on the application of the 
method to such other links with the purpose of com- 
paring predicted fields with observation. The results of 
these comparisons will be presented in a subsequent 


paper. 
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Effect of the Ground Screen on the Field Radiated 
from a Monopole* 
J. R. WAITT 


radiation characteristics of an L.F. vertical an- 

tenna with a radial wire ground system. The in- 
vestigations! have been primarily concerned with the 
input impedance at the terminals of the antenna. This 
was shown to be mainly a function of the number and 
length of radials and the ground conductivity. It was 
assumed in most of this work that the ground wave field 
for a given current on the antenna was not appreciably 
affected by changes in size of the ground screen. Under 
this assumption, the radiation efficiency of the antenna 
is determined mainly by the input resistance. 

In an earlier paper,? an approximate method was 
given which was suitable for estimating the dependence 
of the ground wave field on the size of the ground 
screen. It is the purpose of this note to show that quanti- 
tative results can be obtained which support our earlier 
contention that the ground screen has only a small effect 
on the ground wave field intensity for a specified current 
on the antenna. 

The ground screen is assumed to be a perfectly con- 
ducting disc of radius a lying on a homogeneous flat 
ground of conductivity o. Choosing a cylindrical coordi- 


yee INTEREST has been shown recently in the 


* Manuscript received by the PGAP, December 15, 1955. 

+ National Bureau of Standards, Boulder, Colo. 

1F, R. Abbott, “Design of buried R.F. ground systems,” PRoc. 
IRE, vol. 40, pp. 846-851; July, 1952. 

2J. R. Wait and W. A. Pope, “The characteristics of a vertical 
antenna with a radial conductor ground system,” App. Sct. Research, 
vol. B4, pp. 177-195; March, 1954. 

3]. R. Wait and W. A. Pope, “Input resistance of L. F. unipole 
aerials,” Wireless Eng., vol. 32, pp. 131-138; May, 1955. 


nate system (p, ¢, 2) the antenna is considered to be 
coincident with the positive z axis and the surface of the 
ground is defined by z=0. Denoting H(p, z) as the mag- 
netic field of the antenna in the absence of any ground 
screen and AH,(p, z) as the change due to the presence 
of the ground screen, it follows from an earlier paper,? 
that 


AH 4(p, 2) Bice 
H 4(p, 2) no [o? + 2?]/? 


a ol} if 
i) eens 75) cos 4)p'dp’ (1) 
o'=0 H4”(p, 2) 
where 6 =27/\, \=free space wavelength, 7 = (uw/c)!? 
e'"/4 (surface impedance of the ground), ».=47X10-, 
w =angular frequency, 79 =1207 (intrinsic impedance of 
free space), 0=tan7! z/p, and where J, is the Bessel 
function of the first type. In the above, H,” refers to the 
field of the antenna over a flat perfectly conducting 
ground. This expression for the fractional change of the 
magnetic field is approximate and neglects terms of 
higher order in (n/n). It also assumes that the attenua- 
tion of the ground wave has a dependence with distance 
which is independent of the size of the ground screen. 
AH;,/H, can be regarded as the fractional change of the 
effective height of the antenna due to the presence of the 
ground screen. 

As an example a quarter-wave monopole antenna is 
considered with an assumed sinusoidal current distribu- 
tion. In this case 


180 
by I . 72 D4 
H3°(p’, 0) = e— iB lo” +0\/4)7] (2) 
2p’ 
and since p>; 
Hy(p, z) ~ =i e-eleteT cos (7/2 sin 0) (3) 
2p 


It then follows that the fractional change of the field is 
given by 


AH 4(p, Z) i — Bn 


H 4(p, 2) no € ; ) 
cos ee 


cos 6 


f e Ble"? +0014)"1" 7,(Bo! cos @)dp’. (4) 
0 


The integral, apparently, cannot be evaluated in closed 
form. It isnot too difficult, however, to evaluate it by a 
graphical method. This has been done for the case z=0 
which corresponds to the ground wave field. Letting 


AH 4(p, 0) 
H 4(p, 0) 


with 6= |n/no| , the integrals to consider are 


= 6|X1 + ix,| (5) 


a 
l 


27a! 
= f cos [ip + 93/4)? = w/4] (pap (6) 


2mra/Xr 
X.= f sin [(p* + 09/4)" — w/4] Julp)ap. (0) 


Some numerical values of Xi and X» are given in the 
following Table I. 


TABLE aI 

Qra/d Xy Xe 
0.0 0.000 0.000 
0.5 —0.042 0.040 
Le 0 —0.130 0.181 
15 (0.201 0.417 
2.0 —0.209 0.700 
25 —0.102 0.947 
3.0 0.042 1.093 
3H 0.155 Isley! 
4.0 On 7al Iss 
4.5 Onis 1.178 
5.0 0.050 1.300 
Soe 0.050 1.468 
6.0 0.119 1.612 
6.5 0.205 1.674 


The ratio of the field with the ground screen to the 
field H5(p, z) without the ground screen is then given by 


AH 4(p, Z) om 


Ae’® 8 
FR ue m 


Le 
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Fig. 1—The ratio of the ground wave field strengths of the antenna 
with and without the ground screen for a specified antenna current. 


where the amplitude ratio A is given by 


A = [(1 + 6X1)? + (6X2)? ]12 (9) 
and the phase factor is ® given by 
OX 2 
® = tan-! ———— - (10) 
1+6X. 


Since 6 has already been considered small, these are 
given adequately by 


A~wit+6xX, 
and 
@ ~ 6X_ (radians). 


The factor 6 can be obtained conveniently from the fol- 
lowing relation 


5 = 0.0075 (fme/o)1/2 


where fme is the frequency in mc and cq is the ground con- 
ductivity in mhos/meter. (Since displacement currents 
in the ground have been neglected throughout, the for- 
mulas are valid only when 6?<1.) Using (9), the ampli- 
tude ratio A is shown plotted in Fig. 1 as a function of 
2ma/X, the circumference of the ground screen in wave- 
lengths, for typical values of 6. 

It is interesting to note that for small screens the 
ground wave field strength is actually slightly less than 
it would be in the absence of the screen. As the screen 
becomes larger, the amplitude ratio increases somewhat, 
but is still only a few per cent greater than unity for a 
less than a wavelength. Thinking in terms of the recipro- 
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cal situation where the monopole is regarded as a re- 
ceiving antenna and the transmitter is located at some 
distant point on the surface of the ground, the increase 
of A above unity is characteristic of a “recovery” effect. 
Such a phenomenon occurs in ground wave propagation 
from land to sea.* 

Another interesting feature of the curves in Fig. 1 is 
that the minimum values of A occur approximately 
where the input resistance of the antenna is a minimum. 
It can be generally concluded, however, that the de- 
pendence of the L.F. ground wave field strength on the 


4G. Millington, “Ground wave propagation over an inhomo- 
geneous smooth earth,” Proc. IEE, vol. 96, pp. 53-63; January, 1949. 
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size of the ground screen is of minor significance com- 
pared to the dependence of the input resistance on the 
size of the screen. It should be mentioned that Page and 
Monteath® have used a similar method to calculate the 
radiation pattern of a vertical antenna over an irregular 
ground plane. 
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5H. Page and G. D. Monteath, “The vertical radiation patterns 
of medium-wave broadcasting aerials,” Proc. IEE, vol. 102, pp. 279- 
297; May, 1955. 


Admittance of Thin Antennas* 
GIORGIO BARZILAIt 


ITH REFERENCE to my recent paper,! I 
would now like to report some additional cal- 
culations. They refer to center driven half-wave 


antennas. 
The results are the following: 


60 
Vor — Vor = — (1.19 + 70.23) (1) 
CU Se Y 38 " 
: fi ae 
Vere ion od) 54 (3) 
Pe 23687: (4) 


Result (1) was obtained by introducing the current 
Iy(z) —Io(z) into the integral expression for the vector 
potential. The numerical integration was carried out by 
Simpson’s rule, after dividing the antenna into 16 equal 
parts. 

Result (2) was obtained using the method of moments, 
after dividing the antenna into 4 equal parts. To com- 
pute the moments, 


* Manuscript received by the PGAP, February 14, 1956. 

+ Instituto Superiore delle Poste e delle Telecomunicazioni, Fon- 
dazione U. Bordoni, Roma. 

1G. Barzilai, “On the imput conductance of thin antennas,” 
Trans. IRE, vol. AP-3, no. 1, pp. 29-32; January, 1955. 
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was expressed in terms of the sine and cosine integral 
functions. 
The final results arrived at are 


Q Gi Gs,9 Gow Gok Gs 
10 10.22 9.69 9.61 9.38 9.77 
15 10.22 9.85 9.81 9.60 9.94 
20 10°22 9.94 9.92 9.74 10.02 


It should be noted that in this case, since the current 
is prescribed at the input terminals, the resistance can 
be found immediately from the knowledge of the radi- 
ated power. For half-wave antennas, and in general for 
short antennas where the current can be prescribed at or 
near the input, both components of the admittance can 
therefore be calculated with good accuracy. From the 
values of the conductances recorded above, and from the 
values of the resistances, the following susceptances 
were found 


2 B, Ba Bow Box BS 

10 —5.94 —5.05 —4.90 —4.52 —5.68 
15 —5.94 —5.36 —5.29 —5.16 —5.78 
20 —5.94 —5.50 —5.46 —5.41 —5.82 


Admittances are expressed in millimhos. 
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